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NOTICE TO AUTHORS 


1. Communications.—Papers must be communicated to the Society by a Fellow. They 
should be accompanied by a summary at the beginning of the paper conveying briefly the 
content of the paper, and drawing attention to important new information and to the main 
conclusions. The summary should be intelligible in itself, without reference to the paper, 
to a reader with some knowledge of the subject; it should not normally exceed 200 words 
in length. Authors are requested to submit MSS. in duplicate. These should be 

typed using double spacing and leaving a margin of not less than one inch on the 
: left-hand side. Corrections to the MSS. should be made in the text and not in the 
margin. Unless a paper reaches the Secretaries more than seven days before a Council 
meeting it will not normally be considered at that meeting. By Council decision, MSS. of 
accepted papers are retained by the Society for one year after publication; unless their 
return is then requested by the author, they are destroyed. 


2. Presentation—Authors are allowed considerable latitude, but they are requested to 
follow the general style and arrangement of Monthly Notices. References to literature 
should be given in the standard form, including a date, for printing either as footnotes or in 
a numbered list at the end of the paper. Each reference should give the name and 
initials of the author cited, irrespective of the occurrence of the name in the text (some 
latitude being permissible, however, in the case of an author referring to his own work). 
The following examples indicate the style of reference appropriate for a paper and a book, 
respectively :— 

A. Corlin, Zs. f. Astrophys., 15, 239, 1938. 

H. Jeffreys, Theory of Probability, 2nd edn., section 5.45, p. 258, Oxford, 1948. 


3. Notation.—For technical astronomical terms, authors should conform closely to the 
recommendations of Commission 3 of the International Astronomical Union (Trans. 
1.A.U., Vol. VI, p. 345, 1938). Council has decided to adopt the I.A.U. 3-letter abbrevi- 
ations for constellations where contraction is desirable (Vol. IV, p. 221, 1932). In general 


matters, authors should follow the recommendations in Symbols, Signs and Abbreviations 
(London: Royal Society, 1951) except where these conflict with I.A.U. practice. 


4. Diagrams.—These should be designed to appear upright on the page, drawn 
about twice the size required in print and prepared for direct photographic 
reproduction except for the lettering, which should be inserted in pencil. 
Legends should be given in the manuscript indicating where in the text the 
figure should appear. Blocks are retained by the Society for 10 years; unless the author 
requires them before the end of this period they are then destroyed. 


5. Tables —These should be arranged so that they can be printed upright on 
the page. 


6. Proofs.—Costs of alterations exceeding 5 per cent of composition must be borne by 
the author. Fellows are warned that such costs have risen sharply in recent years, and it 
is in their own and the Society’s interests to seek the maximum conciseness and simplifi- 
cation of symbols and equations consistent with clarity. 


7. Revised Manuscripts.—When a are submitted in revised form it is especially 
requested that they be accompanied by the original MS. 


Reading of Papers at Meetings 
§. When submitting papers authors are requested to indicate whether they will be 
willing and able to read the paper at the next or some subsequent meeting, and approxi- 
mately how long they would like to be allotted for speaking. 


g. Postcards giving the programme of each meeting are issued some days before the 
meeting concerned. Fellows wishing to receive such cards whether for Ordinary 
Meetings or for the Geophysical Discussions or both should notify the Assistant Secretary. 
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MEETING OF 1953 NOVEMBER 13 
Dr J. Jackson, President, in the Chair 


The election by the Council of the following Fellows was duly confirmed :— 

John Grant Davies, M.A., Ph.D., 140 Manchester Road, Wilmslow, Cheshire 
(proposed by A. C. B. Lovell); 

C. Gordon Little, Jodrell Bank Experimental Station, Lower Withington, 
Macclesfield, Cheshire (proposed by A. C. B. Lovell); 

Marjorie Elsie Pillow, B.Sc., Ph.D., A.Inst.P., 16 Redford Avenue, 
Wallington, Surrey (proposed by R. W. B. Pearse); 

Harold Bytham Ridley, 20 Boileau Road, Barnes, London, S.W.13 (proposed 
by J. P. M. Prentice); 

William Marshall Sutherland, 39 Auchinloch Road, Lenzie, Lanarkshire 
(proposed by T. R. Tannahill); and 

Albert William Vince, Wyeross, 17 Stanley Hill Avenue, Amersham, 
Buckinghamshire (proposed by N. J. Goodman). 


The election by the Council of the following Junior Members was duly 
confirmed : 


Roger Derek Wade, 76 Stortford Road, Hoddesdon, Hertfordshire (proposed 
by R. M. Baum); and 

Derek Henry Wilson, M.Sc., 47 Parkfield Road, Ickenham, Uxbridge, 
Middlesex (proposed by F. Sheldrake). 


One hundred and two presents were announced as having been received since 
the last Meeting, including :-— 

E. Guyot, Reprints of astronomical articles (presented by the author); 

The Earl Nelson, Life and the Universe (presented by the author); 

F. A. Paneth, Reprints of articles relating to Thomas Wright (presented by the 
author); and 

M. F. Soonawala, Maharaja Sawai fai Singh II of Jaipur and his observatories 
(presented by Dr A. Armitage). 


45 
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Professor Subrahmanyan Chandrasekhar delivered the George Darwin Lecture 
(see p. 667), taking as his subject “‘ Problems of stability in hydrodynamics and 
hydromagnetics ”’. 


After the lecture, the President delivered a short address, and presented the 
Gold Medal of the Society to Professor Chandrasekhar, to whom it had been 
awarded for his contributions to mathematical astrophysics. The President 
said :— 

Professor Chandrasekhar, if I were to follow precedent it would now be my 
privilege and pleasure to present to you the Gold Medal of the Society. You 
have mentioned that for some purposes 10° years is not a long time, but you 
might have thought that a Society which has existed for more than a hundred 
years would have reached stability in presenting its Gold Medal. However, the 
Chancellor of the Exchequer has introduced a disturbing influence on our 
proceedings by the introduction of purchase tax. The Treasurer of the Society 
is here and so I must be careful what I say or do. Although I am not going to 
hand over the medal to you I do not think there is any objection to your seeing it 
now, while I assure you that in due course it will be forwarded to you. 

With the Gold Medal goes our best wishes for you in the future and on behalf 
of the Society I would like to express the hope that you may have the satisfaction 
of seeing the application of many of your theoretical results to the interpretation 
of the nature of the actual bodies of the Universe. 


MEETING OF 1953 DECEMBER 11 
Dr J. Jackson, President, in the Chair 


The election by the Council of the following Fellows was duly confirmed :— 


Sydney George Brewer, M.A., Fettes College, Edinburgh (proposed by 
D. W. Dewhirst); 

Joseph W. Chamberlain, Ph.D., Yerkes Observatory, Williams Bay, 
Wisconsin, U.S.A. (proposed by S. Chandrasekhar) ; 

George Herbert Avery Cole, B.Sc., Ph.D., 1a Alma Terrace, Allen Street, 
London, W.8 (proposed by H. S. W. Massey); 

Amulya Narayan Datta, M.Sc., Bengal Engineering College, Howrah, India 
(proposed by K. E. Bullen); 

John Edward Geake, B.Sc., M.Sc., A.Inst.P., The Physical Laboratories, 
The University, Manchester 13 (proposed by Z. Kopal); 

Johann Wolfgang Goethe, Dr.rer.Nat., 77 Falkstrasse, Frankfurt-am-Main, 
Germany (proposed by A. Beer); 

Richard Ingram, S.J., M.Sc., Seismological Observatory, Rathfarnham 
Castle, Co. Dublin, Eire (proposed by H. A. Briick); 

Alan Hamer Jarrett, B.Sc., M.Sc., Ph.D., A.Inst.P., The University Obser- 
vatory, St. Andrews, Fife (proposed by A. Beer); 

David Nelson Limber, Princeton University Observatory, Princeton, New 
Jersey, U.S.A. (proposed by S. Chandrasekhar); 

Robert J. Lockhart, B.A., University of Manitoba, Winnipeg, Manitoba, 
Canada (proposed by A. Vibert Douglas); 
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Hari Narain, M.Sc., D.Phil., The University, Sydney, Australia (proposed by 
K. E. Bullen); 

William Newsham, M.B.E., Brancepeth, 22 Sandbanks Road, Parkstone, 
Dorset (proposed by T. K. Metcalfe); 

James Ring, B.Sc., The Physical Laboratories, The University, Manchester 
13 (proposed by Z. Kopal); 

Allan R. Sandage, Ph.D., Mount Wilson and Palomar Observatories, 
Pasadena 4, California, U.S.A. (proposed by E. M. Burbidge); 

Norman D. Sykes, 56 Springfields, Welwyn Garden City, Hertfordshire 
(proposed by H. Wildey); 

Waltham Vaughan, 8 New Dock Road, Llanelly, Carmarthenshire (proposed 
by B. E. Featherstone); 

Joan George Erardus Gysbertus Voite, Daendelsstrasse 46, The Hague, 
Holland (proposed by W. S. van den Bos); 

William Leslie Wilcock, B.Sc., Ph.D., A.Inst.P., The Physical Laboratories, 
The University, Manchester 13 (proposed by Z. Kopal); and 

Robert I. Wolff, Ph.D., College of the City of New York, New York 31, 
U.S.A. (proposed by S. N. Milford). 


The election by the Council of the following Junior Members was duly 
confirmed :— 
Joyce Margaret Blackler, B.Sc., Royal Holloway College, Englefield Green, 
Surrey (proposed by W. H. McCrea); and 
Roderick Vernon Willstrop, Emmanuel College, Cambridge (proposed by 
B. G. Tunmore). 


Seventy-seven presents were announced as having been received since the 
last Meeting, including :— 


United States Hydrographic Office, Sight reduction tables for air navigation 
(presented by D. H. Sadler). 


Mr R. H. Garstang and Mr H. C. King were appointed Honorary Auditors 
of the Treasurer’s accounts for 1953. 








THE LAW OF RED-SHIFTS 


George Darwin Lecture*, delivered by Dr Edwin Hubble on 1953 May 


I propose to discuss the law of red-shifts—the correlation between 
distances of nebulae and displacements in their spectra. It is one of the two 
known characteristics of the sample of the universe that can be explored and 
it seems to concern the behaviour of the universe as a whole. For this reason 
it is important that the law be formulated as an empirical relation between 
observed data out to the limits of the greatest telescope. Then, as precision 
increases, the array of possible interpretations permitted by uncertainties in 
the observations will be correspondingly reduced. Ultimately, when a 
definite formulation has been achieved, free from systematic errors and with 
reasonably small probable errors, the number of competing interpretations will 
be reduced to a minimum. 

The path towards such a definitive formulation is now clear and the 
investigations are under way at Mount Wilson and Palomar. I shall discuss 
this programme in some detail after sketching the history of the law of 
red-shifts as a convenient background. 

There are three phases in the history, namely, the discovery phase which 
ended with a crude formulation in 1928-29, the rapid extension and improved 
formulation out to the limit of the 1o0-inch in 1929-36, and the current 
attempts to reach the limit of the 200-inch with the definite formulation, 
which began about two years ago. 


I. THe First Puase (Discovery) 
A. Radial Velocities 


The discovery emerged from a combination of radial velocities measured 
by Slipher at Flagstaff with distances derived at Mount Wilson. You will 
recall that the first velocity of a spiral nebula was measured by Slipher in 1912 
(—300 km sec for M 31). I mention the incident because the first step in a 
new field is the great step. Once it is taken, the way is clear and all may 
follow. In this case, however, Slipher worked almost alone, and ten years 
later, when he turned to other problems, he had contributed 42 out of the 
46 nebular velocities then available. The list was completely dominated by 
large positive velocities, ranging up to +1800 km/sec for NGC 584. 

Attempts were made, from 1916 onward, to derive the solar motion with 
respect to the nebulae; but the positive signs, indicating general recession, 
proved intractable. Then, in 1918, Wirtz attempted to “ save the phenomena ” 
by introducing a constant A term in the equations. The device improved the 
results somewhat but not enough to render them acceptable. Nevertheless, 

* Editorial Note-—This paper comprises the text of the George Darwin Lecture, which the 
late author had intended to revise before publication. His notes, together with the manuscript 
from which he spoke, made it clear what form he wished the published material to take. 
A reorganization of the original manuscript according to his marginal notes, with the addition 
of a few connecting sentences, was the extent of the editing required.—A. R. Sandage. 





No. 6, 1953 The law of red-shifts 


it did suggest that the next logical step would be to replace the constant K 
with rA, a term which varied with the distances of the individual nebulae. 
This procedure was suggested by the run of residuals in the earlier solutions 
but it was also encouraged by the rapidly developing studies of relativistic 
cosmology. 

Various solutions were made in the early and middle twenties, using 
apparent luminosities and apparent diameters as measures of relative distance. 
They failed but, as later appeared, only because the distance criteria were not 
sufficiently accurate to give reliable individual distances over the relatively 
short range covered by the radial velocities. The uncertainties are evident in 
the triple system comprising M 31 and its two companions, M 32 and 
NGC 205, where luminosities ranged from 100 to 1, and diameters from 60 to I. 


B. Distances from Brightest Stars 

The problem was finally solved when the brightest stars in late-type spirals 
and irregular nebulae were established as suitable distance indicators, and 
were calibrated by Cepheids in a few of the nearest nebulae. 

This new development started in 1923-24 when Cepheids were found in 
M 31 and M 33. The variables indicated the distances of the spirals, 
confirmed their suspected resolution and identified them as extragalactic stellar 
systems. Following up the clue, stars were soon recognized in many of the 
more conspicuous spirals, and the upper limits of brightness of these stars 
appeared to be substantially constant. It was this upper limit for the three 
or four brightest stars in great stellar systems which was used as a measure 


of distance. When, in 1928, it was applied to the nebulae with measured 
velocities, the law of red-shifts, or, as it was then called, the velocity-distance 
relation, emerged at once in approximately its present form. The law was 
established out as far as the Virgo Cluster, or about 7 million light-years, as 
an approximately linear relation, according to which velocities of recession 
increased at a rate of the order of 500 km/sec per million parsecs of distance. 


II. THE SecoND PHAsE (RECONNAISSANCE) 
A. Humason’s Adventures among the Clusters 

‘The second phase was entirely a Mount Wilson project, and lasted from 
1928 to 1936. Humason assembled spectra of nebulae and I attempted to 
estimate their distances. Attention was concentrated on clusters because their 
brightest members could be regarded as distance indicators, and, at the same 
time, offered maximum distances for a given apparent magnitude. Furthermore, 
the populations of great clusters are dominated by elliptical systems, which 
are the most highly concentrated of all nebulae, and hence the most easily 
observed with spectrographs. Faint elliptical nebulae give semi-stellar images 
and, therefore, they respond to the increased light-gathering power of great 
telescopes. 

Humason’s adventures were spectacular. He first observed some of Slipher’s 
nebulae, and then, when he was sure of his techniques, and confident of his 
results, he set forth. From cluster to cluster he marched with giant strides right 
out to the limit of the 1o0-inch. Spectrographs were improved from time to 
time—the names of Rayton and Schmidt and others are associated with the 
cameras——and each improvement was followed by a new surge forward. The 
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limit was bracketed when apparent velocities of 40 000 km/sec were recorded 
in the Bodtes Cluster and in Ursa Major No. 11, and readable spectra could 
not be obtained in the fainter Hydra Cluster (now known to have shifts half 
again as large). And always, throughout the march, the data fell in line with 
the relation as originally formulated. 

When clusters could not be observed to advantage, spectra of field nebulae, 
isolated or members of small groups, were assembled in large numbers. ‘These 
data furnished an account of red-shifts in the nearer, more conspicuous nebulae 
of all types. 

B. Results of the Second Phase 


When the second phase ended, the following results could be listed. 

(1) The law of red-shifts appeared to be linear (within the uncertainties of 
the distance estimates) out to approximately 250 million light-years, and the 
red-shifts, when expressed as velocities (c . d\/A), increased at the rate of about 
530 km/sec per million parsecs. 

(2) The solar motion is largely accounted for by galactic rotation, the 
residual representing the peculiar motion of the stellar system being small and 
uncertain (order of 100-150 km/sec towards galactic latitude about + 30° and 
longitude in the quadrant 50° to 140°). 

(3) The law of red-shifts does not operate within the Local Group. 

(4) Magnitude-velocity relations, fitted to data for clusters, resolved nebulae 
and field nebulae, and, of course, displaced from one another by magnitude 
intervals indicating relative absolute luminosities, could be used to calibrate 
all the relations in terms of brightest stars. Assuming M,, = — 6-35 for Brightest 
Stars, the corresponding values for Field Nebulae and Fifth Nebulae in 
Clusters were —15:20 and — 16-45, respectively. 


III. INTERIM BETWEEN SECOND AND THIRD PHASES 


A. Lists of Velocities Assembled by Humason and Mayall 


Towards the end of the second phase, Mayall entered the field and began 
assembling spectra with a very fine ultra-violet spectrograph on the 36-inch 
Crossley at Mount Hamilton. He paid special attention to the later-type 
spirals and other systems of low surface brightness for which his reflector was 
fully as efficient as the larger telescopes on Mount Wilson. From 1935 onward, 
the data increased in volume until now he has a list of the order of 280 nebular 
velocities as compared with Humason’s list of about 570. ‘There are some 
100 nebulae in common which ensure a reliable comparison of the two systems, 
and the two lists when published will include about 750 individual nebulae. 

These data furnish a wealth of information concerning nebulae of all kinds. 
Although full discussions must await publication of the lists, it is permitted to 
say that the data confirm and refine the results previously mentioned and, 
in particular, establish the order of precision in the measures (probable errors 
of the order of 30 km/sec or less) and of the dispersions among nebulae of 
different categories. Thus, by making reasonable allowance for the spread of 
absolute luminosities, the dispersion among isolated nebulae and small groups 
appears to be less than 150 km/sec. The dispersion among 11 members of the 
Local Group is about 60 km/sec, and, in general, the dispersions within groups 
and clusters increase steadily with compactness of the clusters. For instance, 
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a is 140 km sec for the small loose M 81 Group, 275 km sec for the Large 
Ursa Major Cloud, 700 km sec for the Virgo Cluster, 1050 km sec for the 
compact Coma Cluster, and perhaps 1200 km, sec for the Corona Borealis Cluster 
(from velocities of eight nebulae by Humason). 

These results are important for the present discussion because they indicate 
the reliability of the measured red-shifts for clusters used in the formulation of 
the law of red-shifts. The percentage errors are trivial except possibly for the 
nearby Virgo Cluster where it is about 5 per cent (in addition to the unknown 
peculiar motion of the cluster). With increasing distance the percentage errors 
diminish, and, beyond the Coma Cluster, they can safely be ignored. 
Evidently the uncertainties in the formulation of the law of red-shifts are 
entirely those in the distances. Furthermore, because distances are derived 
from the apparent faintness of objects of supposedly known luminosity, the 
difficulties seem to be largely those of photometry, except for uncertainties in 
the fundamental unit of distance. 


IV. THe Tuirp Puase (ATTEMPTS AT A DEFINITIVE FORMULATION 
WITH THE 200-INCH REFLECTOR) 


A. Extension of the Observed Range 


This was le situation when the completion of the 200-inch, and its 
accessories, initiated the third and current phase of the development. Because 
this phase represents an attempt at the definitive formulation of the law, I will 
be a little more precise in the use of terms. Red-shifts, for instance, will be 
used for the fractions, dA/A, by which details in the spectra are shifted. The 


shifts cannot be distinguished from Doppler shifts; they are constant throughout 
any given spectrum within the errors of measurement, the most reliable tests 
being those made on emission lines in large-scale spectra from H« to A 3727. 


The term “ apparent velocity ”’ will be discarded, and replaced by “ velocity ” 
signifying c.dd/A, or red-shifts expressed on a scale of velocities. The 
procedure is not formally correct but it is convenient. 

The current programme has two parts: one is the extension of the 
observed range and the other is the revision of the distance scale. The first 
part is straightforward, once suitable clusters are found. During the first 
season that a nebular spectrograph could be used on the 200-inch (1950-51), 
Humason obtained spectra in three clusters beyond the limit reached with the 
100-inch. They gave velocities of the order of 50000, 54000 and 
61 000 km;sec, respectively, the last being the Hydra Cluster which had been 
tried unsuccessfully with the 100-inch. Out to these limits, clusters can be 
selected from plates with the 48-inch Schmidt reflector, and the programme 
calls for groups of clusters at successive intervals to furnish mean points for 
the correlation curve. Beyond 60 000 km;sec, clusters cannot be identified 
with certainty on the survey plates. Many possible cases suggested by the 
48-inch are being checked with the larger reflectors, but otherwise we must 
depend upon purely chance finds with the 1o0-inch and 200-inch. Humason 
has now observed groups of three clusters each at 25 000 km,sec, at 40 000, at 
50 000 (representing both galactic hemispheres), and two clusters at about 
58 000 to group with Hydra. He is confident that larger velocities, at least 
up to one-quarter of the velocity of light, can be recorded when such clusters 
have been located. 
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B. Distances of Clusters 


The second part of the programme, namely, the revision of the distance 
scale, lags behind the measurement of red-shifts. There are several steps in 
the programme, namely, the revision of the unit of distance, the setting up 
of faint magnitude sequences in Selected Areas, the actual measurement of 
nebular magnitudes, corrections for effects of red-shift on apparent magnitudes, 
and, finally, the examination of the possibility of errors arising from 
inter-nebular obscuration and systematic evolutional changes in luminosities 
during the travel-time of the light from the nebulae. Let us consider these 
steps in their logical order, beginning with the revision of the unit of distance. 

1. The unit of distance.—A review of the unit of distance has been pending 
for more than twenty years—ever since it was realized that the many novae 
and globular clusters in M 31 were fainter with respect to the Cepheids than 
was the case in our own system. ‘The discrepancy was emphasized when Baade 
resolved the body of M 31 and recognized it as a Type II stellar population whose 
brightest stars were known to be comparable with those in globular clusters. 
Then, recently, Baade and his younger colleagues calibrated the cluster-type 
variables (and hence the brightest stars) in globular clusters by extending the 
colour-magnitude arrays in some clusters down to very faint limits on the main 
sequence, and comparing them with the corresponding arrays for stars in the 
vicinity of the Sun (stars so near that their distances were known from direct 
triangulation). Shapley’s value for the cluster-type variables was confirmed 
within small uncertainties.* 

The results clearly indicated that the classical Cepheids in M 31 are nearly 
1™-5 brighter than previously supposed, and, consequently, all our current 
estimates of absolute distances must be nearly doubled. The revision does not 
affect the estimates of relative distances, of course, but it nearly doubles the 
so-called “‘ age of the universe”. Further study will be required to reach a 
definitive value for the revision but meanwhile the factor 2 can be used in 
discussions. At the moment it seems likely that the final value may be less 
rather than greater than the even number. 

The revision was unexpected because Shapley made his primary calibration 
from parallactic motions of classical Cepheids in the galactic system (in 1918), and 
used it to derive his secondary calibration of cluster-type variables. It now appears 
that the latter is about correct while the former was in error by nearly 1™-5. 

2. Magnitude standards. The faint stellar magnitude standards and colours 
required for photometry near the limit of the 200-inch are being established 
in nine Selected Areas forming a convenient pattern around the sky. The 
first measures, reaching to about m,=18-5, were made by Stebbins and 
Whitford in three of the Areas, using their photo-cells on the Mount Wilson 
reflectors. ‘The programme has been continued and expanded by Baum, who 
has reached the 19th magnitude in all nine Areas and had pushed experimental 
measures in one or two of them down beyond the 21st magnitude. He is 
now putting in operation a new instrument designed to reach the limit of 
the 200-inch. The standards already available are sufficient for the faintest 

* Since this was written, preliminary results by Baum in the globular cluster M 13 suggest 
that grave uncertainties exist in this method of calibration. The colour-magnitude diagrams 
of the globular clusters M 3 and M 13 appear to differ by two magnitudes in the region of the 
main sequence. However, Blaauw’s direct calibration of the classical Cepheids and Savedoff’s 


fit of pulsation theory to the observations leave little doubt that the order of the correction to 
the classical Cepheids is 4M= —1™-5.—A.R.S. 
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clusters in which red-shifts have been measured, and Baum’s new programme 
should soon reach the faintest clusters in which red-shifts can be measured. 

3- Measurement of nebular magnitudes.—The measurement of _ total 
magnitudes of nebulae has been revolutionized since the war, both by the 
photo-cell in the hands of Stebbins and Whitford, Pettit, and others, and by the 
precise photographic techniques of Bigay and Holmberg. All these investiga- 
tions tend to reach general agreement on magnitudes (and colours, where 
measured) for given diameters. However, they all agree on the necessity of 
increasing apertures until no further change in luminosity can be detected. 
The formulation of a practical definition of diameter, however arbitrary, 
remains a problem for the future. The subject is under investigation but, 
meanwhile, the application of uniform principles to apparently similar nebulae 
in different clusters, using apertures adjusted to red-shifts, is expected to 
control the relative errors within reasonable limits (possibly of the order of 
o™-1 either way). 

The character of the energy-distribution curves for E nebulae shows that 
the corrections to observed magnitudes due to red-shifts are smallest in the 
wave-length region centred about 46300 A. With this circumstance, the 
definitive formulation of the law of red-shifts can be most accurately made in 
the red. Unfortunately, the most reliable method of measuring magnitudes— 
the photoelectric technique—does not yet reach the necessary limits in the red, 
although of course it does so easily in the blue and yellow. The compromise 
adopted at Mount Wilson and Palomar is to use a jiggle-camera in all three 
colours, blue, yellow and red, and to use the photo-cell for the blue and yellow. 
This procedure furnishes valuable information on colours as a function of 
distance, and also permits the photo-cell calibrations in blue and yellow to be 
applied to the red as well. One additional safeguard has been added to the 
jiggle-camera programme. Because these images must be large compared with 
the in-focus images, and the size of throw must be adjusted accordingly, it is 
found convenient to make images in all three colours with a uniform, medium 
size, and to repeat in one colour only (generally yellow) with double the throw. 
Corrections for size of in-focus images can then be transferred to the other 
colours on the plausible assumption that colour does not vary much throughout 
the images of elliptical nebulae. 

A jiggle-camera for the prime focus of the 200-inch, giving smooth, square 
images with throws of 0-5, I, 2 and 4mm on a side, was put into operation 
last autumn by my colleague, A. R. Sandage. Plates have now been assembled 
on half a dozen of the fainter clusters, although the delay in getting transfers 
of stellar magnitudes from Selected Areas, due to exceptionally poor winter 
weather, has seriously hampered the progress in final reductions. Baum 
proposes to join in the programme, now that his new photo-cell is ready for 
operation, and he has already contributed some data to the calibration problem. 
Meanwhile, Pettit has measured many of the nebulae in the nearer clusters with 
a photo-cell on the Mount Wilson reflectors in the course of his programme 
of getting magnitudes and colours of all nebulae for which radial velocities have 
been measured (except the very faint cluster nebulae observed with the 200-inch). 

These remarks summarize the state of the programme. They indicate why 
the present discussion must necessarily be a progress report, with perhaps 
some indications of the direction in which any revisions of the earlier results 
are likely to fall. 
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4. Effects of red-shifts on apparent magnitudes.—The total etfects of red-shifts 
on apparent luminosities depend upon their interpretation. The mere fact 
that we observe the shifts clearly indicates that each light-quantum from the 
nebulae reaches us with reduced energy. Regardless of the interpretation of 
red-shifts, we must accept the loss of energy by the individual quanta (or reject 
the fundamental relation E .A=constant) and we must correct the apparent 
magnitudes accordingly. These corrections are routine and have the same 
kind of validity as those for atmospheric extinction or instrumental 
characteristics. The bolometric luminosities are evidently reduced by the 
factors (I+dA/A), but since the red-shifts change the wave-length scales of the 
energy-distribution curves, the effects are selective. If the true, unshifted 
energy curves are known, they can be distorted by the observed shifts, and the 
effects then traced through the atmosphere, the telescopes, and the measuring 
devices to give the desired corrections. This procedure has been possible only 
since Stebbins and Whitford in 1947 derived the energy curve for the elliptical 
nebula, M 32, reduced to no atmosphere. The curve is derived from 6-colour 
photo-cell measures, ranging from A 10300 to A 3530. 

Because elliptical nebulae, which dominate cluster populations, are 
remarkably homogeneous as to spectral types, colours and distribution of 
colour throughout their images, the energy curve for M 32 can be applied 
directly to the estimation of corrections for red-shifts to apparent magnitudes 
of all elliptical nebulae in any region of the observed range in the spectrum. 
Minor uncertainties may arise from the facts that M 32 is a dwarfish nebula, 
that it is superposed on an outer arm of M 31 and, finally, that it is in the low 
galactic latitude of —22°. The curve will be redetermined eventually from 
giant elliptical nebulae in the Virgo Cluster, but the present formulation seems 
to be fairly reliable and it is all we have. 

Luminosities of nebulae are measured through “ windows ”’ consisting of 
the atmosphere, telescope and plate-filter or other receiving devices. By 
means of the energy-curve for M 32, the energy getting through a window from 
a nebula with a known red-shift can be compared with the energy that would 
have got through in the absence of a red-shift. These ratios furnish the desired 
corrections required by the mere observed presence of the red-shifts. These 
corrections I once called the “‘ energy effects’ and, I repeat, they must be 
applied to the measured luminosities, regardless of the interpretation of 
red-shifts. 

5. Progress report.—In eleven clusters blue and yellow magnitudes are 
available for a number of the brightest nebulae, and in two of the faint clusters 
red magnitudes also have been measured. Among these data, corrected for 
galactic obscuration and for red-shifts, the Ist, 3rd, 5th and roth nebula in 
each cluster have been used as distance indicators, and the magnitudes plotted 
against logarithms of red-shifts. The curves suggest that the brightest nebulae 
in the different clusters are about equally luminous, but that the spread between 
Ist and 10th varies with the richness of the cluster or with the size of the 
sample. On this interpretation of the data, the most accurate results should 
be furnished by the roth nebula in the great clusters only. This procedure 
represents a compromise, because a higher ordinal number, say the 25th nebula, 
would be expected to reduce accidental errors but at the price of introducing 
systematic errors due to variations in size of samples. 
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The yellow magnitudes offer the most reliable results, because corrections 
for red-shifts in the blue are large and uncertain while the red measures are 
not yet assembled in sufficient numbers. The favoured correlation, at the 
moment, is 

logyg U=0°2 m,, +116 (10th nebula), 


where v=c.daja. 


This correlation is shown in Fig. 1. Evidently the correlation is linear within 
the uncertainties of the data, and the residuals are surprisingly small. There 
are no indications of absorption in space in these particular directions. As the 
number of observed clusters increases, the residuals should offer a critical test 
of the presence of inter-nebular absorption. Meanwhile the colours offer an 
approximate criterion where spectra have not yet been recorded. 
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Fic. 1.—The relation between velocity and apparent magnitude. The new data obtained 
with the 200-inch are represented by the last four points on the regression line. Humason’s 
red-shifts are expressed on a scale of velocities as c . dij A in km/sec. The photovisual magnitudes 
have been corrected for the energy effect only. They do not include the recession factor. 

Similar correlations are readily derived for other distance indicators in 
clusters, for isolated field nebulae, for isolated pairs, and for brightest stars 
{blue supergiants) in late-type spirals and irregulars with Type I populations. 
The curves are all parallel, and are displaced by 0-2 times the difference in the 
absolute magnitudes of the distance indicators. In this way it is possible to 
calibrate all the correlations in terms of any one, the most accessible, of course, 
being the brightest stars in nearby nebulae. Thus, the M,,, of the roth nebula 
seems to be of the order of —17-4 on the old distance scale, and —18-9 or 
somewhat less on the new. 

Another method is offered by novae and globular clusters in Type II 
populations found in the Virgo Cluster. The rather meagre results now in 
hand appear to be consistent with those from brightest stars in Type I 
populations. When precise photometric data are available, the two methods 
together should furnish a reliable calibration for all the correlations. 

Meanwhile, the relations between red-shift and apparent magnitude 
represent relative distances that should be reliable, and the slope of the curve, 
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logyy v =0-2m+C, indicates a linear relation to a very close approximation. 
It should be emphasized that the magnitudes have not been corrected for 
recession of the nebulae. 

Such corrections are closely represented by the term, Am=dd/A, which 
must be subtracted from the measured magnitudes if the red-shifts are 
interpreted as Doppler shifts. The corrections range up to a maximum of 0™-2 
for the Hydra Cluster, and, if the uncorrected data define a linear relation, the 
corrected data would lead to a non-linear relation in the sense of an accelerated 
expansion. However, the data now available are not sufficient to furnish a 
critical test of the two interpretations; either one is permitted by the uncertainties, 
although it seems unlikely that a decelerated expansion can be represented. 

In conclusion it should be repeated that the discussion necessarily ends 
with a progress report, and not with a definitive solution of the problems 
considered. Nevertheless, it can be stated with some confidence that a definitive 
formulation of the law of red-shifts will be available before long, in the form 
of a relation between red-shifts and apparent magnitudes, out to red-shifts of 
the order of 0-25. 

When no recession factors are included, the law will represent approximately 
a linear relation between red-shifts and distance. When recession factors are 
included, the distance relation is expected to be linear or non-linear in the 
sense of accelerated expansion. The “age of the universe”’ is likely to be 
between 3000 and 4000 million years, and thus comparable with the age of 
rock in the crust of the Earth. 

The implications of the law established over this distance may be traced 
out to nearly double that distance by its effect on the apparent distribution of 
nebulae in depth, N,,, =f (m). 

Thus, if red-shifts do measure the expansion of the universe, we may be 
able to gather reliable information over a quarter of its history since expansion 
began, and some information over nearly a half of the history. 

As for the future, it is possible to penetrate still deeper into space—to 
follow the red-shifts still farther back in time—but we are already in the region 
of diminishing returns; instruments will be increasingly expensive, and progress 
increasingly slow. The most promising programmes for the immediate future 
accept the observable region as presently defined, hope for only modest extensions 
in space, but concentrate on increased precision and reliability in the recorded 
description. ‘The reconnaissance is being followed by an accurate survey; the 
explorations are pushed towards the next decimal place instead of the next 
cipher. ‘This procedure promises to reduce the array of possible worlds as 
surely as did the early rapid inspections of the new territory. And later perhaps, 
in a happier generation, when the cost of a battleship can safely be diverted 
from insurance of survival to the consolations of philosophy, the march outward 
may be resumed. 

For I can end as I began. From our home on the Earth, we look out into 
the distances and strive to imagine the sort of world into which we are born. 
Today we have reached far out into space. Our immediate neighbourhood 
we know rather intimately. But with increasing distance our knowledge fades, 
and fades rapidly, until at the last dim horizon we search among ghostly errors 
of observations for landmarks that are scarcely more substantial. The search 
will continue. The urge is older than history. It is not satisfied and it will 
not be suppressed. 





PROBLEMS OF STABILITY IN HYDRODYNAMICS AND 
HYDROMAGNETICS 


George Darwin Lecture, delivered by Professor S. Chandrasekhar on 
1953 November 13 


I 


Most of what I shall say may be considered as properly belonging to applied 
mathematics ; indeed, in large part they will be based on the works of the great 
applied mathematicians: Lord Rayleigh of the past and Sir Geoffrey ‘Taylor 
and Sir Harold Jeffreys of the present. And applied mathematics though it be, 
I hope I can succeed in bringing some measure of conviction to you that the 
considerations I shall present have some bearing on and some relevance to the 
problems of astronomy. 

Let me begin then by stating in general terms what we may expect to learn 
from studies on stability. 

When we know that an object has existed in nearly the same state for a long 
time we generally infer that it is stable; and by this we mean that there is 
something in its construction and in its constitution which enables it to withstand 
small perturbations to which any system in Nature must be subject. More 
precisely, what we have in mind when we say that an object is stable is this: if the 
object experiences an external perturbation which produces small displacements 
among its component parts, then it reacts in such a way that it eventually restores 
itself to its original state. In practice this definition of the meaning of stability 
requires some amplification: for, though an object may be strictly unstable, it 
may still be that the time which an external perturbation needs to manifest 
itself is so long compared to the total duration of its existence that its potential 
instability has not had time to reveal itself. All that we can conclude, then, 
from the knowledge that an object has endured for a certain period of time is that 
its time of instability is long compared to that period. If this is the case, what, 
one may ask, can we learn from an investigation of the stability of an object? 
It appears that an investigation of the stability of an object can, under favourable 
conditions, provide information of two sorts. First, we can often infer some 
significant facts concerning the construction and the constitution of an object 
merely from the knowledge that it has successfully endured the effects of small 
perturbations for a known period of time; I shall presently give illustrations of 
this. Second, an investigation of stability may disclose the relative importance 
and sometimes the conflicting tendencies of the different forces and constraints 
which are kept dormant in the undisturbed state of the object. Thus, when we 
are confronted with a novel object—and most astronomical objects are novel— 
a study of its stability may provide a basis for a first comprehension. And I may 
confess that for an applied mathematician problems of stability present a 
particular attraction: by their very nature, these problems lead to linear equations 
and linear equations are always more pleasant to deal with than nonlinear ones. 
While to emphasize problems of stability for this last reason may appear as an 
undue concession to one’s limitations, I think it will be admitted that an attempt 
to unravel the wide range of astronomical phenomena in which turbulence, 
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rotation and magnetism all play their parts, is not likely to meet with success 
if one started with the most general equations of hydrodynamics and electro- 
magnetism without discrimination. In any event, I should like to illustrate by 
means of examples as to what we might learn about some of the newer problems 
of astronomy from studies of stability. 


Il 


I said that we can often infer some significant facts concerning the constitution 
of an object from the mere knowledge that it has existed for a long period of time. 
Let me illustrate this by a classical example from the theory of stellar structure 
and then by a more modern version of the same. 

Consider a sphere of perfect gas—a‘“‘ star’’ let us say—in equilibrium under 
its own gravitation. Then a simple calculation shows that its total energy, €, 
including its internal heat energy (@) and its gravitational potential energy ({2), 
is given by 
=f JQ|, (1) 

3(y—1) 
where y denotes the ratio of the specific heats. Now it is clearly necessary 
that € be negative; otherwise the star will expand to infinity releasing energy 
in a time measurable in days. The condition €<o requires that y>}. This 
conclusion regarding a physical parameter describing the star is clearly a 
significant piece of information. 

The meaning of the condition y > becomes clearer when we examine the 
stability of the gas sphere to its own natural modes of oscillation. By considering 
adiabatic radial oscillations, Ledoux (1) showed by a very simple argument 
based on the virial theorem that the formula 


o? =(3y—4)| Q|/J, 


€= 


-M 
where I= | r? dM(r) (2) 
0 


denotes the moment of inertia, gives o the frequency of the oscillation to a 
sufficient degree of accuracy under most conditions. The meaning of the 
condition y > 4 is apparent from this formula: if this condition is not met, the 
frequency of the oscillation will become imaginary and the star will be 
dynamically unstable, i.e. unstable with respect to one of its own natural modes 
of oscillation. All this, of course, is well known. But let us apply similar 
considerations to the magnetic variables discovered by Horace Babcock. 

It can be shown that for a gas sphere of the dimensions of a star in which 
there is a prevalent magnetic field, H, the total energy, €, including the heat (@), 
the gravitational (£2) and the magnetic (#) energies is given by (cf. Chandrasekhar 
and Fermi (2)) 


8 as 
€= — (12 |-M), (3) 


where A= = ||| | Pax, diy dy (4) 


From equation (3) it follows that one of the conditions for equilibrium is 


fi <|Q|. (5) 
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This condition clearly sets an upper limit to the strength of a magnetic field 
which may prevail in a star. Since for a spherical configuration of uniform 
density 2= —3GM?/5R (where G denotes the constant of gravitation, M the 
mass and R the radius of the configuration), we can estimate the limit on the 
magnetic field set by (5) by using this expression for Q in it. In this manner 
we find 

V/|H [Pay.<2-0 x 10° M/R? gauss, (6) 
where M and R are now expressed in solar units. While for most of the magnetic 
variables discovered by Babcock the surface fields measured are a hundred to a 
thousand times smaller than the limit for the root mean square field set by (6), 
there are a few stars—three to be exact—for which the surface fields already 
approach the limit set by (6) to within a factor ten. 

Now the meaning of this upper limit to the magnetic field which may prevail 
in a star can be understood by considering the period of radial adiabatic 
pulsation. By following the method of Ledoux, we can show that in this case 
equation (2) is replaced by (cf. Chandrasekhar and Limber (3)) 


o* = (3y—4)(|2|-M)/T; (7) 
the meaning of the condition #l< |Q| is apparent. From formula (7) it would 
appear that we can make o?-0 (i.e. make the period tend to infinity) by letting 
fl—-|02|. It is possible that this is the explanation of the known relatively long 
periods of the magnetic variables. I am here referring to the fact that while 
the observed periods of the magnetic variables are of the order of 6 to 10 days, 
the period of radial pulsation which we should predict for these stars, were they 
normal, is of the order of $ day. If the suggested explanation of the long 
periods is correct, then the peculiar A-stars for which surface fields of the order 
of 104 gauss have been measured must have internal fields which are a thousand 
times larger; and this is perhaps not impossible. 


III 


I should like to consider next a problem in gravitational stability which 
may have a bearing on galactic structure (cf. Chandrasekhar and Fermi (2)). 
Consider an infinite cylinder of an incompressible fluid of uniform circular 
cross-section of radius R,. Is this gravitationally stable? We can answer this 
question by considering a perturbation of the cylinder which deforms the 
boundary into 
r=R+acos kz (8) 


(where a<R, z is measured along the axis of the cylinder and & is a constant) 
and asking whether the resulting change in the gravitational potential 
energy, AQ, per unit length is positive or negative. ‘The fact that the potential 
energy per unit length of an infinite cylinder is infinite requires that some care 
be exercized in the evaluation of AQ. 

When AQ is evaluated it is found that 


AQ>o for kR,>1-067 (9) 


and AQ <o for kR,<1-067. 


The cylinder is, therefore, gravitationally unstable for all symmetrical deforma- 
tions of the boundary with wave-lengths exceeding Ax =27R,/1-07. This recalls 
to mind a well-known result of Rayleigh’s (4) that a liquid jet is unstable, on 
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account of surface tension, for all symmetrical deformations with wave-lengths 
exceeding the circumference of the cylinder. And as in Rayleigh’s discussion 
of the stability of the liquid jet, we may ask whether there is a mode of maximum 
instability. We find that there is one whose wave-length is approximately 
1°84 times the minimum wave-length, Ay, for instability; and the time of 
instability as judged by this most unstable mode is 


7 = 4:07(47Gp) 2. (10) 


If we substitute for p in this formula the value of the mean density of the 
interstellar matter in the spiral arm in which we are located (namely 
2x 10-*4 g/cm’) we obtain +=10% years. This is so short a time that a 
structure like the spiral arm, if we can idealize it as an infinite cylinder, would 
have disintegrated a long time ago. 

Now the spiral arms of a galaxy show every sign of instability ; but we should 
be happier if the time of instability were of the order of 5 x 10° years. It appears 
that we can make an infinite cylinder gravitationally stable for periods of this 
order by supposing that along the axis of the cylinder there is a magnetic field 
of the order of 

H,=47pRyy G. (11) 


Thus a field H=H, increases the wave-length of maximum instability by a 
factor ten and the time of instability by a factor nine. ‘Table | gives more data 
relating to this problem. From this table we may conclude that if the spiral 
arms can be idealized in the manner we have done, the requirement that they 
be unstable but with times of instability of the order of 5 x 10° years will imply 
that there are magnetic fields of the order of 7 x 10°® gauss along the arms. 
It is gratifying that a number of other independent evidences lead one to 
postulate galactic magnetic fields of this same order. 


TABLE I 
Wave-lengths A, and X,, at which instability sets in and at which it is maximum, and 

characteristic time, 7, needed for instability to manifest itself for case R,=250 parsecs and 
p=2x107*4 g/cm? 

H Ay Nes T 

(gauss) (parsecs) (parsecs) (years) 

° “5 X 10° < 10° 10° 
10° 2X30" 
10° 1°8 x 10° 
10? : 10° 
10 7x 10° 
107 10° 
10° ' ro** 
10° ' 101! 


5x 20~* 10° 
io* 3°3 X 10° 
io 6 10° 
iq ’ 10! 
< 1078 5°2 x10 
<10 
ro * 10° 


Ne Uw dN 


mw 
= po 


IV 


In a general way one may expect that a conducting fluid in the presence of 
a magnetic field and under the influence of Coriolis’ force is subject to conflicting 
tendencies and that it may exhibit unexpected patterns of behaviour. I should 
like to illustrate this by considering a classical problem in the theory of thermal 
instability. ‘The problem I have in mind is the instability of a layer of fluid 
heated below which has been investigated by Rayleigh (5), Jeffreys (6) and others. 
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When a horizontal layer of an incompressible fluid is heated below— 
“underside” as Rayleigh calls it—then on account of thermal expansion the 
liquid above, being colder, will be heavier. This is a top-heavy arrangement 
and is potentially unstable. Consequently, if the temperature gradient is made 
sufficiently adverse, we should expect the instability to manifest itself. 

‘The manner of the onset of convection under these circumstances has been 
the subject of both experimental and theoretical investigations since Bénard’s 
first experiments on the subject in 1900 and 1g01. Those early experiments of 
Bénard established that instability sets in only when the temperature gradient 
exceeds a certain critical value and that when it does set in, it does so as a pattern 
of cellular convection. The correct interpretation of Bénard’s experiments 
was given by Rayleigh in 1916. Rayleigh showed that what decides the stability 
of a layer of fluid heated below is the numerical value of the non-dimensional 
quantity 


R= 2 a, (12) 


where g denotes the value of gravity, d the depth of the layer, 8=|d7 dz], the 
constant adverse temperature gradient which is maintained, and x, « and v are 
the coefficients of volume expansion, thermometric conductivity and kinematic 
viscosity respectively. We shall call R the Rayleigh number. Both experiments 
and theory show that instability must set in when the Rayleigh number exceeds 
a certain determinate critical value: for example, it is 1708 when it is confined 
between two rigid plane boundaries. Accordingly, higher temperature gradients 
can be maintained, before instability sets in, in a liquid of higher viscosity and or 
higher thermal conductivity ; and the reason for this is quite evident. 

Before passing on to anexamination of the effects of Coriolis’ force and magnetic 
field on the onset of convection by thermal instability we may first consider very 
briefly the manner in which one determines the critical Rayleigh number, 
theoretically. One determines this by supposing that an initial constant state 
is slightly disturbed and asking whether there exists any type of disturbance 
which will not be damped. Since any disturbance in the horizontal plane can 
be expressed as a stiperposition of two-dimensional periodic waves, we can 
investigate this problem by asking what the lowest Rayleigh number, R(a), is 
at which a mode of disturbance in the horizontal plane of a given wave number, a, 
when excited, does not get damped. On solving this problem one finds that the 
resulting function R(a) has in general a single minimum which it attains for a 
particular value of a. It is clear that this minimum of the function R(a) specifies 
the critical Rayleigh number at which instability will first set in; at the same 
time the value of a at which the minimum occurs determines the dimensions of 
the cell which will appear at marginal stability. 

We shall now suppose that in addition to gravity the fluid is subject to the 
Coriolis acceleration 2u x Q resulting from the fluid partaking in a rotation with 
angular velocity 2. ‘The discussion of the stability of fluid heated below under 
these circumstances can be carried out and it is found (Chandrasekhar (7)) that 
the effect of the Coriolis acceleration is to inhibit the onset of convection, the 
extent of the inhibition depending on the value of the non-dimensional parameter 

® cost 
4 tt d', (13) 
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where # denotes the angle which the direction of 22 makes with the vertical and 
Q=|Q|. The inhibition factors are of the order of ten for T of the order of 108 
(see Fig. 1). It is further found that 


R->constant T?% and a—constant T!® as T> «~; (14) 


here a denotes the wave number (in units of 1/d) of the cells which appear at 
marginal stability. I might mention here that experiments carried out by 
Dr D. Fultz at Chicago confirm the predicted dependence of R on 7. 

In a general way the reason for the inhibiting effect of 2 is clear: according 
to a theorem of Kelvin and Helmholtz, the vortex lines have a tendency to be 
dragged with the fluid, the attachment of the fluid to the vortex lines being 
stronger the lower the viscosity and the higher the angular velocity; and in the 
limit of zero viscosity the attachment is a permanent one. Consequently as 
© cos # increases and/or v decreases, motions at right angles to the vertical will 
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Fic. 1.—The variation of the critical Rayleigh number, Re, for the onset of convection as a 
function of T(=4Q2d‘ cos? B/v*) for the three cases: (a) both bounding surfaces free, (b) both 
bounding surfaces rigid and (c) one bounding surface free and the other rigid. 


[Reproduced by courtesy of the Royal Society (Proc. Roy. Soc. A, 217, 306, 1953)] 


become increasingly difficult ; and this prevents an easy closing-in of the stream 
lines required for convection. This also explains why instability, when it sets in, 
does so for a value of a which is increasingly large: for, a large value of a means 
that the cells are elongated in the direction of the vertical and motions in this 
direction are not hindered by the component of Q in this direction. 

According to equations (13) and (14), in the limit 7—> ©, the relation giving 
the critical temperature gradient at which convection sets in changes from 


g«B,=constant xv d-4 (15) 
to 


ga«8,=constant «(024 cost #/d4v)"3, (16) 
This changed dependence of 8, on d from a d~*-law to a d *°-law is likely to be 
a decisive factor in determining the character of the convection when Coriolis’ 
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force is acting. If the indications of this theory may be taken as a guide as to 
what may happen in reality, then we should conclude that the effects of Coriolis’ 
force must control thermally induced convection in all atmospheric layers 
exceeding ten metres on the Earth and a hundred kilometres on the Sun. 

We turn next to the effect of a magnetic field on the onset of convection. 
On general grounds, we may expect that the effect of a magnetic field will also 
be to inhibit convection and that this inhibiting effect will be greater the stronger 
the magnetic field (H) and the higher the electrical conductivity (o): for, when 
the field is strong (or the conductivity high), the lines of force tend to be glued to 
the material and this will make motions at right angles to the vertical increasingly 
difficult. Moreover, when cellular convection does set in, we should expect that 
the cells are elongated in the direction of the vertical, the elongation being 
greater the stronger the vertical component of the magnetic field. In the limit 
of infinite electrical conductivity, when the lines of force are permanently 
attached to the fluid, convection in the usual sense will become impossible. 
A detailed theoretical treatment of the problem (cf. Chandrasekhar (8)) confirms 
these anticipations. In particular, it is found that the critical Rayleigh number 
for the onset of instability in the presence of a magnetic field depends on the 
strength of the field and the electrical conductivity o through the non-dimensional 
parameter 

2 2 
on an _ 
pv 
where # now denotes the inclination of the direction of the impressed magnetic 
field to the vertical. It is further found that we have inhibition factors of the 
order of ten for O of the order of a thousand (see Fig. 2) and that in the limit O—- x, 


R->constant Q and a—constant Q'°®. (18) 


All that I have said so far regarding the effect of 2 and H on the onset of 
thermal instability might possibly have been predicted on general grounds. 
But let us now consider what happens when both Q and H are simultaneously 
present. The results I shall describe have been obtained for the case when 
g, H and Q are all co-planar and both the confining boundaries are free surfaces ; 
however, the general character of the solution to be described does not depend 
on these latter r strictions. 

The solution of the problem presents some unexpected features which 
result from the curve R(a) (which represents the lowest Rayleigh number at 
which a horizontal disturbance having a wave number a when excited is 
undamped) having two minima, the one occurring at the larger a giving the lower 
Rayleigh number when Q is small and the higher Rayleigh number when Q is 
large (see Figs. 3 and 4). Thus for 7\(=7/7*)=100000 and Q,(=Q/z*)=50 
the two minima occur at a (the wave number measured in the unit 1/d)=18-6 
and 3-48 where R= 4-03 x 10° and 7-26 x 10°, respectively. On the other hand 
for O, = 100 (and the same value of 7) the two minima occur at a, = 18-2 and 3-37 
where R = 3-98 x 10° and 3°93 x 10°, respectively. Consequently, for a value of QO, 
slightly less than 100 the wave number of the cells which appear at marginal 
stability will suddenly decrease from a= 18-2 to a=3-4. In other words, if we 
start with an initial situation in which 7, has the value ro®° and no magnetic field 
is present and gradually increase the strength of the magnetic field, then at first 
the cells which appear at marginal stability will be elongated; but when the 

46* 
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magnetic field has increased to a value corresponding to Q,= 100, cells of two 
very different sizes will appear simultaneously: one set which will be highly 
elongated and another set which will be much less elongated. As the magnetic 
field increases beyond this value, the critical Rayleigh number will actually 
begin to decrease. However, for sufficiently large Q the inhibition due to the 
magnetic field will predominate and will take control of the situation. at his is an 
unexpected sequence of events and I do not know if one could have predicted it. 





2-0 ! 
1-0 2-0 3-0 4-0 





Fic. 2.— The variation of the critical Rayleigh number, Re, for the onset of instability as @ function 
of O for the three cases (i) both bounding surfaces free (curve marked 1), (ii) both bounding surfaces 
rigid (curve marked 11) and (iii) one bounding surface free and the other rigid (curve marked II). 
The points marked 1, I and III on the R,-axis are the limiting values for the three cases considered 
in the absence of a field. 

[Reproduced by courtesy of the Editors of the Philosophical Magazine (Phil. Mag. (7), 43, 501, 
1952)] 


V 


I shall turn now from these problems in thermal instability to one of a 
somewhat different kind: the problem of the gravitational instability of an infinite 
homogeneous medium. ‘This was first considered by Jeans (g) more than fifty 
years ago. In this problem we start from an initial state of homogeneity and 
rest and consider the velocity of propagation of a density fluctuation through the 
medium. If the gravitational effects of the density fluctuation are ignored, the 
problem is the classical one of the propagation of sound ; and as is well known, the 
velocity of sound in a gaseous medium is independent of wave-length and is 
given by 

c= \(yp/p), (19) 
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Fic. 3.—The dependence of the wave number a of the cells which appear at marginal stability 
under the simultaneous influence of a Coriolis acceleration and a magnetic field. The differeiit 
curves refer to the various assizned values of Ty and varying values of Q;. 
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Fic. 4.—The variation of the critical Rayleigh number R for the onset of instability as « 
function of Q, for various assigned values of T,. 
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where p denotes the pressure, p the density and y the ratio of the specific heats. 
On the other hand, if the change in the gravitational potential, 5V, consequent 
on the density fluctuation, 3p, is taken into account in the equation of motion, 
with 5V related to 5p through Poisson’s equation 


V6 V = — 47Gép, (20) 
then, as Jeans showed, the velocity of wave propagation is no longer independent 
of the wave-length. It is given by 

Vj? =c?—47Gp k?, (21) 
where &( = 277A) is the wave number. The velocity of wave propagation, therefore, 
becomes imaginary if 


aG mc 
pe < EP og yt TO =),? (say). (22) 


( 
c* Gp 


But an imaginary velocity of wave propagation only means that the amplitude of 
the wave will increase exponentially with time. Accordingly, if an arbitrary 
initial perturbation in density is represented by a Fourier integral, then the 
amplitudes of those components in the Fourier representation which have 
wave-lengths exceeding A, will increase exponentially with time; A, is therefore 
a measure of the linear dimensions of the condensations which will form in the 
medium on account of this instability; this is Jeans’s result. 

Now suppose that the medium considered is partaking in rotation and that 
there is a term 2u x Q representing the Coriolis acceleration in the equation of 
motion. Then a reconsideration of Jeans’s problem taking account of this 
additional term shows that there are now two modes of wave propagation. 
If V, and V, are the velocities of propagation of these two modes, then one finds 


that 
poeta QJ cos #\2 _. , 
V2V2=4 (“S*) V;*, 3) 


where # denotes the angle between the direction of Q and the direction of wave 
propagation considered. From this relation between lV, and I’, it follows that 
if V,?<o then either V,? or V2 must be negative. In other words if Jeans’s 
condition for gravitational instability is satisfied, then the propagation by one 
of the two modes must be unstable. Jeans’s condition is, therefore, unatfected 
by the inclusion of Coriolis’ force. There is only one exception to this rule and 
this occurs when we consider a wave propagated in a direction at right angles to Q. 
In this latter case, there is only one mode in which a wave can be propagated and 
its velocity is given by 
ei , 47Gp 


V?=475 +¢c?- Re (24) 


Consequently, in this plane the propagation of a wave will always be stable provided 

Q? >7Gp. (25) 

But even when this last condition is satisfied the propagation of waves in other 
directions will exhibit gravitational instability. 

Similarly, we find (cf. Chandrasekhar and Fermi (2)) that Jeans’s condition 


is unaffected also when a uniform magnetic field is present. In this case it is 
found that there are three modes of wave propagation. One of these is 
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unaffected by gravity and compressibility and represents the usual hydromagnetic 
wave of Alfvén. But the other two modes are affected and their velocities are 
related by 
caer H? cos? } 
VV f= —_—_ V ', (26) 
47p 

where # now denotes the angle between the direction of H and the direction of 
wave propagation considered. From this relation it is apparent that if V;?<o 
then either l’,? or V’,? must be negative; and Jeans’s condition follows. 

Finally, if we consider the case when a magnetic field is present and Coriolis’ 
force is also acting, we find that Jeans’s condition is still unaffected. In this 
case there are three modes of wave propagation all coupled in such a way that 
their velocities V,, 1’, and V, are related by 

H?2 cos? #\2 
VPV2V,?= (ra V;, (27 
47p ) 
where # has the same meaning as in (26). And again if V.;?<o, one of these modes 
will be unstable and gravitational instability will ensue. 


VI 

The foregoing discussion of the various problems of stability draws attention 
to an aspect of theoretical investigations in astronomy and geophysics which is 
not often recognized: it is that matter and motions in the large may exhibit 
patterns of behaviour which one might never suspect if one restricted oneself 
to matter and motions in the small. Several examples could be given to underline 
this. Thus, in most hydromagnetic problems which occur in astronomy one may 
treat the matter (be it in interstellar space, stellar atmospheres or stellar interiors) 
as an infinitely good electrical conductor with the lines of magnetic force 
permanently attached to it; and one may do this even though the electrical 
conductivity of stellar material judged by terrestrial standards is by no means 
extraordinary. And the reason why we may do this is simply the very large 
linear dimensions of the systems one deals with, the large dimensions leading to 
times of decay of a prevalent magnetic field which are very long compared to the 
times in which one is normally interested. This applies even to certain aspects 
of geomagnetism. For example, in considering the secular variation of the 
Earth’s magnetic field we may treat the fluid core of the Earth as an infinitely 
good conductor ; for, while the time of decay is of the order of 50000 years, the 
periods in which the secular variations take place are of the order of a hundred 
to two hundred years. This reference to the Earth’s magnetic field reminds 
us that its existence is the best illustration of what a large mass of electrically 
conducting fluid, in rotation, can exhibit. Indeed, rotation must influence 
fluid behaviour in the large in a manner which we are still very far from 
comprehending. Thus, all the evidence we have tends to support the view that 
solar rotation plays a decisive role in the entire complex of phenomena associated 
with the sunspots and the solar cycle; and yet the magnitude of the solar rotation 
is, by all conventional standards, very small. 

It would appear then that the scope for applied mathematics in astronomy 
and geophysics is a very large one. But in order to make progress it may often 
be necessary to make the most severe idealizations, as will have been evident from 
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the examples I have considered. In defence I may quote from Sir George 
Darwin’s address of welcome to the Fifteenth International Congress of 
Mathematicians when it met in Cambridge in 1909: 

“T appeal, then, for mercy to the applied mathematician and would ask you 
to consider in a kindly spirit the difficulties under which he labours. . . . yet 
they are honest attempts to unravel the secrets of the Universe in which we live.” 
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THE MOTION OF MAGNETIC FIELDS 
F. W. Dungey 


(Received 1953 July 20) 


Summary 
‘The temporal change of a magnetic field cannot be described as a motion, 
if it involves interlinking of the lines of force. It is shown that the idea of 
motion can be retained, if the lines of force are regarded as being 
‘* disconnected ’? on a suitable surface. A process is described by which 
magnetic flux can be generated other than at a neutral line. 





1. [ntroduction.—The idea of the motion of a magnetic field has been found 
to be useful in problems of cosmic electrodynamics, because the motion can 
often be taken to be the same as the motion of the material. Sweet (1) has 
extended this idea and has considered the motion of the field relative to the 
material. He does not point out, however, that in some cases the temporal 
change of a magnetic field cannot be described by a motion of the field. If it 
could, it would follow that the magnetic flux linked by any closed curve moving 
with the velocity of the field would be constant, and in particular the flux linked 
by a line of force would be constant. But the rate of change of the flux linked by 


a line of force is given by the e.m.f. pe . ds round the line of force, and this 


need not be zero. It is shown here that the change in the magnetic field can 
still be represented as a motion except on some suitably chosen surface. 
2. Representations of the field.—We denote the magnetic field by H. Sweet (1) 
uses the representation 
H=FV¢a Vu, (1) 
where ¢ and % are functions of position and time and F is a function of ¢ and w. 
The lines of force are the lines of intersection of the surfaces of constant 
¢ and y%. Sweet shows that the motion of the field is equivalent to motion of the 
surfaces of constant ¢ and y%, F remaining the same function of ¢ and yw. 
Corresponding to the difficulty concerning the motion of the field, there is a 
difficulty concerning the representation (1). 
Lundquist (2) has pointed out that the problem is more complicated than 
appears from Sweet’s paper. Lundquist uses a similar representation 
A=V2+dVu 
and 
H=VaA=Vda Vu. 


He then points out that the flux linked by a line of force is given by 


$A.ds= V2. ds, 


since Vy% is perpendicular to the line of force. He concludes that, if this flux 
does not vanish, « is multivalued. We now show how the representation (1) can 
be modified, so that it can still be used when the flux linked by a line of force 
does not vanish. 
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In this case there must be linked lines of force. The simplest field of this 
type is the twisted toroidal field, in which each line of force lies on a toroidal 
surface and is linked with the lines of force inside this surface. In general each 
line of force is infinite in length and covers the whole of the toroidal surface; 
this is possible because the surface is not simply connected and so the line need 
not cross itself. Then the only surfaces containing lines of force are the toroidal 
surfaces. ‘The representation (1) can still be used, however, if a cut is made, 
so that the surface covered by any line of force is made to be simply connected ; 
on a toroidal surface the cut may be taken to be a closed curve linking the toroid. 
Then, supposing that ¢ is constant on the surface, % may vary, but will be 
discontinuous at the cut; F can be chosen so as to make ¢ single-valued. 

3. Motion of the field—The statement that a magnetic field moves with 
velocity u means 

—e curl E=dH ot=curl (uaH), (2) 


but, if the electric field E is given, u is not completely determined by (2); any 
velocity w such that wa H is irrotational can be added tou. Thus the component 
of u parallel to H is arbitrary, and is conveniently taken to vanish. The 
component of w perpendicular to H must have the form H "Ha VW, where 
H =|H| and W is constant on any line of force, but is otherwise arbitrary. Such 
a velocity w does not change H, but just permutes the lines of force. 

The temporal change of a magnetic field cannot always be described as a 
motion of the field. By integrating (2) over a surface bounded by a line of force 
and using Stokes’ theorem, 


~cf E.ds= 5 | H.dS= (unt) ds (3) 


lor a line of force the right-hand side vanishes and hence the magnetic flux 
linked by a line of force must be constant; this may also be seen by picturing the 
change in the field resulting from any motion. ‘The left-hand side of (3), 
however, does not in general vanish, so that the concept of the motion of the field 
is not always valid. 

This difficulty can be overcome by choosing a surface S which is crossed 
at least once by each line of force and by regarding the lines of force as being 
‘“‘disconnected”’ on S. ‘The lines can be pictured as being cut at this surface, 
the two ends arising from the same cut being allowed to move on S with different 
velocities. ‘Their motion must be such that H remains continuous and hence 
their relative motion must just permute the ends of the lines of force; at any 
boundary of S their relative motion must be parallel to the boundary. S may 
consist of a set of separate surfaces and in general these must be given a motion 
such that S continues to intersect every line of force. This picture is of course 
just a mathematical device and S has no physical significance. 

The motion of the field can be defined as follows. Let 


V= f E . ds, (4) 


the integral being taken along a line of force in the direction of H, starting from S 
and not crossing S again. Then, everywhere except on S, E-—VV is 
perpendicular to Hi and (2) is satisfied with 


u=cH-{E-VV)aH. (5) 
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At S both V and u may be discontinuous, but, on both sides of S, 0H/dt is equal 
to —c curl E, which is continuous. 

If, as in the case of the twisted toroidal field, some lines of force cover a whole 
surface, they will cross the surface S an infinite number of times. A line of force 
which crosses S more than once is cut into sections, which could be labelled and 
followed through the motion. The connections at S between the sections are 
altered by the motion, so that two sections, which at one instant form part of the 
same line of force, may belong to different lines of force at another instant. 
Similarly a line of force, which at one instant crosses S only once, may at another 
instant be just a section of a line of force which crosses S more than once. Thus 
the lines of force do not generally retain their identity during the motion, though 
the sections do so. 

The validity of the description of the change in the field as a motion has some 
useful implications. It is also rather convenient that the change in the field is 
separated into two parts represented by the continuous motion and the 
discontinuity, because the change in the linkages between lines of force depends 
only on the latter. The Hall electric field is perpendicular to H, so that it does 
not contribute to V and it is usually the effect of the Ohmic electric field j o that 
is required. The effect of the Ohmic field is well known, it causes a decay of field 
energy and generally tends to untwist the lines of force. This effect is usually 
very slow in astrophysical systems, but may be increased by turbulent motion 
in the gas. 

Sweet (1) used the idea of the motion of the field to discuss the effective 
conductivity in a gas in turbulent motion. Due to the Ohmic field any surface 
generated by lines of force crosses material at a rate inversely proportional to a. 
This rate may be expressed as a surface integral and is found to be increased by 
turbulence; the effect can then be expressed as a reduction in co. Sweet’s 
discussion is restricted to axially symmetric fields and turbulence, and, though 
this is sufficient for certain applications, an extension to three dimensions may be 
of interest. The results obtained here show that the motion of the field could 
still be used, but V would not vanish as it does in the two-dimensional problem, 
and, if the effect were to be represented as a reduction in conductivity, the 
untwisting of the field would also have to be considered. 

4. The generation of magnetic flux.—Some time ago Cowling (3) used a 
formulation, very similar to that of Section 3, to discuss the possibility of a 
sunspot field being maintained by the mechanism of the self-excited dynamo. 
Alfvén (4) and Bondi and Gold (5) have also discussed this problem using the 
idea of the motion of the field. These discussions are concerned with the 
“creation” and “destruction” of magnetic flux. When the lines of force retain 
their identity, the significance of these terms is obvious and it is easily seen that 
they can occur only at neutral lines of the field. Flux flows into or out of a neutral 
line at a finite rate, since wx H~!. If the limiting lines of force near a neutral 
line are circular, the Ohmic field always destroys flux there. At the other type 
of neutral line the lines of force have the form of an X. Here flux is created by 
the breaking of a line of force to form two lines or destroyed by the reverse 
process ; this has been described in detail by the author (6). 

In the general case, when the lines of force do not retain their identity, the 
term “creation of flux”’ is not so easily defined. It may be noted here that the 
flux through an unclosed surface such as a sunspot can in any case be much larger 
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than the “total flux”’ of the field, if the individual lines of force cross the surface 
many times. We now describe a process, not involving a neutral line, which 
clearly creates flux. An attractive mechanism, in which this occurs, has been 
suggested by Alfvén (4). Suppose that two parts of the surface of a tube of force 
come into contact over a finite area; in practice the field will be continuous at the 
interface. Now choose S to be a surface crossing both parts of the tube where 
they are in contact. It is evidently possible in principle for the ends of the lines 
of force to move, so that the ends which originally belonged to one part of the 
tube become connected to the ends which originally belonged to the other part ; 
if this occurs, two separate tubes are obtained, each containing the same amount 
of flux as the original tube. 

This example shows that care is sometimes needed in applying the picture 
developed here. If S had been taken to cross only one part of the tube, the 
correct result would not have been obtained. The reason is that the relative 
velocity of the ends of the lines of force at the boundary of S would not have been 
parallel to the boundary. Care is also needed in discussing any model containing 
finite regions of zero field, since magnetic flux can cross such a region at a finite rate. 


I would like to thank Dr R. G. Giovanelli for the benefit of his comments. 
This work was carried out during the tenure of a Senior Research Fellowship 
at the University of Sydney. 


School of Physics, 
University of Sydney: 
1953 June 4. 


References 


(1) P. A. Sweet, M.N., 110, 69, 1950. 

(2) S. Lundquist, Ark. f. Phys., 5, No. 15, 1952. 
(3) ‘T. G. Cowling, 1/7.N., 94, 39, 1934. 

(4) H. Alfvén, Tellus, 2, 74, 1950. 

(5) H. Bondi and 'T. Gold, W.N., 110, 607, 1950. 
(6) J. W. Dungey, Phil. Mag., 44, 725, 1953. 








THE MAXIMUM EFFECT OF CONVECTION IN STELLAR 
ATMOSPHERES ON THE OBSERVED PROPERTIES OF 
STELLAR SPECTRA. I 


Antoni Przybylski 
(Communicated by the Commonwealth Astronomer) 
(Received 1953 October 3)* 


Summary 


Computations have been made of model stellar atmospheres corresponding 
to spectral types from about Go to Ao. ‘The atmospheres have been 
constructed with radiative temperature gradients and also with adiabatic 
gradients in those regions where they are less than radiative gradients. 

The calculations have been carried so far as to give the monochromatic 
fluxes and colour temperatures of the stars. It is found that the difference of 
colour gradients between “‘ radiative’ and ‘‘ adiabatic’’ models is easily 
accessible to observation (0°58 in the most extreme case). 





1. Introduction.—In spite of many attempts to solve the problem of convection 
in stellar atmospheres, our present knowledge of its effects on the structure of 
the outer layers of stars is still inadequate. ‘The observations of solar granulation 
and theoretical investigations show clearly that there is an unstable zone in 
stellar atmospheres, in which the transfer of energy is partly due to radiation and 
partly to convection and in which therefore the actual temperature gradient has 
a value intermediate between the two limiting values of the radiative and 
adiabatic gradients; but we have no sure basis on which to investigate 
theoretically the influence of convection on the temperature distribution in the 
convectively unstable zone. Since, however, the energy distribution in the 
stellar spectrum is a function of the temperature distribution in the outer layers 
of stellar atmospheres, it can be hoped that the analysis of stellar spectra may 
lead to an approximate evaluation of the influence of convection on the structure 
of the convective zone, if there are any properties of stellar spectra which are 
sensibly affected by convection and which can be easily measured. ‘The present 
investigation shows that the colour gradient of a star satisfies both the above 
conditions and therefore may be used for that purpose. In subsequent papers 
it is intended to deal with other such properties. 

The present investigations consist in the computation of colour gradient for 
a number of pairs of model stellar atmospheres, the models of every pair having 
the same surface temperature and the same value of surface gravity, but differing 
by disregarding convection in one model altogether, and by taking it at its highest 
amount, i.e. the adiabatic value, in the convective zone in the other; the colour 
gradients of real stellar atmospheres must then lie between the two values found 
in that way. Such pairs of model stellar atmospheres were computed ten years 
ago by Strémgren and his collaborators (1), but they soon became obsolete, 
when the investigation of Chandrasekhar and Breen (2) showed that the old 


* Received in original form 1953 June 22. 





684 Antoni Przybylski, The maximum effect of convection in stellar Vol. 113 


value of the absorption coefficient for the negative hydrogen ion used by 
StrOmgren was only a rough approximation. Since they seemed to be inadequate 
for the study of the effect of convection in stellar atmospheres on the observable 
properties of stellar spectra, a new series of models based on recent values of 
absorption coefficients was computed. 

2. Mean absorption coefficient used.—The natural constants c (velocity of 
light), e (electronic charge), m, (electron mass), / (Planck’s constant), k (Boltz- 
mann’s constant), my, (mass of the H atom) appearing in the formulae for 
absorption were taken from the critical work of Du Mond and Cohen (3), whose 
values differ slightly from those used by Strémgren. In order to achieve 
agreement between the adopted values of the natural constants and the observed 
limits of the spectroscopic series of hydrogen the value of 13°595 volts was 
assumed for the ionization potential of hydrogen ; this seems to be confirmed by 
recent experiments (4). 

As in the work of Strémgren (1), only the continuous absorption coefficients 
of the neutral atom and negative ion of hydrogen were taken into account, the 
effect of other sources of absorption being reserved for further investigations. 
The metal atoms were considered only as affording a supply of electrons, and they 
are very important at lower temperatures, where the amount of ionization of 
hydrogen is small on account of its relatively high ionization potential. The 
logarithm of the ratio 4 between the number of hydrogen particles and the 
number of metal particles is assumed to be 3-8, the composition of the metal 
mixture being the same as in the work of Strémgren (1); hence Strémgren’s 
table of the average ionization degree of the metals XY, was used throughout in 
the present investigations. 

The mean absorption coefficient has been computed simply by weighting the 
coefficients for various wave-lengths with Planck’s function. The Planck mean 
has been chosen rather than Rosseland’s mean or Chandrasekhar’s mean for a 
number of reasons. In the first place, the weights for the Chandrasekhar mean 
are given only to optical depth 2, which is insufficient for our purposes ; and again, 
while there was no special reason to use a mean based on the theory of radiative 
equilibrium in the convective zone, yet it was desired to use the same mean in 
both cases. Finally, theoretical arguments have been put forward by Michard (5) 
to show that the Planck mean gives a good approximation to the temperature 
distribution in the radiative case. Since the Planck mean has been chosen 
rather than the Rosseland mean, the tables of absorption coefficients published 
by Vitense (6) cannot be used. 

2.1. Mean absorption coefficient for neutral hydrogen.—The continuous 
absorption coefficient for a neutral H atom corrected for stimulated emission 
is given by the formula 

_ 6474m,e 


—hwl/kT,. —\XgQlkT 
= 3 ach (1 —€ ) e ° D, 


Ky 


where D= 2 Bn xgin'kT 
_aee 


’ 


xu being the ionization potential of hydrogen in the ground state and g, the 
Gaunt factor for the mth quantum state. Defining the Planck mean by 


a(H)= | ag dv, 
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De du. 


O m,e'” | 
0 


we find a(H)= i. —— os 
V3 e4ch(kT)s 


Introducing further the numerical values and the reciprocal temperature 


5040 


6 T 


we obtain a(H) = 3°7429. 10° 46% e ies De“du. 
0 


If D is constant in an interval, that is, if the Gaunt factors are put equal to 1, 
the contribution of this interval to the Planck mean can be found algebraically. 
In the present investigations only the initial terms were summed in D and 
the remainder with »>9 were approximated by the corresponding integral, 
giving for all terms 
er 2uUH 
‘The Gaunt factors computed by Menzel and Pekeris (7) were taken into account 
for n=2, 3, 4,5. ‘The Lyman continuum (”=1) was omitted altogether, since 
it can be expected that the radiative flux in it is negligible. The contributions 
of the four next continua were evaluated numerically with the aid of Cotes’ 
formulae, the corresponding intervals being divided according to need at 
five or more equidistant points; and those of further continua were found 
algebraically. In all, eight continua were computed and the effect of the rest 
was evaluated approximately by extrapolation. ‘The contribution of the further 
continua to the Planck mean is generally higher than their contribution to 
the Rosseland mean. 

The Planck mean a(H) was computed from @=0-18 to 6=0-72 at intervals 
of 0-02 and from @=0-72 to 6=1-08 at intervals of 0-04, and then interpolated 
for intermediate values of @. This interpolation was performed with the aid 
of the function {De~“ du, which varies only slowly in the whole interval of @. 
The results of the computations are given in Table I. 

2.2. Mean absorption coefficient for negative hydrogen ion.—The absorption 
of the negative hydrogen ion is the main source of opacity in lower-temperature 
stars, especially in dwarfs. The numerical value of the coefficient «, of absorption 
of the negative hydrogen ion H_, referred to one neutral hydrogen atom and 
unit electron pressure p, (in dyncm *) and corrected for stimulated emission 
was tabulated by Chandrasekhar and Breen (2). Defining its mean again as 

- 
aH~)= | 


. 
Kk, — dy 
Jo 'S 
and using equation (1) we find 
15 re u® e u 
ajH-)= | « = du. 
( ) a4 Jo v oe “ 
For the sake of convenience (easier interpolation for intermediate values) 
the computations were carried out with values of «, referred to one negative 
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hydrogen ion and then referred to one neutral hydrogen atom by multiplying by 
the factor (0) = 4-158 x 10-10 952 , el 7266, 


giving the number of negative ions present for one neutral atom. ‘The numerical 


TABLE I 


Mean absorption coefficient of neutral hydrogen per atom 


(The quantity tabulated is the common logarithm of a( A) in cm*) 
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evaluations between A=912A and A=22800A were made with the aid of 
Cotes’ formulae, using the same points of division as for neutral hydrogen; 
five equidistant points were used between u=o (A= 0) and the value of wu 
corresponding to A=22800A. Computations were made between 6 =0-50 
and @=0-72 at intervals of 0-02 and for 9 >0-72 at intervals of 0-04. The results 
are given in Table II. 

3. Construction of models.—Models were computed for five dwarfs with 
logyy g=4'5 and 6,=0°6, 0-7, 08, 0-9, 1-0, and for three giants with logy) g=3°5 
and @,=0-6, 0-8, 1-0, 0) being the reciprocal temperature of the stellar surface. 
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The tabulated results of the computations (Tables V to XII) give the values 
of 0, p and p, for a certain number of suitably chosen values of the mean 
optical depth 7. 

Two models for 6)=0-6 were computed in order to extend on a uniform 
basis the available data, though probably they do not correspond to real 


TABLE II 


Mean absorption coefficient of negative hydrogen ion per neutral hydrogen atom and unit 
electron pressure, in 10~*8 cm* dyn™} 





a(H-) a(H-) a(H-) 





0615 2809 : 8-288 
0665 . 2°959 8-620 
0'717 f 3°114 8-961 
0-772 3°275 9°313 
0830 3°442 9°670 
o'89gI 3°614 10°041 
0956 3°794 10°422 
1°024 3°979 10°814 
095 4°170 11215 
170 ‘ 4°369 11-628 
‘250 ‘ 4°575 12°051 
"332 : 4°789 12-485 
“419 5 ‘orl 12°933 
"509 ’ 5°239 13°389 
‘604 5°477 13°859 
‘703 5°720 14°342 
‘806 5°971 14834 
‘914 6°231 15°344 
6-499 15 863 
6°775 16-396 
7061 16°944 
7°353 17°506 
7°656 18-079 
7°967 


— 


I 
1 
I 
I 
I 
I 
I 
I 
I 
I 
2° 
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stellar atmospheres, since at higher temperatures the electron scattering is 
not negligible. It is intended to include the effect of this source of opacity in 
further investigations. 

3-1. Radiative models—The computation of model stellar atmospheres in 
radiative equilibrium consists in the numerical integration of the equation 


dr=-Kdp, (2) 


(1-Xy){a(H) + a(H-)p,} (3) 


I 
where k= — 
My 
is the absorption coefficient per unit stellar mass. The ionization degree Xq 
of hydrogen and the electron pressure p, must be found in each step of integration 
from the value of 6 given by Milne’s formula 


0 =0,(1 + 3r)-44; 
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Milne’s law of temperature distribution was used instead of Chandrasekhar’s 
more recent “fourth approximation” (8) in order to have a common basis for 
comparison of our models with those of Strémgren, who also used Milne’s 
formula. The values of X, and p, were computed with the aid of the formulae 


logyys = —0°477 — 13°5950 + logy, 7, 


?.= PF Xy-s+ {(2j%u-s) + ps (: + yy. 
1—Xy=P,/(p.+5), 


the ionization degree of metals Xy being taken from the tables of Stromgren. 
No use was made of Strémgren’s tables for X, and p, on account of the change 
of ionization potential from 13°53 to 13°595 volts. 

‘The integration of equations (2) and (3) can be effected by choosing p as 
the independent variable and computing 7 from 


I 
&my 





[ @-Xu)faH) + aH p.} dp, (4) 


where the integrand depends only on p and @ (since p,, Xy, a(H), a(H~) are 
functions of p and @, or even @ only). As @ does not vary much for the first small 
values of 7 it is possible to assume 6 =6), to compute the values of the integrand 
for a number of equidistant values of p and then, using Cotes’ integration 
formulae, to find sr. When @ begins to deviate from @, its correct value at each 
step of the integration is found from Milne’s equation, in which the extrapolated 
value of 7 is introduced, and then the corresponding value of the integrand in (4) 
is calculated. From this the integral Av, and thus 7, are found. Should the 
value of z found by the integration differ appreciably from the extrapolated one, 
the corresponding step of integration must be repeated; but by choosing the 
interval of integration reasonably small it is possible to extrapolate + correctly 
and thus avoid the repetition of computations. 

The numerical integration of equation (4) was made with the aid of a Cotes 
formula extended over a number of steps including the new step and several 
previous steps. The integrated value of 7 added to its value at the beginning 
of the interval of integration gives the new value of 7 at the end of the interval. 
Generally each new value of 7 was found by using two different Cotes formulae 
applied to two different intervals of integration; the agreement between the 
values thus obtained is a check on the computations and indicates that the Cotes 
formulae used were sufficiently accurate. This method of integration seems 
to be more convenient than the use of other quadrature formulae, in which it is 
sometimes difficult to include the differences of higher orders. Moreover the 
use of a Cotes formula reduces the effect of summation of errors, which enters 
the other quadrature formulae more seriously. 

The interval of integration, which was chosen small for the initial steps of 
integration, was increased for moderate values of 7, at which it was possible to 
extrapolate with more confidence. Since, however, the increase of temperature 
with the increase of optical depth affects seriously the integrand in equation (4), 
it was finally necessary at larger values of 7 to reduce the interval of integration 
again, usually by halving it when necessary. 
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The method of integration described above gives the optical depth 7 as a 
function of total pressure p. For the sake of convenience the final tables of 
model stellar atmospheres are arranged so as to give 6, log,, p and log,, p, as 
functions of r. 

3.2. Adiabatic models.—T he stellar atmosphere computed under the assumption 
that the temperature is controlled by radiative equilibrium becomes unstable at 
a certain depth r,, at which the radiative temperature gradient 


d logy T oe i he ¢ 
G login P rad ™ 8 ria sid = " (5) 


becomes greater than the adiabatic gradient, defined with sufficient accuracy 
by Unsold’s formula 
, Iy Xy(1 — Xp) ' —. 
(doe r) _ Stdet 34e* ee -Xy)(3 ig er] 
lo adi: 5 Xy(1-X 2° 
Pte P7ome 3a +Xq+ 5Xut yaa +Xy(t—-Xn) ( > xr) 
(6) 

By computing both values of the gradient for several depths in the radiative 
model one easily finds the level p, at which both gradients become equal. The 
stellar atmosphere below this level is convectively unstable and its real temperature 
gradient lies between the limiting values (5) and (6). Since the radiative 
value (5) is valid for the radiative model, in order to find the possible maximum 
effect of convection on the structure of stellar atmospheres, an additional 
‘adiabatic’? model was computed for every radiative one, assuming that the 
temperature gradient in the unstable zone is adiabatic, i.e. defined by equation (6). 
The computation of these adiabatic models consists in numerical integration of 
equation (6), log,)p being the independent variable. 

The process of integration was similar to that for radiative models, the 
intervals for log,,p being chosen, according to circumstances, between 0-01 
(for 6,=1-0) and o-10 (for 6,=0-6). The optical depth in these models was 
computed additionally with the aid of the formula 





logie P 
+ ————— 
gmy loge. 10810 Pe 
which can be derived from equations (2) and (3) and in which 7, and p, signify 
the optical depth and the total pressure at the top of the convectively unstable 
zone. ‘The use of equation (7) instead of equation (4) for the adiabatic models 
is due to the fact that the computations of adiabatic models give the values of T 
(or 0), p, and Xy as functions of log,yp, thus making the formula (7) more 
suitable in practice. 

4. Monochromatic fluxes.—As a preparatory step for the evaluation of mono- 
chromatic fluxes* F,, necessary for computing colour gradients, the monochromatic 
depths 7, were found as functions of total pressure p for two wave-lengths, 
A=4235A and A=6252A. These particular wave-lengths were chosen because 
they are division points of the numerical integration for a(H) and a(H_) in the 
visible part of the spectrum, and therefore their monochromatic coefficients of 


*The symbol F, for monochromatic flux is as used by Woolley and Stibbs in The Outer Layers 
of a Star (Oxford, 1953). In his book Radiative Transfer (Oxford, 1950) Chandrasekhar uses 7F, 
for this quantity, and does not use the symbol ¥, which is employed by Woolley and Stibbs where 
Chandrasekhar uses Fy. 


7(log19P) = 7 (1-Xy)ia(H) + aC )p.jpdlogiop, (7) 


47* 
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absorption «,(H) and «,(H~) were immediately available without extensive 
computations. ‘The monochromatic depths for the radiative models were 
found by integrating the equation (4) rewritten in the form 


I Pp 
r= =| (1 Xu) fes(H) +H ).} dp, 
&My! 0 
and those for adiabatic models with the aid of formula (7) rewritten as 


I 1ogi0 P ; 
Letatias, +i] 1—~Xy){x,(H)+«,(H-)p, }p dlog,, p, 

7 (logiop) =7-(v) gmy loge a, aik(H)+«(H )p, jp dlogiyp 
in which 7,(v) and p, signify the monochromatic depth and the total pressure 
at the top of the convectively unstable zone. 

The monochromatic fluxes F(v) were found by numerical evaluation of the 
integral “84 
F(v)=27 | B,EAr,) dr, 

0 
with the aid of a Cotes eight-point formula applied to the intervals from 7,=0 
to r,=2°8 and from 7,=2-8 to r,=8-4. The contribution of deeper layers to the 
monochromatic fluxes was omitted since it is negligible. Special tables of B, 
were constructed for both wave-lengths in question. 
5. Colour gradients.—The colour gradient is defined by the equation 
d 
-— — 5 Ft 

d= d(a)A) {Ind F;,}. 

From that equation we find as an approximation over a finite interval A, <A <A,: 
_ _ InQy°F(A2)) — In Ay) FA) 





or, introducing F,=(c/A*)F,, 
i 3 In (Ag/Ay) + In (F(v2)/ Fs) 
ei I/A, — 1/A; . 


Introducing further the numerical values A, =4235A and A,=6252A used in 
computation of monochromatic fluxes we find 
$ = 3°0231(0°5074 + logiy F(v2)/ F(¥1))- 

6. Results—The model stellar atmospheres are given in Tables V-XII. 
The upper part of each table gives the reciprocal temperature 6, the total 
pressure p and the electron pressure p, as functions of the optical depth 7 
above the top of the convectively unstable zone; the last entries in this part 
give the critical values 7,, 6,, logigp, and logy)p,.. The lower part gives, 
side by side, the structures of the radiative model and of the adiabatic model. 
The actual computations of models were made to one or two places more than 
are given in the tables. 

The monochromatic fluxes for A=4235A and A=6252A are given in 
Table ITI. 

The colour gradients of the computed models are listed in Table IV. 

For adiabatic models the increase of gradient with the reciprocal temperature 
is quite uniform ; they show also only slight differences in colour gradient between 
dwarfs and giants. 

These regularities do not exist for radiative models. The change of gradient 
with the reciprocal temperature is less regular than in the adiabatic series of 
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models, though the interpolation of gradients for intermediate values of 9%, 
which will be necessary for further investigations, is quite possible. The giants 
have lower gradients than the dwarfs, with the exception of the models with 
) =0°6, when the dwarfs and the giants have the same gradient, corresponding 
to an atmosphere in which neutral hydrogen is practically the only source of 
opacity. 

The differences in colour gradients between the corresponding radiative 
and adiabatic models are easily accessible to observation, especially in hotter 
stars, in which the unstable zone begins at low optical depths; in giants, for 


TasB_e III 


Monochromatic fluxes F, in units of 10~° erg cm~* sec 
v J 












































Radiative model Adiabatic model 
G, logio £ l 
A=4235 A A=6252A A=4235A A=6252 A 
o°6 3°5 97°62 60°98 50°09 47°42 
06 4°5 66°56 60°32 49°69 47°50 
o'7 4°5 55°75 42°01 28-9 32°26 
od 3°5 31°02 25°98 16°64 21°63 
o's 4°5 26°42 24°55 16°99 21°61 
°o'9 4°5 13°37 16:09 IOI! 14°73 
1°O 3°5 TT 10°70 5°94 9°97 
1°O 4°5 7°44 10°63 6-22 10°16 
| 
TABLE IV 
Colour gradients 
Dwarfs (logy) = 4°5) Giants (logy) ¢ = 3°5) 
é< jative f jab: ic “ox adiative d iabatic ° 
6 | Radiative | Adiabatic _— 8, Radiative Adiabatic om 
| model model model model 
| 
1°o | 2°00 2:18 o18 1-0 1:96 2°21 0°25 
og | 1°77 2°03 0°26 
o8 | 1°44 1°85 O'4! o's 1°30 1°88 o'58 
oF 1°10 1°67 O'57 
06 | o'92 1°45 0°56 o°6 0-92 1°46 0°54 
| 





























the same reason, these differences are higher than in dwarfs, with the exception 
of the hottest stars, in which they seem to reach an asymptotic value (about 0-56) 
corresponding to an atmosphere in which the contribution of negative hydrogen 
absorption to the opacity is negligible. 

The calculations shown here would go some way towards deciding the amount 
of convection in actual stars, by referring to measures of their colour gradients, 
if the connection between 4, and spectral type were well established. Unfortunately 
this connection is not well enough known at present to enable us to make a 
decisive comparison with the observations. It is therefore proposed to calculate 
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the strength of the calcium lines in the present models, in order to provide a 
criterion of spectral type. It is intended to construct in this way a diagram 
connecting two observable quantities, spectral type and colour gradient, for the 
two cases maximum convection and zero convection. Comparison with 
observation will then, it is hoped, give a good indication of the amounts of 
convection present in actual stars. 

7. Comparison with the models of Strémgren.—Comparison of our radiative 
models with Strémgren’s work (1) shows that both kinds of models have a 
similar structure. The increase of total pressure and electron pressure with 
optical depth in Strémgren’s models is however steeper, since his values of the 
mean absorption coefficient are lower; for models with 6) between 0-8 and 1-0 
the values of p and p, are higher by a factor of roughly 1-3 than in our models, 
this factor being almost constant throughout the whole atmosphere; for 6) =0-7 
this factor is higher and its constancy is less striking. Another feature in 
common between Strémgren’s models and ours, resulting from their similar 
structure, is that the unstable zones begin at almost exactly the same optical 
depth in both kinds of model; the differences in 7, between our models and 
Strémgren’s do not exceed 0-03. Further, in our giants as well as in Strémgren’s 
the top of the convective zone lies at lower optical depth than in dwarfs. 

8. Acknowledgment.—The present investigations are due to the initiative and 
were made under the supervision of Professor R. v. d. R. Woolley, to whom 
the author is grateful for many valuable discussions and much advice. 
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TABLE V 
Model stellar atmosphere: 0,=0°6, logieg=3°5 





























T 6 logiop loi, 
0°00 0-600 
ool 0°598 '°92 0'92 
0:02 0596 1°50 1°07 
0°03 0°593 1°61 oe 
0°04 o's5gI 1°69 12 
0°05 0°589 hey 1°29 
0:06 0'587 1°80 1°34 
0:07 0585 1°84 1°38 
0:08 0°583 1°88 1-4! 
0°09 o°581 1g! 1°44 
o-10 0579-194 1°47 
o'12 0°576 1°99 .*s9 
O°122 0°575 1°99 ¥°§2 
Radiative model Adiabatic model 
a 
0 logioP logo p, 6 logio P 10810 Pe 
O'14 0°572 2°03 1°56 0°572 2°02 1°56 
o'16 0569 2:06 1:60 0°569 = 206 1-60 
018 0°565 2°09 1°64 0°567 2°09 1°63 
0:2 0562 2'11 1°67 0565 2°12 1°66 
0°25 0°554 2°17 1°74 0°560 2°17 hee | 
0°30 0°547 2‘21 1°79 0°556 2°22 £7 
0°35 0°540 2°25 1°84 0°553 2°25 1°81 
0°40 0°533 2°27 1°88 0°550 2:28 1°84 
0°45 0°527 2°30 1°92 0°547 2°31 1°87 
0°50 O°§22 2°32 1°95 0°545 2°34 1°90 
0-60 O'5il 2°36 2‘01 0°541 2°38 1°95 
0°70 O'501 2°40 2°05 0°538 2°42 1°99 
080 0-493 2°43 2°09 0°535 2°45 2°02 
090 0°485 2°46 2°13 0°533 2°48 2°05 
1‘o 0°477 2°48 2°16 0°'531 2°51 2°08 
ae: 0464 2°53 2°21 0°527 2°55 2°13 
I°4 0°452 2°57 2°26 0°524 2°59 2°17 
16 0°442 2°61 2°30 o°521 2°62 2°20 
1°8 0°433 2°65 2°34 0'518 2°65 2°23 
2'0 0°424 2°68 2°38 0°516 2°67 2°26 
2°5 0-406 2°76 2°46 O'512 2°73 2°33 
3'0 0-392 2°83 2°52 0508 2°97 2°36 
3°5 0379 = _ 289 2°58 0504 2°81 2°40 
4°0 0°369 2°94 2°64 0°502 2°84 2°44 
4°5 0-360 = 299 2°69 0499287 2°47 
5'0 0°35! 3°04 2°73 0°497 2°90 2°50 
6-0 0°337 3:11 2°81 0°493 2°94 2°54 
7:0 0°3260=—_ 318 2°88 0'489 = 2208 2°58 
8-0 0°316 =. 3-24 2°94 0°486 8=3:-01 2°62 
g'0 0°307, = 3-30 3°00 0°484 304 2°65 
10°0 0300 =. 3°35 3°05 0-482 3°07 2°68 
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Tas.e VI 


Model stellar atmosphere: @,=0°6, logy ¢=-4°5 








T 6 logiop logy, 

0-00 0-600 

o'o! 0°598 1‘98 1°42 
0°02 0°596 2:21 1°59 
0°03 0°593 2°34 1°68 
0°04 O'591 2°43 1°75 
0°05 0589 2°49 1°80 
0:06 0°587 2°55 1°85 
0°07 0°585 2°59 1°89 
0°08 0°583 2°63 1°93 
0°09 o'581 2°67 1°96 
o'10 0°579 2°70 2°00 
o'r2 0°576 2°76 2°05 





0°137 0°573 2°79 2°09 




















e Radiative modei Adiabatic model 
6 logioP 10810 Pe 6 logioP logo p, 

O'14 0°572 2°80 2°10 0°572 2°80 2°09 
0-16 0569) = 283 2°14 0569 = 282 2°53 
o'18 0°565 2°86 2°18 0'566 2°86 2°17 
0°20 0°562 2°89 2°21 0°564 2°89 2°20 
0°25 0°554 2°94 2°28 0°559 2°94 aa 
0°30 0547 298 2°35 0°555 2°99 2°31 
0°35 0°540 3°01 2°41 O°551 3°02 2°35 
0-40 e533 3°03 2°46 0549 = 305 2°39 
0°45 0°527 = 305 2°50 0546 3:08 2°42 
0°50 O°522 =—3 07 2°54 0544 3:11 2°45 
o'60 O'sil 3°11 2°61 0°540 3°15 2°50 
0°70 O'501 $33 2°67 0°537 3°18 2°54 
0°80 0493 3°15 2°71 0°534 = 3-21 2°57 
o'9g0 0°485 4°17 2°75 0°532 3°24 2°60 
1'0 0°477 319 2°79 0530 89 3°27 2°63 
r°2 0°464 3°22 2°85 0°526 aay 2°68 
1-4 0452 3°24 2°89 0°523. 3°34 2°72 
16 0°442 3°26 2°93 0°520 3°37 2°75 
1°8 0°433 3°29 2°96 0°518 3°40 a4 
2'0 0°424 $35 2:98 O°515 3°42 2°81 
2°5 0-406 3°36 3°04 O'511 3°47 2°87 
ao 0-392 «34! 3°10 O'507 = 352 2°92 
3°5 0379 = 345 3°14 0504 3°55 2°96 
4°0 0369 = 349 3°19 o'501 3°58 2°99 
4°5 0360 = 353 3°23 0-499 «361 3°02 
5°0 O351 3°57 3°27 0497 93°63 3°05 
6:0 0337 = 3°04 3°34 0493 = 3°67 3°09 
70 0-326 ©6370 3°40 0-490 86371 3°13 
8-0 0-316 3°76 3°46 0488 = 3°74 3°17 
9°0 0307 3°81 3°51 0485 3°77 3°20 
10°0 0-300 8=—- 3.86 3°56 0483 3°79 3°23 
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TasBLe VII 
Model stellar atmosphere: 0,=0°7, logy g=4°'5 





























F 7 6 logiop logio Pp, 
0'00 0°700 
oo! 0°697 2°74 1°16 
0°02 0°695 2°99 1°34 
0°03 0'692 3°12 1°43 
0°04 0-690 9°42 1°50 
0°05 0688 3°29 1°56 
0:06 0685 3°34 1°61 
0°07 0683 3°39 1°65 
0:08 0680 3°43 1°68 
0°09 0°678 3°47 1°72 
o'10 0676 3°50 1°75 
O'12 0°672 3-65 1°81 
O'l4 0-667 3°59 1°86 
0°16 0663 3°62 I'gl 
o-18 0659 = 3°65 1°95 
0°20 0656 3°67 1°99 
o'25 0646 3-92 2:08 
0°256 0°646 73 2°09 
| 
‘ Radiative model Adiabatic model 
6 logioP login p, 6 logiop logo, 
0°30 0638 3°76 2°17 0639 3°76 2°16 
0°35 0°630 3°78 2°24 0633 3°79 2:22 
0°40 0-622 3°81 2°31 0628 3°81 2°a7 
0°45 o'615 3°82 2°37 0624 3°83 2°31 
0°50 0'609 3°84 2°43 0'620 3°85 2°34 
0°60 0°596 3°86 ° 2°53 0615 3°88 2°40 
0°70 0°585 3°88 2°62 o'610 3°90 2°45 
080 0°575 3°89 2°70 0606 3°93 2°49 
0°90 0565 3°90 2°77 0603 3°94 2°52 
ae) 0°557 3°91 2°84 0600 3°96 2°55 
1°2 0'54I 3°92 2:96 0°595 3°99 2°61 
1°4 0°528 3°93 3°05 O'591 4°02 2°65 
1°6 0°516 3°94 3°14 0°587 4°04 2°69 
1°8 0°505 3°94 3°20 0°584 4°06 2°72 
2'0 0495 3°94 3°25 0582 4°07 2°75 
2°5 0474 3°95 3°35 0576 4°11 2°81 
3°0 0457 3°96 3°44 O571 4°14 2°86 
3°5 0443 3°96 3°51 0568 4:16 2°90 
40 0430 =. 397 3°56 0564 = 4°19 2°93 
4°5 0°420 3°97 3°60 0'562 4°20 2-96 
5'0 O'410 3°97 3°62 0559 4:22 2°99 
6-0 0394 83°98 3°65 0555 4°25 3°04 
7:0 0°380 3°99 3°67 O°551 4°28 3°08 
8:0 0369 4°00 3°68 0548 4°30 3°11 
9°0 0°359 4°01 3°79 0°545 4°32 3°14 
10°0 0°350 4°02 3°71 0°543 4°34 3°17 
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Tase VIII 
Model stellar atmosphere: @)=0°8, logy g=3°5 





logiop logo p, 





2°66 0°43 
2°88 0°56 
3°01 0°65 
3°09 o-71 
0:76 
3°21 o'81 
3° 0°85 
0°89 
0°92 
0:96 
1°02 
1-07 
2i2 
E-t7 
I 
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Radiative model Adiabatic model 
0 logiop logiop, 0 logioP logo P, 





‘720 3°61 “48 ; 3°61 
"FEI 63 55 ‘ 3°63 
793 65 ‘ : 3°65 
696 ‘66 ‘ ‘ 3°66 
‘681 3°69 
‘669 4°93 
‘657 3°73 
646 3°74 
636 3°75 
618 so 

603 3°80 
"589 3°81 
‘577 3°83 
"566 3°54 
"542 3°87 
o"5 3°89 
0°506 3°91 
0°492 3°93 
0°479 3°95 
0469 3°96 
0°450 3°98 
0°434 4°00 
0°421 4°02 
0°410 4°04 
0-400 4°05 
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TABLE IX 
Model stellar atmosphere: 0@,=0°8, logio g=4'5 





logiop logy p, 





3°41 0:82 
3°62 0°94 
3°73 1-02 
3°81 ‘08 
3°87 ‘14 
3°92 ‘18 
3°97 ‘22 
4°00 26 
4°03 "29 
4°06 33 
4°10 38 
4°14 44 
4°17 49 
4°19 ‘$3 
4°22 58 
4°26 ‘68 
4°30 Be | 
4°32 "85 
4°35 


mm tt 


‘Q2 





4°36 ‘98 

















Radiative model Adiabatic model 
6 logioP logo Pe 6 logioP logio p, 





0°703 4°36 1°99 4°36 
0696 4°38 2°05 4°38 
0681 4°40 2°17 4°41 
0669 4°42 2°28 4°43 
0657 444 2°37 4°44 
0°646 4°45 2°46 4°46 
0636 = 446 2°54 4°47 
0618 4°47 2°68 4°49 
0603 4°48 2°80 4°51 
0589 449 2°91 4°53 
°'577 4°49 3°00 4°54 
0566 = 4°50 3°09 4°55 
0542 4°50 3°26 4°58 
O°522 4°51 3°41 4°60 
07506 «4°51 3°53 4°62 
0-492 4°51 3°62 4°63 
0479 94°52 3°71 4°65 
0469 = 4°52 3°78 4°66 
O'450 4°52 3°89 
0434 «64°52 3°97 
O'42I 4°52 4°03 
O'410 4°52 4°08 
0°400 4°52 4'1I 


‘99 
“ 
14 
‘21 
- 
oe 
43 
49 


57 
‘61 
‘68 
‘74 
‘79 


®NNNNNN NNN NNNNNNNNN DN 
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TABLE X 
Model stellar atmosphere: 9,=0°9, logy) g=4°5 





logioP logo p, 





3°79 0°36 
3°97 0°50 
4°08 0°58 
4°15 0°65 
4°21 “70 
4°25 ‘75 
4°29 ‘79 
82 
‘89 
04 
‘99 
04 
09 
13 
23 
31 


“46 
53 
‘60 


° 
wu unm rm OOOO 0 0 6 Oo 


02000000000 





71 

















Radiative model Adiabatic model 
6 logiop logo p, 6 logiop logo Pp, 





0-767 ; 1°73 0767 = 4°73 
0-752 1°84 0°754 = 4°75 
©°739 ; 1°94 ©0°743 4°77 
0°727 ; 2°03 0735 44°78 
0-716 _ 2°12 0728 §=6. 4°79 
0696 : 2°27 0-717. «4°81 
0°678 : 2°41 0:708 =. 483 
0663 8: 2°53 0-700 4485 
0649 . 2°63 0694 84°86 
0636 . 2°73 0688 4°87 
0-610 . 2°94 0°677 4°89 
0°588 ‘ 3:11 0'669 4°91 
0569 3°25 0662 4°93 
0°553 : 3°37 0656 = 494 
0°539 3°47 o651 = 4°95 
0°527 ; 3°57 0646 89496 
0°506 3°72 0639 «= 498 
0489 ; 3°85 0°633 500 
0°474 : 3°95 0-628 5 ‘Ol 
0-461 : 4°03 0623 5°03 
0°450 ; 4°10 0619 5°04 


"73 
$3 
‘91 


N NNN NKRN KN NKR DNDN DN & & 
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TABLE XI 
Model stellar atmosphere : 945=1°0, logi9 g=3°5 





logio P logio Pe 





3°37 
3°54 
3°64 
3°70 
3°76 
3°80 
3°84 
3°87 
3°90 
3°92 
3°96 
4°00 
4°03 
4°05 
4°08 
4°13 
4°17 
4°20 
4°22 
4°25 
4°27 
4°30 
4°32 





4°32 

















Radiative model Adiabatic model 
6 logiop logio P, 6 logiop 10810 P, 





o°821 4°34 I'I4 0822 4°34 
0808 = 4.36 1°24 o-'811 4°36 
°°795 4°37 1°33 0802 4°37 
°°773 4°39 1°50 0-788 =. 439 
©0'754 4°40 1°64 0777.0 4°41 
0736 = 441 1°77 0768 4°42 
0-721 4°42 1°89 0'760 4°44 
0707 4°43 2°00 0°754 4°45 
0677 4°44 2°22 O'741 4°47 
0653 4°45 2°40 O'731 86449 
0632 3= 44°45 2°56 0°'722 4°50 
o'615 4°45 2°70 o'716 8 =§=64°52 
0599 84°46 2°81 O'710 89453 
0586 = 44-46 2°92 0705 4°54 
0562 4°46 3°09 0696 §=4°56 
0543 4°46 3°24 0690 = 44°57 
0527 446 3°35 0°684 4°59 
o'512 4°46 3°45 0679 = 460 
0500 | 446 3°54 0674 84°61 


14 
‘22 
‘29 


>> 
1 CO 


N NNN SH toe 
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TaBLe XII 
Model stellar atmosphere : 945=1°0, logy g=4°5 
T 6 logio P logio p, 

0°00 1°000 

o'ol 0996 3°92 o'l2 
002 0993 4°08 0°27 
0°03 0989 4°18 0°35 
0°04 0986) = 44:25 0°42" 
0°05 0982 4°30 0°47 
006 °°979 4°34 O°51 
0°07 0975 438 O°55 
0:08 0'972 4°41 0°58 
0-09 0969 44°44 0-61 
o'10 0966 4°46 0°64 
O'12 0'959 4°50 0°68 
O14 °0°953) 454 o'72 
0°16 0948 4°57 0:76 
o'18 0°942 4°60 0°80 
0:20 0°937 4°62 0°83 
0°25 0'923 4°67 o'90 
0°30 o-9gll 4°72 0°97 
0°35 o°900, 4°75 1°03 
0"40 0889 = 4478 1°08 
0°45 0879 4°81 I-14 
0'50 0869 4°83 I‘19 
‘60 0°852 4°36 1°30 
0°70 0836 =. 489 1°39 
0:80 o'821 4°92 1°49 
0877 | o81r 4°93 1°56 























" Radiative model Adiabatic model 
] logiop logio P, 6 logioP okie Pe 
o'9 0808 = 494 1°58 0808 = 4494 1°58 
1'0 °'795 495 1°66 ©7907 = 495 1°65 
1-2 0'773 4°98 1°82 0'779 4°98 1°78 
1°4 0°754 5°00 1°96 0°766 5°00 1°88 
1°6 0°736 5°01 2°08 0°755 5°02 1°96 
1°8 o°721 5°02 2°20 0°746 5°03 2°02 
2°0 0°707 5°03 2°30 0°739 5°04 2°08 
2°5 0°677 5°05 2°53 0°724 5°07 2°20 
3°0 0°653 5°06 ray 0'713 5°09 ae 
i 0632 5°06 2°87 0°704 5°10 2°37 
4'0 0615 5°07 3°01 0696 5°12 2°43 
4°5 0599 507 3°13 0690 = 513 2°48 
5°0 9°586 3=5, 08 3°24 0685 5°14 2°53 
6-0 0562 5:08 3°44 0675 5°16 2:60 
70 0°543 5°08 3°56 0668 “a7 2°67 
8-0 0°527. 5508 3°68 0661 5°19 2°72 
9:0 O'512 5°09 3°79 0656 8=5"20 2°77 
10°'0 0°500 §= 509 3°88 0°651 5°21 2°81 
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Summary 

A successive perturbation theory is established for the structure of stars 
in uniform rotation with a known law of energy generation. The first-order 
structure is computed for a rotating Cowling model with a law of energy 
generation expT™, 

The results from a first-order model indicate that with the distribution 
of rotation velocities observed, rotation can produce a spread in luminosity 
of order half a magnitude in a mass-luminosity diagram, and of order half 
a magnitude for a given spectral type in a colour-magnitude array, this figure 
being for the most rapidly rotating Be stars. 





1. Jntroduction.—The effect of rotation on a star of given mass is not known 
completely, even to the first order in the ratio of centrifugal force to gravity. 
The rotating polytropes of Chandrasekhar (1) give some idea of the redistribution 
of the mass of the star, together with the distortion of the equipotential surfaces ; 
the work of Slettebak (2) on O, B and Be stars includes a discussion of the relative 
variation of effective temperature over the surface using von Zeipel’s theorem, 
while recently (§) he has investigated the dependence of axial rotation on 
spectral type and luminosity among early-type stars and found a difference in 
axial rotation between the main-sequence stars and the giants and sub-giants. 
The distribution of observed values of wv sin 7 for these two groups appears to be 
quite different, the more luminous stars having generally smaller rotational 
velocities. G. R. and E. M. Burbidge (12) have constructed models of rotating 
Be stars with surrounding electron-scattering atmospheres; Opik (3) and 
Sweet (4) have considered the circulation currents which would be produced by 
rotation in a star of uniform composition. 

One effect outstanding at present is the change in luminosity and mean 
effective temperature of a star of given mass due to its rotation. For this purpose 
a knowledge is required of the absolute changes in temperature and density at 
given points in a star of given mass, rather than the relative distributions. It is 
well known that the central temperature and density of any model of known 
chemical composition cannot be computed from the mass until the source of 
energy is known. Normally the central conditions can be assessed if the radius 
of the star is known, but this procedure is not available when considering the 
effect of rotation, since the change of radius is itself one of the unknown quantities. 
For this reason the changes in central conditions due to rotation cannot be found 
without knowing the law of energy generation. The changes in central density 
of Chandrasekhar’s rotating polytropes (x) are derived from an arbitrary 
polytropic relation, and cannot be used for a real star. 
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The object of the present paper is to provide a model of reasonable accuracy 
from which these changes can be derived. The effect of rotation on a given 
star depends on the ratio «=equatorial centrifugal force/surface gravity. An 
analysis of « for the various spectral types is given in Table I. 


TABLE I 
Ratio «= equatorial centrifugal force/surface gravity 

Spectral type Caverage max 
Os-9 0:06 0°36 
B2e—se 0°32 0°89 
B2-5 0°23 0°42 

V o'12 0°53 
aatom { I-IV 0-04 0-31 
A2-5 0°05 0°27 
Fo-2 oO'o! ols 


T he values given apply to the equatorial region of the surface of the rotating star ; 
thus % = yv»?/Rg,(0), where vp is the equatorial rotation velocity, R and g,(0) are 
the radius and surface gravity of a non-rotating star of the same mass, while x is 
a correction factor derived later in the paper, included to allow for the changes 
in equatorial radius and surface gravity due to rotation. ‘The values of vp are 
derived from the values given by Slettebak (2, §) and Miss Westgate (6), and the 
masses and radii involved are based on the work of Kuiper and of Stebbins and 
Whitford (7). It appears later that the first-order model breaks down when 
7 >0°29; the table therefore shows that a first-order perturbation theory is likely 
to be sufficient to describe the average O, A and F types but that the theory 
should be extended to the second order for B stars and the most rapidly rotating 
O and A types. Such an extension is needed to investigate limb darkening in 
these stars, since von Zeipel’s relation H «xg does not hold in the second order ; 
the analysis made by Slettebak (2) would therefore not be valid. It is also 
required in order to examine more carefully the question of possible circulation 
currents due to rotation. Opik (3), for example, appears to find a singularity 
at the surface in the second-order term in the radial component of the circulation 
velocity, while Mestel (8), using a first-order theory, has shown that the forces 
which were previously regarded as producing circulation currents are insufficient 
to overcome the choking effects of the non-uniformity of molecular weight 
produced by hydrogen consumption in the core. 

The second-order theory is rather complex; a method will therefore be 
established for the first-order structure in such a way as to facilitate the passage 
to the second order. For this purpose the initial zero-order model must be as 
self-consistent as possible; a Cowling model of high accuracy is available from 
Miss Gardiner’s integration (9), and this is the one adopted. 

2. The structure equations.—The equations of state, hydrostatic equilibrium 
and Poisson’s equation in a uniformly rotating star are: 


P=p&T/m, (2.1)* 
grad P=p grad yf, (2.2)T 
V2b= — 407Gp + 202, (2.3) 





* Radiation pressure is neglected. 

+ From the results of an earlier paper (4) it appears that the meridian-plane component of the 
circulation velocity would have a negligible effect on hydrostatic equilibrium as compared with 
that of rotation. 
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where P, p and T are the pressure, density and temperature, m is the mean 
molecular weight in units of the mass of the hydrogen atom and @ is the angular 
velocity, while # and G are, respectively, the gas constant and the constant of 
gravitation. « is the total effective potential and is given by 


p=h+$rPQ{1— P,(u)I, (2.4) 
where ¢ is the gravitational potential, r and @ are the radial distance and angular 
distance from the axis of the star and P,(u) is the Legendre coefficient, 
». denoting cos @. 

In the convective core the adiabatic equation 
px 7% (2.5) 
holds, while in the envelope the equation of radiative transfer 

4aT? grad T= —xpH/c (2.6) 
is satisfied, where a is the coefficient of Stefan’s law, « is the opacity, c is the 
velocity of light and H is the radiative flux. Hi must also satisfy an equation of 
thermal equilibrium. In a non-rotating star this is simply L = 47r?H, where L is 
the luminosity. When the star is rotating, convection currents may be set up 
and some of the energy transferred by these. If v is the meridian-plane velocity 
of these currents, the modified thermal equilibrium condition can be written 


L=| [H+ {iv?+ $r2Q%(1— P,)+4-+c,8T/m)py] . dS, (2.7) 


where S is an equipotential surface, %=constant in the envelope. This 
expresses the fact that the total rate of outflow of energy from an equipotential 
surface is contributed by the radiative flux H and the kinetic, potential and 
thermal energies }v? + 47°02?(1 — P,), d and c,4T/m respectively, where c,, is the 
specific heat at constant volume. The contribution }v? is neglected for reasons 
similar to the neglect of the effect of the meridian-plane circulation velocity on 
hydrostatic equilibrium. By using the relation in (2.4) equation (2.7) can be 
expressed in the form 


L =| (H +(u+¢e,AT m)pv] . dS. (2.8) 


The total rate of outflow of material from S must vanish, hence 
| pv .dS=o. (2.9) 
7§ 

Remembering, moreover, that ¢/ is constant over S, (2.8) can be further reduced to 


L=| [H+c,#Tpv m] . dS. (2.10) 


Equations (2.9) and (2.10) hold for all equipotential surfaces in the radiative 
envelope, and, together with (2.1), (2.2), (2.3), (2.5) and (2.6), constitute the 
structure equations of the star. 

Finally, an energy equation is required; in a point-convective model this is 


L= | pe dr, (2.11) 


~ core 
where « is the rate of generation of nuclear energy per gram. In the present 
paper the power law «xp7™ is adopted. Hence 
L=| p*T™ dr. (2.12) 
core 


45 
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3. The successive perturbation theory.—In what follows a star of fixed mass 
and uniform molecular weight m is considered. At any given point the structure 
quantities are therefore functions of 2. In order to work with a dimensionless 
quantity a parameter A is introduced in place of 2 and given by 

A= ?/27Gp,(0), (3-1) 
where p,(0) is the central density in the non-rotating star. The method adopted 
consists of expanding the structure variables as power series in A; equations for 
the coefficients are then derived from the structure equations and their 
derivatives with respect to A. This procedure is equivalent to successive 
perturbations; it is, however, more straightforward than the direct method, 
particularly when deriving the boundary conditions at the core surface and at 
the surface of the star. 

Thus, if Q is any quantity depending on A, Q is expressed in the form 
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ns A" 
O(r, 0; A) = _ O,(7, @) n! ’ (3.2) 
where 7 
O,, =[0"Q/0A"],-o- (3-3) 


The functions Q, are determined by taking A=0 in the structure equations, and 
are, of course, the quantities arising in the ordinary Cowling model. Having 
determined these, Q, can be derived by differentiating the structure equations 
with respect to A and putting A=o. ‘The boundary conditions for the Q, at the 
centre, at the core surface and at the surface of the star are derived by 
differentiating the general boundary conditions with respect to A and putting A=o0. 
Thus, conditions of the type Q=o at the surface of the star are written 
O(r,, 8; A)=0, where r=r,(0,A) defines the surface. On differentiation this 
leads to the condition QO, + (dQ,/dr)r,,=0 at r=ry,, where ry, = R the radius of the 
non-rotating star and r,,=[0r,/0A],-9. Similarly, the type of condition where 
Q is continuous at the core surface r=r,(0, A) makes QO, + (dQ,/dr)r,, continuous 
at r=Vro,, where 7, is the radius of the core in the non-rotating star. In this way 
the first-order model tabled in this paper was derived. The higher-order models 
may be derived similarly by continuing the process. 
Dimensionless variables are introduced by the following substitutions : 


r=r,x, where r,2=5#T(0)/8tmGp,(0), | 

P= PA), 

p=p-(A)o, ' (3-4) 
T=T.AA)t, 

= 5AT,(A)(q—1)/2m +4,(A), 





where P(A), p,(A), T.(A) and (A) are the central pressure, density, temperature 
and total potential in the rotating star, and p, o, ¢ and g are the new dimensionless 
variables. 

The structure equations (2.1), (2.2), (2.3) and (2.5) then reduce successively to 


Pp = ot, (3-5) 

grad (x, )P = 30 grad, a4» (3-6) 
TAA) oo PAA) _ e 

To) ¥ @0t 5 (yo (3-7) 

o=t*? in the core. (3.8) 
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VI 


Equations (2.9) and (2.10) can be written 


* or, 
| pr. al 50 -— (3-9) 
Ory 
=47 fr rH, — <? Ho+ 1 3G pey,AT m\dpu, (3.10) 


where r=r,(6, A) defines an aaamneiae surface and (H,, H,) and (v,, v) are 
the components of H and v evaluated on r=r,,. 
‘The boundary conditions are: 











o=t=q=1; grad o= grad t=grad g=0 at x=0, (3.11) 
o, t, g; grad o, grad ¢, grad q continuous on x=.x,, (3-12) 
e=t=0 68 c=, (3.13) 


where x=.x,(6, A) defines the interface between the core and the envelope and 
x=x,(0, X) defines the surface of the star. ‘The boundary conditions for g at the 
surface are derived from the condition that ¢ and its derivatives are continuous 
with the gravitational potential outside the star for a star of the same mass. 

4. The zero-order model.—On putting A=o the structure equations, with 
Kramers’ opacity, reduce to 





Aces) - dqy 
= 39 “dx ’ 
ax 
I és a 
x dx ” 
dt, a ko,* 
dx — t,I3'2x2” 


where k is a constant. These are the equations for the normal Cowling model, 
and it is convenient to rename op, ¢, and dq,/dx in order to conform with the 
notation used in Miss Gardiner’s numerical solution (g). Thus, take ¢,=y, 
oj=s and —x*dq,/dx=w. 

5. The first-order model.—On differentiating with respect to A and putting A=o 
the structure equations become 


pi=sht+yor, (5-1) 
Op, /Ox = —50,%/2x? + §20q,/dx, (5.2) 
Op, /00 = $s 0q,/00, (5.3) 
V°q, +0,=1-Xsz, (5.4) 
where 
X =p,'(0)/p,(0) — T.’(0)/ T.(0). (5.5) 


The suffix fig. 1 denotes the operation [0/0A],_9 as in (3.3). The above equations 
hold throughout the star. In the core (3.8) yields the further equation 

o, = 3y!*t, in the core. (5.6) 
In the envelope (2.6) gives 


1s TO) pe(0) 1s 
H,=H,| ¥ 7 T(0) ~ po) *? 








o #,. 
3 - 2 “| + w = grad f, (5-7) 
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where Kramers’ law is taken for the opacity. Again, from (2.9) and (2.10) the 
relation 


“1 
L(0) =47"9| A,, du (5.8)* 
0 


can be derived. By substituting for H,, from (5.7), noting that Hy =[ L(0) 477°, o], 
the above equation reduces to 











wh 6% Ot, /Cx > a ey 
j (3 7 "s = dy/ dx dnl, (5-9) 
where L'(0) To) 0) 

_ LO) 5 Te(0) , , e'(0 ; 

~ Lio) 2 To) +2 (0) ’ (5-10) 


(5.9) is the counterpart of the equation of radiative transfer in the Cowling model. 
The boundary conditions, derived from (3.11), (3.12) and (3.13), are: 


o,=t,=9,=0; 00,/0x=0t,/0x=09,/0x=0 at x=0, (5.11) 
dz +2 dy w ] 
yt FT Mes fy “ter M1 33 1c) ; 
dx * dx continuous 
iia d*3 3. 2: d*y 30,/2 ee at (5-12) 
Co, OX+ F3* t,/Ox+ — Tak Ste 11) x + (2t/x? — 3)x4,; las 
00,00, Ot,/60, 0q,/00, J 


where xj, is the core radius in the normal Cowling model and +,,(9, A) =[éx, CA] ,_93 
and 
a, + (dz/dx)x,,=0; t,+(dy,dx)x,,=0 at v=Xp,, (5.13) 
where xy, is the radius of the star in the Cowling model and x,,=|0@x, 0A};<». 
From the Cowling model it is seen that dz dx, dy dx and w are continuous at the 
interface; (5.12) therefore shows that 
1, ty, Gy, Cg, Ox are Continuous at x= X%p,. (5-14) 
d*z/dx* and d*y dx* are not continuous in the Cowling model, hence the fourth 
and fifth members of (5.12) only serve to derive x, at the end of the problem. 
Similarly, dz dx=o and dy/dx+o at the surface of the Cowling model, hence 
the first member of (5.13) gives 
o,=0 at X=X%,, (5.15) 
while the second member will be used to determine x,, at the end of the problem. 
The boundary conditions for g, are derived from the boundary conditions 
for g referred to at the end of Section 3, and are found to be 
o 14 es -, 














on x2(r — m ) 

ae T.(0) ae (5.16) 
0q, _ T,(0)w a6 IP, ‘ 

x 7." T (0) * +4x°(1- ’ (5-17) 


at x=Xy,, where 





z- 2mp,(0) TT. (0) _ vel an 
5AT (0) T.(0) 
and 4 is an unknown constant. 
1 
*If the integrand on the R.H.S. of (2.10) is denoted by J, then dL /dd=4n[ dI/dd. dp. 
~ 4 
Now I=I(r,, 0; A), hence dI/dA=0I/@A+ él /éry. ér/2A. (5.8) follows on noting that Hog=o, 


el 
r°Ho,=constant, Vg=0 and] vi~du=o, 
~ 9 
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The solution of the first-order structure equations can now be undertaken. 
(5.16), in conjunction with the structure equations, shows that g, can be 
expressed in the form g, = a(x) + 6(x)P,, with similar expressions for p,, o, and f,. 
It is found convenient to write a(x) and b(x) in the forms 

a(x) = —}Xwx+)x?+}A(x), b(x)= —)x?k(x), 
thus . . ‘ 

: gi= —$Xw/xt dx? + bA— beRP,. (5.18) 
By substituting these expressions in (5.1) to (5.4) it follows, after some elementary 
algebra, that 














Pi=32(4 — $1 (x), (5.19) 
x22’ ? xtz 
a= — == (@- We) +), (5.20) 
xy’ < - 
= — = @- a(x) -8 2 PW), (5.21) 
where \, [ and & satisfy 
YP xis’ 2! a 
avs 2a nut os oe Ser 9, (5.22) 
w w u w 
k Ck - = k=0, (5.23) 


dashes denoting inicsaitsie en respect to x. ‘These last two equations are 
satisfied throughout the star, hence k, which occurs on its own, can be determined 
immediately the appropriate boundary conditions are known. (5.18) in 
conjunction with (5.11) shows that kx?--0 as x->0. ‘The two solutions of (5.23) 
near the origin are k= constant and kx x °; at the origin the appropriate boundary 
condition is therefore that k should be finite. The second boundary condition 
necessary is obtained from (5.16) and (5.17). By equating coefficients of P, and 
eliminating the unknown 4 these give 

5(k—1)+xk’=0 at x=Xp,. (5.24) 
The solution was then evaluated numerically. One further equation is necessary 
in order to determine A and I’. In the core this is provided by (5.6) which 
reduces to 7 ‘ 

I’ =o in the core; (5.25) 
thus I'=constant in the core. Noting that g,=p,=0 at the origin, (5.18) and 
(5.19) show that T=A=o0 at x=0; therefore [=o throughout the core. 
(5.22) thus reduces to 

. 2. eee xis’ . 
A” + . A’ - A= a the core. (5.20) 
xX v 








Since x*s' w is finite at the origin the two complementary functions behave as 
A=constant and Ax.wx! at the origin. The second of these is unsuitable since 
A=o at x=0; the solution required is therefore the particular integral 
commencing with (=o at x=0. This solution was derived numerically. In the 
envelope .\ and I’ are governed by (5.9), which, after a little reduction, gives the 


relation 

xs -y\. yy, xy’ xy’ xs 

ae oe ome fe OY 4 ee As — (F—A)+ —[2- — 
w -* y wy y 





alt 
a 


bwy’ = 
=—- 1X+Y). (5.27) 
xy —gow 9 ) (5 4) 
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The conditions (5.12) show that A, A’, I and I” are continuous at x=). The 
required boundary conditions at the core surface are therefore given directly 
by the core solution at x=%p,. (5.22) and (5.27) cannot however be integrated 
directly from the core surface because the constant (}X + Y) in (5.27) is as yet 
unknown. By substituting A unr Y)§ w/x+land [=(3X+ Y)By+h, (5.22) 


and (5.27) reduce to xis! 
= <1 +5 — =. (A I+ = —W = ss (5.28) 
” 2 = x3 17- , " ial es — » 
h + “ “2)\n- . + — — (A—l)= ¥(= 2). (5.29) 


The above equations have a singular point at the BAS where the four 
complementary functions behave as /~[(1/z), 1, y, vz]; A~((1/yv2?), I, v, y?3]. 
The boundary conditions (5.13) and (5.16) show that / and h are finite at the 
surface, thus the solutions /~(1/z), h~(1/yz") are rejected. An inspection of 
the equations shows that one solution is (/=1, h=1); the two remaining linearly 
independent complementary functions (/,, 4,) and (/,, h.) were computed by 
inward integration from the surface, using suitable starting series in the variable 
£=x,x—1. For this purpose the Cowling variables (w, y, z) were also expanded 
in series. A particular integral (l,, 4,) was computed in a similar manner. 
A and T can then be season in the forms 
=(4X+ Y)8 (w/x) +1, +21,+Bl,t+y, 
a VY) y+h, +h, + Bhaty. 
The constants [(}X + Y), «, 8, y] are then determined from the four boundary 
conditions at the core surface. In particular it is found that 
($X + Y)= +0°312. 

The solution (A, I’) is then determined in the envelope. 

The constants [X, Y, L’(o)/L(0), T.’(0)/T.(0), p,(0)/p(0)] may now be found 
in the following manner. The substitution of g, from (5.18) into the boundary 
condition (5.17) shows that 

P on, ew . 
A, = =| T.(0) -1x| ’ (5.32) 

where A,’ is known from the <i solution (A, ['). Thus (5.5), (5-10), (5-31) 
and (5. 32) are four equations in the five unknowns required. The additional 
relation necessary is derived from the energy equation (2.12). By differentiating 
this logarithmically with respect to A and putting A=o it can be shown, after a 


little reduction, that oie 20, , 17ty 
, , ' * gtylt ( <2} 2) x2 deed 
Lo) _, (0) _ 7:0) _ Joly *? ( sy )s wie 


(5.30) 


(5.31) 


Oy 























; 735 ~*-— : (5.33)* 
L O (O gf O om 
(0) p-(9) (0) |  ayitaAds 
Jo 
By substituting for o, and t, from (5.20) and (5.21), and using the core relations 
3'/3 =3y'/2y= — 3w/2xy, (5.33) can be reduced to 
-“ 
L'(0) _, pe(0) T.'(0) 20] yl — LXxw + gx? A + grt) de ‘is 
“ 5+. 


—2 _ 17 = 
L(o) pO) T,{0) ee 
J, yx2dx 





* It has been assumed here that the energy production at the core surface in the non-rotating 
star is negligible. 
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With the aid of the core solutions for \ and the Cowling variables, the necessary 
integrals on the right-hand side of this equation may be computed numerically. 
The resulting relation is 
Lo) _, ec(0) _ 1» Te(0) _ 
10) ~~ po) *7 To) ~ 
Equations (5.5), (5.10), (5.31), (5.32) and (5.35), when solved, give 
L'(o)/ L(o) = — 22:14, 


—1-402X + 1-472. (5-35) 





X= 5°842 
T,(0)/ T.(0) = — 1°427, (5.36) 
y 


— 2°609 


Pe (0) p.(0)= + 4°415. 
The functions g,, t; and o, may then be computed; the resulting model is 
tabulated in the Appendix. 
The figure of the star is governed by x,,; the second member of the boundary 
conditions (5.13) gives 
tis 


a (5.37) 


- 


Here ¢,, can be obtained _ (5.21), thereby showing that 


. 2: Te , 
\,-T.)+} ‘of 
MT) +4 


X},= 


X11. r y 


*os 





Se vo pX+4 +3— + RP 


1 
%, j%R 
= pe ty (5.38) 
The surface gravity due to the material of the star is given by 
g,?= grad ¢,. grad ¢,. 
On differentiation with — to A this relation shows that 


3 
a| $= “ue 4 “0 (P,—A,) +4 (2h, + bp hy —1)Py 


= — 29°44+ ona (5.39) 
where the chain of relations (2.4), (3.4) (last member) and (5.18) was used to find 
the derivatives of ¢, and (5.38) provided the on for x,,. Similarly, 


dg ft dg x, 
1 8 — “ss os 7 
808 dn 4 m =fos dr .. “7 tha Py) 


= —66:20+ a (5.40) 
where g,"" is the effective surface gravity. 
A similar a >i for H, the surface flux of — shows that 


1 Hor See 
| See] = GO ena BE a 4d eh ba DP, 


= — 45°74 +78 ‘ah (5-41) 

equations (5.7), (5.8) and (5.21) being used for the derivatives of H,. 
6. The interpretation of the model.—In order to apply the model to stars as 

observed, a parameter ~,, defined by 
R8Q? _, p,(0) 
GM ~*p (0) oe 
is introduced, where p,,(0) is the mean density of the non-rotating star; «, is 
therefore the ratio of centrifugal force at a distance equal to the radius of the 

















= 
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non-rotating star, to the surface gravity of the non-rotating star. In the Cowling 
model p,, p= 0°0272, thus 
A=o-01814%). (6.2) 


Asummary of the results of the previous section, in terms of z,, is given in Table IT 
below. 
TaBLe II 
Summary of results to first order in a, 
L=L(o) [1—o-4014,] 
Myo = Moi) -++-0'569a, (rotation axis L, line of sight) 
—0427a, (pole-on view) 
Pe=Pe(0) [1 +0°0799%] 
T= To) [1 —0°0259%,] 
H,=H,(o) [1—0°828a, + 1-4200,P3] 
£.=5(0) [1 —0°5330, + 0°7540,P 2] 
gett = g.efl(o) [1 — 1-2000, + 1°4200,P 2] 
rg=R [1 +0°2140, —0°3430,P 9] 
—o0'065a, (rotation axis L, line of sight) 


T,=logy T 
logie Te= loge 7 (0) — 0-040, (pole-on view) 


The values for M,,, and log 7, require some comment; owing to the temper- 
ature variation over the surface the star will appear brighter and hotter when 
viewed pole-on than when viewed at right angles to the rotation axis. The 
results were obtained by integrating the emergent intensity in the appropriate 
directions over the appropriate parts of the surface, assuming a normal coefficient 
of limb darkening 0-6. 

The parameter « tabled in Section 1 from observations of stars is given by 


a — 7(7/2,A)/R 


% B42, A)/ Bog 
=1+1-2952, (6. 3) 


a 
to the first order in ~,, where 7, R and g,/g,, are obtained from ‘Table II. Further, 


remembering that v)=7,(7/2,A)Q, the correction factor y used in calculating « 
from observation is therefore given by 


X=14+0°523%, (6.4) 
to the first orderin«,. ‘Thus « and «, may be expressed in terms of the observable 
quantity z’ =v »?/Rg,, by 

a=o'(1+0°523a’), (6.5) 
%,=«'(1+0°7722’) (6.6) 
to the first order in «’. 

The relative displacements in a H—R diagram for the limiting case «=1 as 
derived from ‘Table II are shown in Fig. 1.* The mean effective temperature 
of a star is reduced as seen from both pole-on and axis perpendicular viewpoints. 
‘The apparent luminosity is reduced when the star is viewed perpendicular to the 
rotation axis and, although the total light output is reduced, the apparent lumino- 
sity is increased when the star is viewed pole-on. It can be seen that the main 
sequence is brightened by about 0™-5 for a set of stars viewed pole-on, with a 
smaller increase for stars viewed at right angles to the rotation axis. Since a 


* The diagram indicates only the order of magnitude of these displacements. 
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** pole-on” star would show no apparent rotation as judged from line broadening, 
no definite separation of stars according to rotation would be expected. On the 
observational side there seems to be no conclusive correlation between luminosity 
and rotation velocity. Eggen (10) appeared to find a brightening with rotation 
for the late A and early F types in the Hyades and Pleiades, with the reverse effect, 
however, for early Atypes. The limiting spread of half a magnitude is approached 
only by the most rapidly rotating Be stars. ‘The observed spread for B-type stars 
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Fic. 1.—Limiting displacements of star of given mass and resulting spread in H-R diagram. 


of luminosity class V in clusters, as illustrated for example in Roman and Morgan’s 
work (11) on the Perseus cluster, appears to be of the order of a magnitude. 
It seems unlikely that rotation can account directly for the whole of this, unless 
rotation were more rapid in the interiors of stars than in their outer regions. 

Table II also shows (a) that the surface oblateness for small angular velocities 
is 0°514%,, Which is close to the figure 0-5, for the Roche model, and to the figures 
obtained by Chandrasekhar (1) for rotating polytropes, and (+) H, xg,“ to the 
first order in %,, as expected from von Zeipel’s theorem. 
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The distribution of p and 7 within the star can be computed to the first order 
in x,. Working in terms of (v,@), where v is the proportion of the total radius 
x,(0, X) of the star at an angle @ to the rotation axis, the substitutions (3.4), expanded 
in the forms defined by (3.2) and (3.3), give the first-order expressions 


p(vs05 X)= pote) [ 2+ 9 Se) a {21°} 4 = Su 2h], 2 


pO) ass Xos < 


tT. tz 1s ty > 
70; N= TI) [1 +4 OS | fat Sts} | 68 


i, A A Se ee ee ee me ieee i 


te ee ewes 
— 











- 
ae oad 


aS 











RPh aS = a ff 
0 0.1 0-2 O38 O04 05 06 0-7 08 09 1-0 
1/1, ——> 
Fic. 2.—First-order density changes. 
--- polar axis. 
equatorial plane. 





where the right-hand sides are evaluated at x=vxy,. ‘The results are presented 
graphically in Figs. 2 and 3 in the form of the percentage changes per unit ~, in the 
corresponding values at the same v in the non-rotating star. It can be seen from 
the figure that the surface density at the equator becomes negative for ~, >0-22, 
corresponding to «>0-29; the model formally breaks down at this value of ~. 
A first-order model will therefore not give accurate results for surface phenomena 
such as figure of surface, variation of effective temperature over the surface and 
surface gravity in the more rapidly rotating stars. In discussing these stars 





No. 6, 1953 The structure of rotating stars 713 


von Zeipel’s relation H, « g,"" therefore cannot be used as in Slettebak’s discussion 
(2) of limb darkening in O and B stars, since it only holds to the first order in «. 

It can be stated in conclusion that (i) the direct effect of rotation, apart from 
questions of mixing, may account for part of the observed spread in luminosity 
for a given spectral type in the case of B and A stars in clusters, and (ii) a theoretical 
study of limb darkening and also of circulation currents in rapidly rotating stars 
requires a second-order model. 
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Fic. 3.—First-order temperature changes. 
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The following table gives the quantities ¢,, ¢; and o, associated with the 
effective potential, temperature and density, respectively, from which the first- 


order model illustrated in Figs. 2 and 3 was derived. 


x 
° 
o'2 
o'-4 
°°6 
o'8 
1‘o 


1 
° 


—0'032—0'026P, 
—0'126—0-102P, 
—0'274—0°224P, 
—o-465—0°386P, 
—0-684—0°'582P, 


ty 
° 
—0'032—0'026P, 
—0'126—o0:102P, 
—0'274—0°224P, 
—o'465—0-386P, 
—0'684—0°582P, 


Go 
° 
—o0:048—0-038P, 
—o'186—o0-150P, 
—0°398—0°326P, 
—0'660—0°549P, 
—0'943—0°802P, 





— 0907 —0°792P, 


—0°907 —0°792P, 


—1°207—1:054P, 





N NNN FS HS 
DARN OC OAS DN 


801P, 
034P, 
-271P, 
-502P, 


o- 
—o'9rb—0-: 

—1°143—1° 

—1°347—1 

—1°513—1 

—1°631—1°722P, 
—1-698—1-928P, 
—1°714—2°122P, 
— 1'683—2°305P, 
—1-609—2°481P, 
—1°498— 2°653P. 
—1°356—2°824P, 
—1°188—2-999P, 
—o'991—3°180P, 
—0'786— 3:369P, 
—0°557—3°560P, 
—0°312—3°780P, 
—0o0'054—4:006P, 
+0:219—4'246P, 
+0:505—4'501P, 
+0:806—4°772P, 
+1'119—5:060P, 
+1°446—5°364P, 
+1°788—5-686P, 
+2°143—6-024P, 
+2:°511—6-380P, 
+2:°895—6°752P, 
+3°293—7°141P, 
+3°705 —7°547P2 
+4°132—7:°969P, 
+4°574—8-406P, 


—o-g10—0-796P, 
—0'999—0°907P, 
—1:068—1-016P, 
—1°104—1°119P, 
—1°122—1'216P, 
—1‘192—1°308P, 
—1:078—1°394P, 
—1-019—1°478P, 
—0'938—1°560P, 
—o'840—1°643P, 
—0o'724—1°729P, 
—o'597—1'818P, 
—0°454—1°914P, 
—o°308—2:016P, 
—0o'149—2°126P, 
+0:017—2:'244P, 
+0:'192—2°372P, 
+0°374—2°509P, 
+0:'561—2:°656P, 
+0°755—2°813P, 
+0°959—2-980P, 
1°168—3°156P, 
~1°385 —3°345P2 
+ 1°608—3°544P2 
+ 1°837—3°753P2 
+ 2°073—3°972P, 
+2°315—4:200P, 
+2'563—4°430P, 
+ 2°818—4-687P, 
+3°079—4°945P2 


—1°421—1-008P, 
—1°554—1°403P, 
1°764—1-°658P, 
-831—1'808P, 
-764—1°852P, 
*599— 1°806P, 
-377—1°696P, 
1°132—1°545P, 
0'894—1-376P, 
-0°679—1°203P, 
—0°495 — 1°037P, 
—0°345 —0'884P, 
—0'226—0-747P, 
—o-'138—0°625P, 
—0'073—0°519P, 
—0:027—0°427P, 
+0:004—0°349P, 
+0:022—0:281P, 
+ 0:033 —0°224P, 
+0:036—0:'176P, 
-+0°036—0°'135P, 
+ 0°'032—o'101P, 
+0°027—0:'073P, 
+0°021—o0'051P, 
+0:015—0'033P, 
+0-010—0'020P, 
+0:006—o-o10P, 
+0°002—0'004P, 
—o-ooiP, 





Xog=7°01168 


+ 4*600—8-432P, 


+ 3°094—4°960P, 
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ROTATION AND STELLAR EVOLUTION 


L. Mestel 
(Received 1953 September 2) 


Summary 


Sweet’s discussion (1) of rotational currents in a Cowling model star 
is reviewed. It is shown that the non-spherical distribution of matter set up 
by the currents themselves tends to choke back the motion, so that no 
rotationally stable star can suffer any considerable mixing between core 
and envelope, if the angular velocity is nearly constant throughout the star. 
If, however, the rotation rate increases substantially inwards, enough mixing 
can occur to influence seriously the evolution of the star. 





1. Introduction.—Recent research on stellar structure has shown that the 
observable parameters of a star depend strongly on the composition of the 
stellar material, and on inhomogeneities in its distribution. Thus the 
mass-luminosity and mass-radius relations for homogeneous stars (2) show that 
the luminosity increases strongly and the radius moderately as the mean 
molecular weight increases, leading to a slight increase in effective temperature. 
Further, it is found that the radii of some models with inner zones of high 
molecular weight are large enough to account for the most diffuse giant stars 
observed (3, 4,5). Itis also known that certain non-homogeneous stellar models 
depart from the static state (4, 6); the evolution of such stars is likely to differ 
profoundly from that of static stars. 

Many theorists believe that it is impossible to explain the evolution of the 
stars without invoking some interaction between the stars themselves and the 
interstellar medium. Thus Hoyle and Lyttleton (7, 8) believe in accretion of 
interstellar matter, both on observational and theoretical grounds; on the 
other hand, Struve (9), Fessenkov, Ambartsumian (10) and others propose 
equally strongly that stars emit “‘ corpuscular radiation ’’, so losing substantial 
quantities of their mass at a rate comparable with their rate of evolution by 
consumption of hydrogen. Whatever the arguments for or against either of 
these mechanisms, it seems desirable to work out carefully the course of 
evolution of a star when left to itself. When this has been done, appeal to 
observation may be able to decide unambiguously what further hypotheses are 
necessary to avoid discrepancies between theory and observation. 

A star of mass well above that of the Sun generates its energy by means of 
the Bethe cycle; as shown by Cowling (11), the high temperature-sensitivity 
of the cycle causes the innermost regions to be in convective equilibrium. Such 
a star has had sufficient time since the birth of the Galaxy to undergo considerable 
evolution through the conversion of hydrogen into helium in the convective core. 
The convection currents are sufficiently rapid to keep the core uniform in 
composition, but the radiative zone is highly stable against turbulence, and 
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other causes of mixing must be sought if material of high molecular weight is 
to penetrate into the radiative zone. 

In a recent paper (x) Sweet has discussed the suggestion that Eddington’s 
meridional currents, due to stellar rotation, will provide the necessary stirring 
in a short enough time. Sweet’s corrections reduce Eddington’s value for the 
vertical velocity of the currents by a large factor; in particular, he shows that 
mixing between core and envelope in the Sun (assumed to rotate throughout 
with its surface angular velocity) must be negligible. Early-type stars, however, 
tend to have much higher surface angular velocities, and Sweet’s figures lead 
him to conclude that such stars, rotating as rigid bodies, will remain homogeneous 
in structure, except for the effect of accretion. 

Except for one detail, discussed in Section 2 and the Appendix, this paper 
does not criticize Sweet’s analysis as far as it goes, but it will be shown that the 
conclusions he draws as regards stellar evolution need considerable revision. 
The Eddington currents are due fundamentally to the departure from spherical 
symmetry caused by centrifugal force. It is impossible, in a uniformly rotating 
star, to satisfy both the hydrostatic condition and the equation of radiative 
equilibrium ; heating occurs at the poles and cooling at the equator of each level 
surface, and a system of meridional currents is set up, with velocities just 
sufficient to preserve thermal equilibrium through the energy they transport. 
Similar results hold for a star with a more general angular velocity field. 

The inadequacy of Sweet’s analysis arises from the neglect of the departure 
from spherical symmetry brought about by the currents themselves. Material 
is transported out of the core at different epochs, and the value, at time t, of the 
mean molecular weight » at a particular point of the radiative envelope is the 
value of ; in the core at an earlier time (t—#), where ¢ is the time of flow from the 
core to the point in question. As ¢ is a function of position, in general a 
non-spherical y-distribution will be set up, which will react back strongly on 
the thermal equation, and hence on the forces driving the circulation. It is the 
object of this paper to combine with Sweet’s theory a discussion of this 
“feed-back” process. 

The mechanism of mixing by rotational currents is very delicate, and the 
detailed theory discussed below is relevant only on the assumption that the 
radiative envelope does not suffer any local turbulent mixing, tending to smooth 
out the distribution of matter set up by the currents. The strongly subadiabatic 
temperature gradient in the envelope is stable even in a star with considerable 
non-uniformities of angular velocity (12), and most workers assume that no 
convection occurs within the radiative zone. An exception is Alfvén (13), who 
suggests that magneto-hydrodynamic waves emanating from the core may bring 
about local convection even in the stable envelope. However, he brings forward 
no quantitative arguments, and the possibility will be disregarded in this paper. 

In Section 2, Sweet’s theory, in which centrifugal force is the only 
perturbation, is summarized, and in the Appendix it is extended so as to avoid 
the singularities that appear at the surface of the star, near surfaces of 
discontinuity within it, and near convective regions. In Sections 3, 4 and 5, 
the effect on the radiative equation of a non-spherical distribution of matter is 
determined, and the results applied to a rotating star possessing initially a 
gravitationally stable »-gradient. The rest of the paper is devoted to an analysis 
of the evolution of a rotating Cowling model star. 




















718 L. Mestel Vol. 113 

















2. Summary of Sweet’s theory.—We consider a Cowling model star rotating 
as a rigid body with angular velocity (2; the theory is easily extensible to other 
cases. The equations of mechanical equilibrium, and Poisson’s equation, are 


oP od — 
zz 7% + prQQ? sin? 4, (1) 
oP 0 : 

a ee 2¢)2 2 
30 =? 30 + pr?Q? sin 6 cos 0, (2) 


where P, p and ¢ are the perturbed pressure, density and gravitational potential, 
and r,@ are respectively central distance and colatitude. Each function is 
expanded as a series in the perturbing parameter; to the first order, 
P= P,(r)+ P,(r, 6) . Q3, 

P=PpolT) + pilt, 9) . Q-, (4) 
f=4(7) + 4,(7, 8) . Q?. 

Similarly, the temperature 7 is written as 
T=T,(r)+ T,(r, 8) . Q?. (5) 
P is eliminated from (1) and (2), and dp,/00 is substituted from (3). 
x = 0¢,/00 is then given by 
x =1r(A—1) sin @ cos 6, (6) 
where h satisfies a homogeneous second-order linear equation. y must vanish 
both at the star’s centre and at infinity, while both y and @y ér are continuous 
at the surface. In the present problem, 6?y/ér? is also continuous there because 
of the continuity of density, in contrast to the problem of a gravitating liquid. 
From (3) and (2), and the equation of state, dp,/00 and 07,00 are now found, 
and substituted in the thermal equation 
, v, dS ; 
Pot o 1 eo — div H. (7) 
Here S, is the entropy per gram-mole in the unperturbed state, v, the vertical 
component of the circulation velocity and H the radiation flux; only first-order 
terms have been retained. This equation asserts that the rotation causes at each 
point a net influx of energy of amount (—div H), and that the rate of circulation 
must be such that just this amount of energy is convected away. 
‘Together with the condition of zero outflow over a sphere, 


le, sin 0 dd =0, (8) 
0 


v, is given by (6) as 
v,=f(r)P,(cos 0)Q?(LR*/M?), (9) 

where f(r) is a complicated function that must be tabulated*, and L, R, M are 
luminosity, radius and mass respectively. The horizontal component velocity 
is given by the equation of continuity. 
The perturbing terms ¢4,, p,, etc., contain also parts independent of @. ‘Thus, 

for example, 
$,= —4r(h—1)P.+h,(r). (10) 


* I am indebted to Dr Sweet for informing me that the “‘ 3’ in equation 3.8 of his paper (1) 
should be a “‘ 4’’. 
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The function k, must vanish at infinity, and again, together with its first two 
derivatives, is continuous at the edge of the star. Also, k,’ and hence k, are zero 
outside the star; this follows from the constancy of the stellar mass under the 
perturbation. The parts in (1) and (3) independent of @ give equations to the 
perturbed mean state. In the radiative zone, the part of (div H) dependent on 0 
causes the circulation; when this is averaged out, as in (8), we are left with the 
perturbed mean state energy equation. ‘Together with the perturbed adiabatic 
condition for the convective core, these equations suffice to determine the 
change in the mean state brought about by the rotation. ‘To the first order, 
the equations giving circulation can be separated from those giving the perturbed 
mean state; in a higher-order theory, however, coupling necessarily arises. 

Sweet points out that the difficulties found by previous workers on the 
problem result from demanding the too stringent condition that the vertical 
circulation velocity should vanish at the star’s surface. He argues that it is 
sufficient that it be finite at the surface, as the zero density implies that the 
vertical transport of matter vanishes. But this cannot be the final answer, for 
the equation of continuity then demands an infinite horizontal component at 
the surface; in order that the horizontal component be finite, the vertical 
component must vanish as least as rapidly as the depth beneath the surface. 
A theory which gives an infinity must be incomplete in some detail. 

A similar difficulty arises in a star with a discontinuity in , such as is assumed 
to exist in current giant models, with the conditions of continuity of pressure 
and temperature being preserved by a jump in density. Such a surface cannot 
be distorted without violating the hydrostatic equation, and so bringing into 
play restoring gravitational forces. ‘hus meridional currents cannot cross from 
one zone into the other, but are deflected horizontally. But the heating and 
cooling due to rotation must still occur, and at an infinitesimal distance from the 
discontinuity a finite vertical velocity field exists. ‘The only way of fitting into 
Sweet’s quasi-static theory a zero vertical velocity at the barrier is to let the 
horizontal component become infinite as the barrier is approached from either 
side. Physically, this means that the energy available to drive the circulation is 
not used up, as the motion is constrained by the barrier to be horizontal. There is, 
to the first order, no horizontal entropy gradient, and so this energy accelerates 
the material. ‘The only steady state allowed by Sweet’s theory unmodified is 
one with an infinite horizontal component along the barrier. 

In fact, the material cannot be accelerated indefinitely ; the viscosity of stellar 
material is finite, and sooner or later the velocity will tend to a finite limit. In the 
ultimate steady state, the viscosity terms in the equations of motion so modify 
the temperature field that heating and cooling over the barrier is eliminated. 
The theory of this process is given in the Appendix. It should be noted that 
this effect is due to viscous force, which depends linearly on the velocity 
components, and not on the effect of viscous dissipation on the energy equation, 
which is quadratic in the perturbation, and can therefore contribute only in a 
second-order theory. For the same reason, the non-linear inertial terms in the 
equations of motion are of no help; in the expansion of the circulation velocity 
in powers of {2?, each coefficient of the vertical component must vanish at the 
barrier. 

Viscosity is small enough to have negligible effect on the circulation at 
non-singular points; it is only when the velocity gradient becomes large that it 
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plays a decisive role. Near a convective zone the polytropic index m approaches 
the adiabatic value $. In Sweet’s theory the vertical component « 1/(m— 3) 
and the horizontal component « 1/(m — #)?, and so again a large horizontal shearing 
motion is set up, and the viscosity prevents the velocities from becoming infinite. 
In all cases the layer in which viscosity acts is very narrow, and, as in hydrodynamic 
flow theory, will be termed a “‘ boundary layer”. 

3. The effect of a non-spherical y-distribution on the radiative equation.— 
Consider now the modification of the rotational currents brought about by the 
unsymmetric nature of the ,-distribution to which the rotational currents 
themselves lead. The distribution in »=,(r, @) may be regarded as a perturbation 
from a spherically symmetric distribution p,(r). A convenient parameter 
measuring the perturbation is 


Kr, 9) — Hol?) = Xr, 0). (11) 
Bo(7) 
Here jx(7) can be defined simply as the average of y(r, @) over a sphere r = const. ; 
or it could be an original spherical distribution perturbed by rotational currents 
to »(r, @). The only condition imposed is that A should be small, so that 
first-order perturbation theory is applicable. Then to the functions py, Py, bo, Ty, 
there will be added perturbing terms due both to rotation and to the perturbation 
in «. Since only first-order terms in the two perturbations are retained, the 
two effects may be simply superposed. Hence, in considering the immediate 
effects of the perturbation A(r, #), we shall regard the star as non-rotating. 
Let p,, P,, 6, and 7, respectively be the perturbing effects of A. Then the 
equation of hydrostatic support and Poisson’s equation 


grad P=(grad 4)/p, (12) 
V°d= —47Gp (13) 
show that p,, P,, ¢, are spherically symmetric. From the equation of state 
R 
P=— pT I 
ta (14) 
it follows that 7/ can be written as 
{To(r) + T(r)}/H0(7)s (15) 
leading to ¥ 
T(r, 9) =A(r, 0)To(r) + T(r). (16) 


The function T is the mean of 7;(r, 6) over a sphere only in the case when the 
mean of A vanishes. 
The first-order effect of the perturbation on the thermal equation is now 
required. We have 
—-divH 1 0 /rT*0T s 277°. oF 
a 55 (FS) + a(S sine Z). (17) 
$ac r-or\ xp or r? sin 6 00\ Kp 00 
where « is the opacity. It is supposed that all the energy liberation occurs 
within the convective core, where the present theory does not apply; then the 
unperturbed equation of heat flow is 


dT, oPoL. I 


oe 


7" @ 4nr? * $acT,3’ 


(18) 
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where L is the luminosity. Also, « is taken as a function of p and 7 independent 
of , as is approximately the case. ‘Thus substitution of the perturbed functions 
in (17) yields, after simplifying by (18), 


divH 1 f0fa tT \. 6 T° 
a. 4nrt BE (In =) Pe Apo (In =) (or) 


an T1 1 Of. aa 
* = (T+an'+ $F 1.) } + 73 aaa sine) | (19) 


The spherically symmetric part of this expression corresponds to an 
expansion or contraction as a whole of the sphere considered, and is of no interest 
for the problem of meridional circulation. ‘To the present order of working, 
the spherical part of a perturbation alters the structure of the star, while the 
rest gives rise to circulation, as discussed in the rotational case. ‘Thus (cf. Sweet’s 
equation (3.6)) the perturbation in y leads to a motion with a vertical component 


v*, where 


(n— 3) Lfofoa Ts ur ) OA 1 
ni 2 ee ee) a 1 
"(m+)" ons 477° Lean ‘“ KoPo amen , r} J 


Trt af a 
+ Fe sin 6 96 \5I9 9 5G) |» (2° 


£o= —% is the local value of gravity in the unperturbed star and n is the local 
value of the polytropic index. The horizontal component 7," is then found from 
the equation of continuity 


: a ao ae 
uy sr pr |e sin 6 dé]. (21) 


The term dp/dt should strictly appear in the equation of continuity; however, 
if the changes in p are due ultimately to transport of matter by currents, this 
term will involve v . grad », where v is the current velocity. As v is of order A, 
this term is of second order, and so is neglected. 

In the general case currents arise both from the non-spherical ,.-distribution 
and from the rotation. The total circulation is Sound by adding vectorially the 
velocity v“ just found to Eddington’s velocity v® 

4. Variation of » with time.—The equation of convection in the radiative 
zone 1s 








du oO 
- = = 5 + (v" +v"). grad p=0 (22) 


It is to be noted that the theory ceases to be linear at this point, for the changes 
in » are brought about by the transport of variations in by currents which 
themselves depend on the variations in yi. 

A is now expanded in a series of Legendre functions: 


A=a)+a,P,+a,P,+..., (23 


the odd terms being omitted, as we are mainly interested in perturbations 
symmetric about the equator as well as about the axis of rotation. 

Hf, as is assumed, {2 is uniform throughout the star, the vertical component 
of v® contains only a P, term. However, v,” in general involves an infinite 
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number of Legendre functions. In a star with opacity obeying the unmodified 
Kramers law, 7," is given by 


(n+1) S. L 6a, T> 15 Jat h.; 20% 
oy TP ia, — =» ao — — — 
Fr (n—3)° Pogo 4ur P, - we 2 * Ty a2 T,’ a2 


20a, T 15 Ty \' of Oo 
nae ((B) aH}. 


with corresponding terms in vg“ involving (P3,;— P2,,4)/sin 0. 
The resultant vertical velocity function is written as 


U,=foPotfyPyt “— (25) 
In the present case f,(7) is the sum of an {-velocity and a y-velocity, while the 
other terms are just the respective parts of (24). Then, from (22), the rate of 
change of the P, component of A is given by 


(pr? )H19@q = -{ = (pr*Wode.fos) (48 + )} . (26) 


‘The equations for dy, ay, etc. are much more complicated. It is reasonable, 
as a first attempt, to neglect all components but a, and a,. The equations 
governing the variation in A throughout the motion are then 





5H4oGo(pr®) = — {(pr*)forp@2}’s (27) 
Ho @2 = ~{f, ee + aq))' + F(Ugae)’ ] +3 {2H0) at) pr*f.) yt. (28) 


However, whereas (27) may be expected to give calle qualitatively correct 
and of the right order of magnitude, use of (28) may in some cases lead to 
misleading results. ‘The last term in the brackets in this equation is proportional 
to —vw,, and indicates that a positive horizontal component of velocity tends to 
increase a», i.e. to increase the variation of 4 over a sphere r=const. Simple 
physical considerations show, however, that a sufficiently great horizontal velocity, 
whatever its sign, tends to reduce aj. This indicates the importance of the terms 
neglected in (28); in fact, in Section 10, this equation is shown to imply an 
incorrect sign for a,’ in part of the star. However, the principal results may be 
obtained by use of (27) alone. 

5. Direction and magnitude of the y-currents.—The first question of interest 
is whether the ,.-currents choke back or augment the Eddington circulation. 
The function f, _— the vertical motion can, by (24), be written as 


(n+1) pe ss 6 Ty d 1 , 15 | b 
f{r)= =Pr)+ 3 r? (n—3)° Poo r2 = ~~-z ; — Pay’ } (29) 


where f(r), proportional to 2?, is given by Sweet’s theory, and the remaining 
terms give the u-currents. We have to determine the sign of the a,-term in this. 

Now py <0, since the »-gradient is due to the progressive transmutation of 
hydrogen into helium; a positive j9-gradient is in any case unstable, as follows 
from a simple application of Archimedes’ Principle to a small vertical disturbance. 
The a,-distribution is due to the distortion of the surfaces of constant jz, by the 
circulation, bringing heavy material up near the axis of rotation; thus a, >0. 
Accordingly, since 7,’ <o, the first term in the brackets in (29) is negative, and 
its effect is to choke back the rotational currents. The other terms in the brackets 
involve a,’ ; thus, if a, does not vary too greatly from point to point, the effect of 
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the -currents is to reduce the circulation. The term —‘2a,' increases the 
choke-back effect if a, increases outwards; the remaining term is smaller unless 
a,’ alters very rapidly. 

These conclusions can be put in another way. Suppose that locally a,~7)"; 
then the a,-terms in (29) produce a reduction in the circulation if 


. _6/f T, \3 
1b) < — 0 : 
(7 + 2 < r2 ( Jo ) ° (30) 


The expression on the right has been calculated for a Cowling model star, using 
Miss Gardiner’s tables for the model (14). The range within which » must 
lie for a reduction in the circulation is shown in Table I. 


TABLE I 
Range of 7 
—8°32 <n < 0°82 
—7°82 <n < 0°32 
—7°65 SN S0'153 
—7°57 <n <0°066 
7 (edge of star) —7'50<n<0 


mupPwWNnN A 


‘Thus a reduction in circulation is secured if a, increases at any reasonable rate, 
or decreases fairly slowly, as r increases. 

As an example, consider a star which possesses initially a negative .)-gradient 
(without inquiring as to how the gradient has been built up). ‘Then except near 
a convective or other singular region both vertical and horizontal components 
of velocity are of the same order; hence if |» | is not too small, equation (28) 
reduces to 


Hos = — Ho fay (31) 
where /,(7) is given by (29). Hence, after time f, 
A= —(Hg /Mo) - P(r). t (32) 


approximately. ‘The function p(r) is slowly varying, and so if 1,’ does not change 
too rapidly with distance, the distortion of the surfaces of constant jy, by the 
()-currents induces p-currents which oppose further distortion; if |j9’| is large 
enough, the motion may be damped out altogether. Equation (29) shows that 
the ,.-currents reduce the circulation very considerably if 


a~p). Tad (- 3). 
OL ° (n+) poo To 
On substituting numerical values, this gives a,~1°5 x 10-°(Q)?, where Q is the 
angular velocity in solar units. (The solar angular velocity is taken as 
2°85 x 10 ® radians sec.) With normal angular velocities this is quite small; it 
corresponds to a distortion of the surfaces ~,=constant through a distance 
AP 34:9 | Ho |, which is also small if | j49'| is at all appreciable. This is sufficient 
to show how profound is the effect of asymmetries in « on the thermal equation. 
The jo-distribution through the star may be such as to damp out the 
rotational currents through a part of the star only. Thus, for example, consider 
a star with a j9-gradient which is both large and slowly varying throughout a 
large zone, but which gradually tapers off to zero on both sides. The circulation 
in the zone of varying jy rapidly distorts the surfaces of constant jo, and leads 
to an a,-distribution which, like the »-gradient, is large and slowly varying. 


(33) 
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Such an a,-distribution will be able effectively to stop the circulation in the zone 
of varying #49; but away from this zone the p-velocities become progressively 
weaker, until the (-currents cease to be inhibited at all and Sweet’s theory 
applies. Hence there must be critical radii between which the circulation 
effectively stops, while outside them the Q-currents sweep away the material 
with small initial j»-gradient. Material brought along by the {-currents in 
the zones outside cannot flow past the barriers at the critical radii, and so must 
be deflected horizontally. ‘Thus outside the barriers Sweet’s currents circulate, 
with boundary layers reducing the vertical velocities to zero at the barriers; the 
star settles down into a quasi-static state, with the region of high p,.-gradient 
separating the two zones of circulation. 

The discussion based on (31) breaks down near the convective core, where 
the horizontal component of velocity is large compared with the vertical, and 
where material of high , is flowing in from the core. However, some information 
as to the a,-distribution set up near the core may be gained by integrating (27) 
from a radius r to the surface R: 


R 
pratt =5 noitaprtdr. (34) 
rT 


This equation may be interpreted as indicating that the time-variation of yu 
outside the sphere of radius 7 is due to the radial flow, proportional to /,, carrying 
heavy material across this sphere. Near the convective core the Sweet velocity 
is large, and f, can be expected also to be large; thus a,=0 at the boundary of 
the core, and increases outwards. Thus unless a,” is very large and positive, 
the y-currents set up by the material flowing out of the core may be expected 
to impede the flow near the core. 

The results of this section and of the Appendix show how powerful is the 
influence of previously existing non-uniformities of ,« on the meridional 
circulation ; a naive application of the results found by Sweet and Opik is liable 
to be misleading. Thus, for example, in a recent paper (15) the present writer 
has considered the distribution of accreted matter through a white dwarf. It was 
shown that the simple Eddington currents are negligibly slow, and it was 
concluded, correctly, that mixing is therefore negligible; but, in fact, it is now 
seen that the circulation must be divided into two separate zones by the 
discontinuity in composition. 

6. Rotational currents in a Cowling model star; the equations of circulation.— 
The most important application of the theory of this paper is to the ;.-distribution 
set up by rotational currents within a Cowling model star. In the convective 
core « increases with time; the rotational currents drag material of high ,« into 
the radiative zone, and also slow down the increase of « within the core by 
replenishing it with hydrogen-rich material. But as explained above, the 
inevitable asymmetry of the p-distribution must affect the driving velocity, so 
that both qualitatively and quantitatively the resulting circulation will differ 
from that predicted by Sweet. 

In this paper the star is at first taken as rotating uniformly; later, the 
possibility is mentioned that the star may consist of a uniformly rotating bulk, 
with outer layers rotating considerably slower. In either case the angular 
velocity field is taken as constant in time, so that the rotational current velocities 
also remain constant in time. The conditions for this to hold are discussed in 
Section II. 
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The equations to the circulation are (27) and (28), with the unperturbed 
mean molecular weight p, a constant. These equations then become 
: i. ae ; 
54> — Ade (pr*fa(7)a2), (35) 
, ea | a 
d= — 4 fol r)(dy + Fag’) +7 pr dr (pr*f2(7)) Ps (36) 
where /,(r) is defined in (29). Equations (35) and (36) can be combined to give, 
as an alternative to (36), 
dy— Fay = —fo(r)(Qy + 742’). (37) 
In order to isolate the driving parameter {2?, and to simplify the analysis, 
the following transformations are made. ‘The non-dimensional radial coordinate 
used by Cowling (11) is introduced by the substitution 
r=€3, (38) 
2 =(5RT. SzuGp,). (39) 
p, and 7. are the unperturbed central density and temperature respectively. 
The quantities y, w and m are defined by 


where 


P=PcJs 
T=T.w, (40) 
M(r) = 47p,<3m. 
The driving velocity p(r) is then written as 
p(z) . (L Q2 R5/M3), (41) 


where quantities with bars have been expressed in solar units, and f(z) is a 
tabulated function of order 10 !° cm/sec through most of the star. In addition, 
a and a, are put equal to yA, and yA, respectively, where 


_ [107*p,?Get LQ2R5 ” oi O2R3 
y= ( Tol ).( WE )=727 x10 ( uM ). (42) 


Then (35) and (37) reduce to 





5A A= — —, {8(F + CA, + DA,’ + EAQ")}', (43) 


V2* 


(A,—-#4,A= — S (F + CA,+ DAy! + EAg")( Ag’ + 44s’), (44) 


n+1) 
ih a (45) 


p= 


, MB. 10 ah wi carta 
5A = 5€ SS = 4:00 x I 20 3 RIEL Q?, ¢ 
fos yee «eee (46) 


F=10"p(z) . (n—3)ym (n+ 1), 
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The factor 10! in F, y and A is chosen so that F, C, D, E are all of order unity 
through the bulk of the star. ‘The expression (F+CA,+ DA,’ +EA,") 
proportional to the radial part of the vertical component of the velocity. 

It is assumed explicitly through the following work that the stellar functions 
y, m, etc., and hence also the circulation functions F, C, D, E, remain unchanged 
functions of z as changes in » are brought about by the circulation. This is 
certainly only an approximation, as the altered structure of the star will react 
back on the driving velocity. However, it is shown below that if circulation 
is able to persist at all it will tend to bring the star into the Hoyle—Lyttleton giant 
state, which can be looked upon as effectively a Cowling model star with a 
highly distended envelope superimposed. Hence the circulation in the inner 
zone, which is our main concern, will not alter sensibly on account of changes in 
the structure of the star. 

7. Quasi-static circulation.—The (-current streamlines calculated by Sweet 
traverse the whole star from the convective core to the surface. If the feed-back 
terms are neglected, the variations in through the star are due simply to the 
different times of flow from the core along the (-current streamlines : 


H(z, 6) =p,,(t ie i(z, 6)) (48) 


for t>t. Here the suffix c refers to the core and ¢(z, 6) is the time for matter to 
reach (zs, @) from the core along the relevant streamline. The luminosity changes 
only slowly with time, so du,/dt is a constant. Provided, then, that ¢ is finite 
for all (z, #), once matter has had time to circulate once round the star, du /dt 
will be constant everywhere. ‘This implies that A, is a constant, and A, etc., are 
zero. Such a circulation will be called “‘ quasi-static” 

On Sweet’s theory, a star will remain ‘“‘ homogeneous”’ if the variations in sj: 
brought about by flow along the {2-current streamlines are not large enough to 
have any sensible effect on the structure of the star; otherwise it will become 
‘“‘non-homogeneous”’. The present theory shows, however, that the streamlines 
themselves may be profoundly modified by the distribution of matter, and it 
may not be assumed a priori that matter from the core will reach the edge of the 
star even if sufficient time be allowed. It is the principal aim of the rest of this 
paper to estimate as a function of {2 the proportion of the stellar mass through 
which currents from the core may circulate continuously. 

Consider now the flow due to the joint effect of Q-currents and ,-currents. 
Equation (48) still holds, but #(z, 6) must be measured along the joint stream- 
lines; and as p(z, 0) changes, so will ¢ in general. We are interested, however, 
not in the details of the initial motion but in the ultimate state to which the 
circulation approximates. One possible assumption is that the flow settles down 
into a quasi-static state as defined above, except that the zone of circulation is 
not assumed to extend throughout the star, but only up to a radius 2’. This 
implies again that within this zone A, is a constant, and A,, Ay, etc., zero; as 
the driving velocity depends only on A,, etc., such assumptions are self-consistent. 
The interior of the sphere of radius 2’ is called the “‘ mixing zone”. The further 
assumption is made implicitly that ¢ is everywhere finite within the mixing zone; 
this will be reviewed later. 

Since we are interested less in the details of the circulation than in the effect 
of the currents on the structure of the star, the quasi-static case is the one best 
worth considering. For suppose that after a certain radius 2, A, drops off 
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sensibly from a uniform value. Relation (48) must still hold for a region connected 
with the core within a finite time, and so it follows that the average of dt/dt over 
a sphere must be positive after z. Thus the time taken by matter to reach points 
beyond 2 increases steadily, and this can only mean that the currents themselves 
build up such a p-distribution as to cause them to decay after the radius z. The 
value of 4, beyond z will therefore decrease still further with time, so that the 
decay progressively increases. Hence it is only in regions of quasi-static 
circulation that rotational currents can be expected to have an important 
influence on the evolution of the star. , 
Equations (43) and (44), with A, constant and A, zero, give 


BA, F + CA, + DAy’ + EAg"| =(5A9A)(m' — m), (49) 


A,y2* 
i © in» sqummganaite 50 
(4, (m’ ee m) ) bad (5: ) 


The constant of integration m’ is the non-dimensional mass of the mixing zone; 
the equation of mixing, (49) states that across any sphere of radius z within the 
zone the net flow of molecular weight is sufficient to increase A, at the rate 4, 
within a mass (m’ —m). , 

The first problem is to determine the relation between m’ and (5.)A), i.e. to 
find for what values of the spin quasi-static flow is possible within a sphere of 
mass m’. Before (49) can be solved, however, the boundary conditions must be 
determined. At the surface of the core A, must vanish, as the horizontal 
component of the velocity becomes large there compared with the vertical. This 
condition is consistent with (49), as the term 1/8 0c(m—$) and so vanishes at the 
core surface. The condition at 2’, however, requires more thought. 

Consider first values of m’ less than the stellar mass. Then at 2’ either 4,=0 
or (F+CA,+DA,'+EA,”")=0. (It is easy to verify that solutions behaving 
near 2’ like (z’— 2)" with o<7<1, do not exist.) Now imagine a barrier placed 
at s’ so as to prevent flow beyond this sphere—either a discontinuity in yx at 2’, 
or a region beyond 2’ in which an A,-distribution cancels out the Q-currents. 
In a star with infinite spin the ,.-current terms disappear, and there is no question 
of any choking of the circulation; the currents flow up to the barrier and are 
deflected there. The vertical velocity component remains finite at the barrier, 
and so from (49) A, vanishes there. 

Now let the spin of the star be gradually reduced, so that the parameter 
(54,) increases. The feed-back terms begin to modify the circulation, but 
initially there will be solutions of (49) satisfying the conditions of vanishing 4, 
at each end of the range. Sooner or later, however, the equation will cease 
to have a solution satisfying these conditions, implying that 4) cannot remain 
constant throughout the zone. It is this limiting solution in which we are 
particularly interested. The physically relevant solutions of (49), for given m’, 
must form a continuous sequence in (5A,A), and so the condition A,=o is the 
correct one to use. Near s=2’, A, behaves as k(z’ — z), and it is possible, a priori, 
that the constant of proportionality k may tend to infinity as the limiting value of 
(5A,A) is approached, so that it is in fact the vertical velocity and not A, that 
vanishes at 2’ in the limiting solution ; but if this is so, the solution of the equations 
will produce this result, and no a priori assumptions about the velocity at 2’ are 
needed 
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However, it must be remembered that the problem just summarized is a 
problem with a constraint, namely a barrier at 2’. A Cowling model star has no 
barriers except those built up by the circulation itself, and so interpretation of 
the results will need care. This is further discussed in Section 9. 

If m’ is the mass of the whole star, the same argument applies, and A, =o 
is again the correct boundary condition to use. In this case, of course, there is 
no need to introduce an artificial barrier. From (47) and (49), it follows that 
A,=k(z' — 2)?" near the edge of the star—the feed-back currents completely 
dominate the (2-currents, and give rise to an infinite positive vertical component 

of velocity. ‘The constant k remains finite for all finite values of (5 A,A), and so 
in this case we may anticipate that as the limiting value of (5.,\) is approached, 
the velocity remains non-zero (and, in fact, infinite). 

It has already been remarked in Section 2 that any infinities given by the 
theory can and must be eliminated by the action of viscosity within a boundary 
layer. However, the zone in which the simple theory breaks down is very narrow, 
and this correction is negligible when considering the properties of the 
circulation as a whole (cf. Section 10). 

The parameter characterizing the particular quasi-static circulation is (5. 4d), 
which by (42) and (46) can be expressed as 


10 (Z) | me (51) 


where M and R are mass and radius in solar units and q is the ratio of centrifugal 
to gravitational force at the equator on the surface. Use has been made here 
of the fact that in the conversion of hydrogen to helium, ;; of the mass is 
available for liberation as energy. The quantity m’ appears explicitly because 
the rate of increase of « in the core depends on the mass of material circulating, 
that is, on the mass of the mixing zone, and not just on the mass of the convective 
core. The angular velocity appears raised to the fourth power. ‘This is due to 
the non-linearity of the problem; for if 22 is doubled, for example, not only is 
the driving velocity doubled, but the feed-back terms are a as only half the 
time has been available for changes in « to take place in the Core. 

It is seen from (51) that the theory is not quite homology-invariant, in that 
the value of q giving the same amount of mixing varies slightly with mass: with 





RxM®®, q varies as M®!), ‘This is because the driving velocity, which is 
proportional (16) to 273 GM(r), is not itself homology-invariant. 

According to the theory without j:-currents, the efficiency of mixing between 
core and envelope is measured simply by the velocity of the (-currents, which 
by (41) is proportional to 2?R°L/M%. From this it follows that if a star in the 
course of its evolution changes its dimensions at all, the mixing efficiency remains 
effectively the same; for (2?R* must remain nearly constant, by conservation of 
angular momentum, and so the velocity is proportional to the radius, implying 
that the time of flow of the {-currents is unaltered. On the present theory, 
however, the efficiency of mixing is measured by the size of the parameter 
(5A4,A), each value of which is shown to define a mixing zone. ‘This parameter 

varies as RM? QR, which in a star of given angular momentum varies as R. 
Thus if the star contracts during the course of its evolution, (5: 4,A) drops, and the 
mixing efficiency of the currents increases. The physical reason is that, from (24), 
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v, xJR*, whereas v2 oR; thus as the star contracts, the feed-back terms decrease 
more rapidly than the driving velocity, so that mixing becomes easier. 

8. Solutions of the equation of mixing.—The simplest case to discuss is the 
limiting one with the mixing zone coinciding with the convective core, and so 
containing 0-15 of the stellar mass. It is found convenient to introduce here 
the function x defined by 


xF = F+CA,+ DA,’ + EA,". (52) 
v is just (v2 +¢7,")/v,2, and so measures the effect of the .-current terms on the 


current velocity. 
Then from (49) and the two boundary conditions, 








i (5A,A) ss m ys? - 
a i Wie ila dz*(n+1)3?° F J’ (53) 
where =’ coincides with z,, the radius of the core. Then 4, and 4,’ both vanish, 
and 
Nien ase dn m y (5.4,A) ‘ 
A dz" (nt+1)'F° x * (54) 


From (52), 
2En'my (5A,\) (55) 
(n+1)F° x ’ " 


SEn'my : 1/2 : 
x=} [ a {1 - (Gar sao} | (56) 


where values at z, are used for the quantities in parentheses. ‘There is no 
ambiguity of sign in the solution, for as 2+ «0, x must approach unity, as the 
p-currents then become negligible; to ensure this, the positive sign must be 
taken in (56). ‘Thus (56) provides an immediate discriminant for material to be 
able to escape from the core. It is to be noted that in the limiting case the vertical 
component of velocity is reduced not to zero but to just half its rotational value. 

Inserting numerical values in (56), it is found that (5.4,A)=0-05733. Taking 
as a typical case a star with M=5, (51) gives that g=0-96, implying that on the 
surface the equatorial value of centrifugal force nearly equals gravity. 

This rather startling result merits some comment. ‘The velocity of the 
rotational currents at any point of the star is proportional (16) to (2?r? GM(r), 
and with £2? constant this decreases by a factor of about 30 from the surface to 
the core. But at the surface Q?R*/GM(R) =q is limited to be <1 by the condition 
that the star be rotationally stable, so that the (2-current velocities near the core 
are correspondingly limited, and just at the region where the choking effect of 
the p-currents is strongest. The result just found means that no uniformly 
rotating star can be rotationally stable and also suffer any mixing between core 
and envelope. 

In this paper attention is confined to Sweet’s driving velocity function, 
which assumes {2 constant; it is quite permissible, however, to imagine the inner 
portions of a star to be rotating uniformly at a high enough rate for some mixing 
to occur, while the outer regions rotate slowly enough to ensure stability at the 
surface. In the following work the parameter g, representing conditions at the 
surface of a uniformly rotating star, will still be used, but it must be borne in 





xF=F+ EA," =F- 


and so 
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mind that if the theory is to be physically relevant, the high value of the angular 
velocity must be confined to the inner regions. ‘The conditions for this to be 
possible will be discussed in Section 11; meanwhile, this first result shows very 
strikingly the difficulties that are encountered by the theory of rotational mixing 
between core and envelope. 

The next cases considered were those with 2’ = 1-3 and 1-5; the mixing zones 
then contain respectively 0-18 and 0-26 of the stellar mass. For each value of 2’, 
different values of (5A,A) were assumed and the equation of mixing (49) was 
integrated from both ends of the range simultaneously, using various values for 
xatz,and 2’. Ataconvenient point within the range two graphs of A, against A,’ 
were plotted, the variable parameters being the two values of x. _ 

For zs’ = 1-3, the two graphs were found for small enough (5.4,A) to intersect 
at two points (corresponding roughly to the positive and negative signs in (56)). 
The critical value of (5A,A) is given by the coincidence of the two solutions. 
One sequence of solutions is again rejected, as it gives x tending to zero instead 
of to unity as (5A,A) tends to zero. In each member of the physically relevant 
sequence, x decreases monotonically from z, to 2’. The remarkable result is 
found that the critical value of (5.4,A) gives from (51) a value of g smaller by a 
factor of about 0-96 than the limiting value of g found above for 2’ = z.,. 

In the case 2’ = 1-5 the details have changed somewhat. ‘The two graphs now 
have three intersections for small enough (5A,A), but two of them may be 
rejected, again because of their behaviour as q varies. The third is the physical 
solution ; it merges into one of the others at the critical value of (5.4,A), which 
this time gives for g a value slightly greater—by a factor of about 1-025—than 
that found for z’=<,. The function x now behaves differently; in the critical 
case it begins by decreasing from about 0-65 at z, and then increases to about 
2°5 at 2’. 

For larger values of 2’ results hold that are similar to those for z’=1°5. 
In each case there exists a critical value of (5.A,)\) which is considerably smaller 
than the value 0-05733 found for z’=z,, but when the larger mass of the mixing 
zone is taken into account, as in (51), the critical value of g is found to be only 
very slightly larger than 0-96. The minimum value of q increases monotonically 
with m’; for m' =m, (the non-dimensional stellar mass), Cowling estimated that 
the minimum value exceeds 0-96 by a factor of about 1-2. In all cases the 
variation in x follows the same pattern—a drop from x, to a minimum value near 
z=1-4, and then an increase outwards. 

Further integrations for z’ near 1-4 show that the critical value of g reaches 
again the value 0-96 for z’~1-46. Also, it is for a value of z’ in this neighbourhood 
that the function x ceases to decrease monotonically outwards from the core, 
but acquires a minimum in the neighbourhood of z= 1-46. 

The general behaviour of the solutions is due to the form of the driving 
function p(r). This is effectively the product of 1/(m—#) and Q?r3/GM(r), the 
first of which falls off rapidly near the core from infinity to a nearly constant value, 
while the second increases monotonically outwards. ‘Together these functions 
account for the minimum in f(r), which in turn produces the minimum in x. 
The extreme sensitivity of the mass of the mixing zone to values of g greater than 
0-96 is again due ultimately to the functional form of p(r), in that the »-currents 
induced by p(r) lead to an enhancement of the velocity on the other side of the 
minimum. 
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g. Interpretation of the results—Consider first a star with the parameter g 
smaller than 0-96, the critical value found above for the case m’=m,. From the 
form of the solution (56) it is seen that the value of (5.4,A) must begin to decrease 
a small distance away from the core, and, as pointed out in Section 7, such a 
decrease corresponds to a decay of current velocities which becomes progressively 
more pronounced with time. Thus for g less than the critical value 0-96, the 
{2-currents are checked as soon as they emerge from the core: the induced 
-currents are strong enough to choke the circulation, and a barrier is built up 
which prevents any further exchange of matter between core and envelope. 
The circulation in the envelope remains of Sweet’s type except near the core, 
where the barrier diverts the descending flow horizontally. Hence the only 
hydrogen available for energy liberation in this stage of the star’s evolution is 
that contained in the core. 

Now imagine g at the critical value 0-96. It was found above that with this 
value quasi-static circulation can occur in any zone between 2’ = z, and 2’ = 1-46. 
Thus matter will be able to stream out of the core and reach z’=1-46. The 
feed-back terms cut down the velocity, but are not strong enough in this range 
to choke the flow; the numerical results mean physically that the ,.-currents, 
which depend, roughly speaking, on differences between yz at pole and equator, 
do not for this value of g become too powerful anywhere, provided the motion 
is restricted to a sphere with z' <1-46. A previously existing barrier at ’ = 1-46 
would divert the flow horizontally, and so prevent the -currents from becoming 
too strong. But a Cowling model star contains no barriers initially, and the flow 
can be constrained only by barriers that it builds up for itself. Also, at z’ = 1-46 
the matter has a non-vanishing vertical component of velocity, and so in the 
absence of a barrier the initial matter from the core will try to flow further out. 
But it is known that for g=0-96, no quasi-static circulation can occur in a zone of 
mass greater than m (1-46), and so it is concluded that the condition A, = constant 
must break down somewhere in the range, leading to the decay of the motion, 
and the formation of a barrier region. 

It is difficult to determine rigorously where exactly this decay sets in, but it 
is certainly plausible that it should occur near the radius where the y-currents 
are already exerting their maximum choking effect, i.e. where the function wv, 
defined by (52), attains its minimum value. For g=0-96, it was found that 
x decreases monotonically outwards, and so is smallest at s=1-46. The flow 
for this value of g is therefore pictured as follows: the initial flow from the core 
is able to reach z= 1-46 without the u-currents becoming strong enough to choke 
the motion, but after this radius the motion decays, leading to the building of a 
barrier region. Further material coming up from the core is able to reach this 
barrier without being choked, but is diverted at the barrier, so that quasi-static 
circulation proceeds within the radius s=1-46. About one-quarter of the star’s 
mass is therefore available for energy liberation at this stage in its evolution. 

With g slightly larger, but less than about 1-15, quasi-static circulation is 
able to proceed within a larger mass, but not through the whole star. The mass 
of the largest mixing zone allowed increases monotonically with g. The argument 
of the last paragraph applies again; in the absence of an already existing barrier 
the initial currents will attempt to flow towards the edge of the star, but as 
quasi-static circulation through all the star is impossible, choking of the motion 
must set in somewhere in the range. For these values of g, however, the 
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function x behaves differently: it decreases monotonically from the core to 
21-46, and then increases monotonically. Again we expect the choking of the 
circulation to take place at the minimum of x, so this time the barrier is built up 
by the flow not at the edge of the maximum mixing zone allowed, for the given 
value of g, by the equation of mixing (49), but at the “ bottle-neck”’ near z= 1-46. 
Thus again only about one-quarter of the star is mixed, as the barrier at the 
“bottle-neck’’ separates the part of the circulation linked with the core from the 
part in the rest of the envelope, which has no influence on the evolution of the 
star. The solution of the equation of mixing describing the asymptotic state is 
that with m’ = m,/4, and with (5.4,) found from (51) with this value of m’ and the 
given value of g. 

In the absence of an equation to replace (44), it is difficult to estimate 
accurately the time taken by the currents to decay once quasi-static circulation 
is no longer possible. However, it is plausible that when choking of the motion 
starts, 1, should approach its asymptotic value at roughly the same rate as Ap, 
in which case it follows that the currents decay, and the barrier is built up, in a 
time comparable with the time for the matter to make one circulation, which is 
about 74, of the time required for all the hydrogen in the mixing zone to be 
used up. 

With g=1-15, the equation of mixing was found to have a solution with 
m' =m,, so that it appears that the whole star can be mixed by rotational currents. 
But a rotationally stable star cannot rotate in its outer regions with the high 
value of the angular velocity given by g=1, and so the driving velocity function 
p(r), which is proportional to 0?73/GM(r), must stop increasing well before the 
surface. ‘The function F of Section 8 will therefore decrease more rapidly than 
that calculated using Sweet’s function, and this will have some effect on the 
solutions of the mixing equation. However, it is found that with g=1, and 
Q? a constant, the function Q2r3/GM(r) remains less than *, up to 3, which 
contains 80 per cent of the mass, so that the bulk of the star can rotate with the 
high angular velocity corresponding to g=1 without (?73,GM(r) becoming 
dangerously large anywhere. ‘The integrations of the equation of mixing show 
that the extent of the mixing zone is controlled mainly by conditions near 
the ‘“‘bottle-neck”” at z=1-46; it is therefore unlikely that a decrease in Q 
after z=3 will do more than increase very slightly the value of g required for 
complete mixing. 

The present work has been based on Sweet’s first-order (2-current theory. 
In Opik’s treatment of the rotational problem (16), one term of the second order 
in Q? has been retained, as it involves the factor 1/p. Opik’s expression for 7,° is 


enn s1(e)(2—)]- 


The singularities given by this formula at the edge of the star are of no 
significance, as they can be eliminated by a technique similar to that developed 
in the Appendix. What is important is that for high enough 22? the vertical 
component is reduced to zero and then reversed before the edge of the star is 
reached. With the parameter g=1, v,® from (57) vanishes when M(r)~0-9M. 
With this formula used in computing F of Section 8, and the new function F 
used in the solution of (49), it is probable that the maximum allowed mixing zone 
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will contain less still of the star’s mass. This result will have considerable effect 
on the position of the star in the Hertzsprung—Russell diagram in the course of 
its evolution ; but it is desirable that the formula (57) be verified by an extension 
of Sweet’s analysis to include second-order terms. 

To sum up, there exist two critical values of the parameter g. In order 
that the material may be able to escape from the core at all, g must exceed the very 
high value 0-96. With g at this value, the mixing zone has a mass of about //4, 
and it retains this value as q is increased until the second critical value, which is 
only very slightly greater, is reached. According to the first-order theory, the 
whole star will be mixed once g=1-15; while if Opik’s second-order term is 
accepted, the maximum fraction of the star that will be mixed by any value of 
the angular velocity is about 3. 

Quite apart from the physical assumptions explicitly made, such as the 
form of the angular velocity field, the present results do depend on the validity 
of certain mathematical approximations. First, all but the P, terms have been 
neglected in the treatment. By considering a quasi-static solution of (26) for 
s’=s,, it is easy to show that the higher Legendre coefficients will choke the 
circulation even more near the core, so that the minimum value of g found above 
cannot be reduced by such refinements. A more serious criticism can be made 
of the step from the results of Section 8 to those of Section 9g. In Section 8 it is 
shown that for a given value of g there exists a series of mixing zones with masses 
up to a limit. ‘Thus, for example, the value g=0-96 was shown to allow quasi- 
static flow through a zone of mass up to m, 4, and the conclusion was drawn that 
a star containing just so much angular momentum will develop a mixing zone 
of mass m, 4, the best possible case. But the initial j.-distribution set up by the 
first currents from the core will be highly asymmetric, and it is possible that the 
p-currents in the early stages of the flow will exert a stronger choking influence 
on the Q-currents than the symmetrical ,-distributions discussed in Section 8, 
so that the barrier is built at a radius nearer the core than ’=1-46. This doubt 
can be settled only by tackling the formidable problem of following the motion 
from the first flow of matter from the core. However, with the assumed angular 
velocity field, the results are so sensitive to slight changes of q above the critical 
value 0-96, that errors arising from this cause are of little consequence. ‘The main 
results hold: a high angular velocity is needed in the inner regions for matter 
to be able to escape from the core, and a star with a spin slightly greater will 
possess a mixing zone containing a substantial fraction of the star’s mass. 

10. Flow near the edge of the mixing zone.—The condition A,=0 at the edge 
of the mixing zone implies a finite vertical velocity field there, and so an infinite 
horizontal component. According to the picture so far, the matter reaches the 
barrier with a finite vertical velocity, moves with infinite horizontal velocity to 
the other point where its streamline intersects the barrier, and then descends 
again. ‘The way the material in fact avoids this unphysical behaviour near the 
barrier has already been pointed out in Section 2: the viscosity of the stellar 
material prevents the matter from being accelerated horizontally beyond a certain 
high velocity, and the flow settles down into a state with the vertical component 
reduced to zero in a distance of order one kilometre, and the horizontal component 
reaching a value about 10° times the normal speed of the currents. Thus 4, 
does not become exactly zero, but certainly small enough for the condition 
A,=0 to be adequate for a discussion of the large-scale properties of the circulation. 
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The A,-distribution in the boundary layer may be determined as follows. 
As in (52), we write as 
xF = F+CA,+DA,' +EA,’. (58) 
Just outside the boundary layer x must take on the value at the edge of the mixing 
zone, written as x), found from the integration of (49); but within the layer the 
changed conditions will cause it to vary rapidly. ‘The vertical velocity due to 
both Q- and ,-currents will have multiplying it a further factor g, this being due 
to the action of viscosity in reducing the velocity to zero. ‘Thus the theory of 
the Appendix must be applied to the function (xg)F'; as F is slowly varying, the 
product (xg) is given by w(t —e*"-4), (59) 


where || is of order 10°. ‘The equation of quasi-static mixing then becomes 





AgxB F(1 — e*°’-*) =(5A,A)(m' — m), (60) 
so that A, is given by 
_ (54d) ys? _ (8-3) (8 = 83) 
A,= BxyF . (1 — eh?’ »)) =F (1 — eh *)) (01) 


say. ‘The feed-back terms are then determined; the first differs from zero by 
a term of order 7/||, while the other terms make large positive contributions 
to x of order 7 and y|k| respectively. ‘The functions x and g, giving respectively 
the influence of the feed-back terms and the viscosity on the flow, can then be 
found from (58) and (59). 

In a similar way, if the mixing zone should extend to the edge of the star, 
the positive infinite vertical component of velocity given by the solution of (49) 
must be reduced to zero by the effect of the stellar viscosity on the thermal 
equation, so that the horizontal component is finite. A modification of the theory 
in the Appendix will be required, and the 4,-distribution may then be found as 
above. 

It is to be noted that strictly speaking, quasi-static flow in the boundary layer 
is an asymptotic state, requiring an infinite time to be reached; for with the 
velocity behaving like (z’—) near z’, the time function ¢(z,@) becomes infinite 
at 2’ like —In(z’—:s). This result is in fact a simple consequence of demanding 
that the horizontal component at the barrier should be finite: for this to be the 
case, the vertical component must decrease at least as rapidly as (s’ — 3), so that ¢ 
must become infinite at the barrier. A similar result holds for flow near the edge 
of the star. 

It is in this connection that the dubious nature of equations (28) and (50) 
appears. With the quasi-static assumptions, A,’ is given by (50). If the vertical 
component is made to vanish at 2’, then A, and A,’ must remain finite, leading to 
a positive infinity for Ay’, whereas the positive infinity in the time function ¢ 
should produce a negative infinity in 4)’. Thus (50), and hence also the equation 
(28) from which it is derived, are quite unreliable near the edge of the mixing 
zone. Equation (28) has therefore not been used above except in the form (31), 
where the doubtful terms are omitted. 

11. The angular velocity field —It has been assumed above that the mixing 
zone continues to rotate uniformly throughout the star’s lifetime, so that the 
field of the driving velocity, and hence the function F, is also constant in time. 
For this to be the case, some azimuthal force must act on the fluid elements so as 
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to nullify the changes in angular velocity brought about by the transport of 
angular momentum by the circulation. Stellar viscosity is much too small to 
exert much influence, but rough uniformity of angular velocity is preserved if 
the mixing zone possesses a meridional magnetic field with lines of force threading 
the whole zone, including the convective core. ‘The high conductivity of stellar 
material causes the lines of force to be effectively frozen into the matter (17), 
so that a meridional magnetic field will have its lines of force distorted and 
stretched by non-uniform rotation. ‘The resulting azimuthal component of 
the ponderomotive force reacts back on the star, and sooner or later the non- 
uniformity of angular velocity along a line of force is reversed. ‘Thus angular 
momentum is propagated to and fro along the lines of force; if the lines pass 
through the convective core, which necessarily rotates as a rigid body, then the 
angular velocity of the whole mixing zone will oscillate about a constant mean. 
The period of the oscillation is roughly the period of a magneto-hydrodynamic 
wave along the lines of force, i.e. R/H(«/47p)'*, where H is the magnitude of the 
magnetic field, . the permeability, and p the mean density. This is small 
compared with the time of circulation of the currents if H is larger than 10-4 
gauss. ‘Thus quite a small magnetic field is both adequate and necessary to stiffen 
the angular velocity field of the mixing zone, so that the theory of this paper is 
relevant. 

However, it must be noted that a general magnetic field of this magnitude 
threading the whole of the star will destroy any chance of rotational mixing, not 
because of the very slight resistance it offers to flow across the lines of force (1), 
but simply because of the action of the field on the rotation of the star. Consider 
a star which begins by rotating rapidly enough in its inner regions for mixing to 
occur, but with the spin decreasing outwards so that rotational stability is 
preserved near the surface. Any meridional magnetic field present will redis- 
tribute the angular velocity by transporting angular momentum outwards. The 
effect will be to bring about rotational instability at the surface, so that the star 
must lose mass and angular momentum until the centrifugal force at the surface 
ceases to be of the same order as gravity. ‘Thus the star must reach a state with 
the magnetic field preserving a mean uniformity of rotation throughout, with the 
value of the angular velocity well below the limit found earlier for mixing to 
occur. A magnetic field greater than 10-* gauss will suffice to produce this effect 
within the time-scale for consumption of the hydrogen within the mixing zone. 

If any magnetic fields present are too small to influence the angular velocity, 
and if the inner regions of the star rotate rapidly enough, then circulation will 
take place, but the currents themselves will tend to set up a state of uniformity 
of angular momentum within the mixing zone. Also, equalization of angular 
momentum tends to be brought about by the braking effect of the radiative flux, 
as discovered by Jeans (18). The viscosity is negligible except near the axis of 
rotation, where strict conservation of angular momentum by the circulation 
would lead to infinities in the azimuthal velocity and in its gradient normal to 
the axis. ‘Turbulent viscosity will ensure that the core rotates as a rigid body. 

It is known (16) that rotational currents flow vertically upwards when the 
change of gravity due to rotation is positive; hence in a star in which (27? sin? @ 
is constant throughout most of the mixing zone, currents will flow out of the core 
at the equator, and in at the poles. ‘This ensures that matter leaving the core 
has a finite quantity of angular momentum to distribute through the mixing zone. 

50 














736 L. Mestel Vol. 113 


The details of the (2-currents for this case may be obtained by Sweet’s general 
method, and then the feed-back currents determined by methods similar to 
those of this paper. The treatment is likely to be much more complicated because 
of the necessity of introducing more than one Legendre component to describe 
the field adequately. A more tractable problem would be the circulation in a 
star obeying the law Qr?=constant; the Q-currents in this case again have a 
P,-distribution, so that the theory of this paper is immediately applicable once 
the function F of (47) has been recomputed. Such a law of angular velocity 
would give some idea of the effect of the uniform distribution of angular 
momentum, (27? sin? @=constant, without introducing too great complication. 

Some qualitative results may be predicted. First, the driving velocity will 
decrease outwards, instead of having a minimum, because of the monotonic 
decrease of the ratio of centrifugal to gravitational force. This will probably 
lead to a stronger choking effect by the p-currents, so that the results will not 
show the striking sensitivity of the mass of the mixing zone to small changes in 
angular momentum which was found for the case with (2=constant. The 
extent of the mixing zone will probably be limited by the amount of centrifugal 
force the inner region can safely bear. There will again exist a “ bottle-neck”’, 
but it must be much further out than that found above, for the total vertical 
component of velocity will not begin to increase until the large positive feed-back 
terms DA,’ and EA,” come into play at low densities far out into the envelope of 
the star. 

12. Evolutionary consequences.—The applications of the work of this paper 
to stellar structure and evolution inevitably depend on what assumptions are 
made initially as to the properties of massive stars at the time of their birth. 
Thus it follows, for example, that a star rotating with roughly the same angular 
velocity throughout, or possessing a general magnetic field threading the whole 
volume, must begin its evolution by transmuting just the hydrogen in the 
convective core. 

The results show that a star with its inner regions constrained by a magnetic 
field to rotate uniformly, and possessing enough internal angular momentum, 
will suffer some mixing; the zone will contain either one-quarter of the stellar 
mass, or, if the angular momentum be very slightly greater, the bulk of the star 
will be mixed. However, it seems unlikely that any general magnetic field 
possessed by the star will be restricted to just the inner regions, while if it threads 
the whole star, no mixing is possible. We shall therefore assume that if a star is 
able to mix at all, it is effectively free from magnetic fields. ‘This implies, as 
discussed in Section 11, that the Q-currents due to the new angular velocity field 
within the mixing zone will have different properties, and, in particular, the mass 
of the mixing zone is likely to be far less sensitive to increases of angular 
momentum. Thus it is probable that we cannot expect more than one-half of 
the star to be mixed, even under the most favourable conditions ; although further 
numerical results are desirable, for the moment we shall assume this to be the 
case. 

The general conclusion to be drawn is that rotational mixing is far more 
unlikely than has been supposed in much recent work on stellar evolution, and 
in the absence of any alternative mechanism of mixing, it is seen that the present 
results will tend to eliminate some stellar models as being inaccessible in the normal 
course of evolution. Thus it is known from the mass-radius-composition and the 
mass-luminosity-composition relations for a homogeneous star that as pz increases, 
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such a star describes in the Hertzsprung—Russell diagram an ascending path which 
diverges considerably to the left from the main sequence. Doubt has been felt 
recently by observers as to whether any such stars exist, and it is satisfactory that 
there should be strong theoretical objections to a star’s remaining homogeneous. 

Consider now the evolution of a star with a mixing zone extending into the 
envelope. Then within the zone » gradually increases, while the rest of the star 
is unaffected by the consumption of hydrogen in the core, so that the star will 
approximate very closely to the Hoyle-Lyttleton giant model (3, 4, 5). From 
the Bondis’ detailed integrations and graphs, it is possible to follow the track of 
the star in the Hertzsprung—Russell diagram, once the mass of the mixing zone 
is given. As y in the mixing zone increases, at first the luminosity increases 
strongly, while the radius increases only moderately ; but when the mixing zone 
contains already a high proportion of helium, the radius increases spectacularly 
for small increases in j., especially if the mass of the mixing zone is about half the 
total mass. Thus in the Hertzsprung—Russell diagram the star will at first ascend 
along a path diverging slightly to the right from the main sequence, but later will 
swing sharply to the right and move along an almost horizontal path, giving a 
large drop in surface temperature, accompanied by a slight rise in luminosity, 
as demanded by the integrations. Such a path is strongly suggestive of the 
main ‘Type I giant sequence (9). 

The way suggested by Hoyle and Lyttleton for the origin of red giants built 
on their model was the accretion of a hydrogen envelope by a completely mixed 
star which has had sufficient time to burn a large fraction of its hydrogen, and 
so has become a homogeneous star of high ». The present work shows at once 
that accretion is not necessary, as a star will almost certainly evolve towards a 
non-homogeneous state of its own accord. The Bondis’ graphs show that the 
most extended giant models are obtained if the masses of the zones of high and 
low p are roughly equal; if, therefore, not more than half of the star can be mixed, 
accretion of extra hydrogen will diminish the extension of the envelope. Thus 
accretion will certainly tend to blur the evolutionary tracks in the Hertzsprung- 
Russell diagram, but it will not make any qualitative changes in the structure of 
the star. 

It is known that some Hoyle—Lyttleton giant models, obeying a Kramers-type 
opacity law, are unable to remain in the static state—it is impossible to satisfy 
the conditions of continuity of energy flux everywhere, so that the star begins to 
gain energy (4). It is difficult to determine rigorously the future behaviour of 
such a star, but it is possible that it may get rid of surplus energy by performing 
non-linear oscillations (19), and so provide an explanation of some of the observed 
long-period red variables. If nearly one-half of the star is mixed, and if the 
opacity law obeyed is of the right sort, such a non-static star must necessarily 
arise in the normal course of evolution of a main-sequence star. 

We turn now to stars which suffer no mixing in their Cowling model stage, 
and so simply burn out their cores. As y in the core increases, the luminosity 
increases slightly, and the radius increases gradually up to ten times its original 
value (20). When the core is completely burnt out, the convection currents die 
out, and the star contracts until the temperature in the shell surrounding the 
core reaches about 15 million degrees, when the Bethe cycle can start up there. 
A small convective zone will surround the shell, so that the barrier built up by 
the -currents in the Cowling model stage, which prevented mixing between core 


and envelope, will be swept away. 
% 


50 
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It is well known that unless the mass of the core is a small fraction of the 
stellar mass, this “‘shell-source’’ stellar model cannot exist in the static 
state (6)—again, the star cannot transmit to the surface the energy liberated 
within it, and so must begin to gain energy. A plausible path for the first stages 
of the star’s evolution has been put forward by Sandage and Schwarzschild (21), 
with the core slowly contracting and the envelope slowly expanding, giving in 
the Hertzsprung—Russell diagram a nearly horizontal curve in the direction of 
decreasing surface temperature. This line of evolution must ultimately be 
modified, however, by the onset of electron degeneracy at the centre, together 
with the possible liberation of energy by the synthesis of heavy elements from 
helium, which starts at about 2 x 10° deg. K (22). The authors estimate that 
if helium burning starts, the rapid envelope expansion will cease, and the star’s 
luminosity will increase sufficiently to account for the bright stars in Type II 
systems. However, they are not certain that the central temperature reached in 
the first stage of contraction will be high enough for helium burning to start, 
and suggest as an alternative that some rotational mixing may start at some point 
on the star’s track. It is at this point that the present work may be relevant. 
It was noted in Section 7 that whereas the simple Q-current theory predicts 
the same efficiency whatever the star’s dimensions, the theory with ,«-currents 
results in greater mixing efficiency after contraction. Thus a star which in its 
Cowling model state suffers no mixing may, after contraction to a shell-source 
state, possess a mixing zone extending beyond the convective zone surrounding 
the burnt-out core. Such a zone would probably slow up the envelope expansion, 
because of the gradual decrease in the discontinuity between core and mixing 
zone, and at the same time would increase the luminosity. Further integrations 
of composite shell-source models are needed to compare the theory with 
observation. 

Once the mixing zone has consumed all its hydrogen, the mass of the burnt-out 
core effectively increases, so that the star will contract further. Electron 
degeneracy will probably play an increasingly important part in the central 
regions of the core, and we should expect both radius and luminosity to decrease 
quite sharply. It remains to be seen whether the theory can account in this 
way for the horizontal portion of the Type II giant sequence. It is in this part 
of the Hertzsprung—Russell diagram that the cluster-type Cepheids are found; 
it may be possible to find a criterion determining when a non-static star, which 
is continually gaining energy, prefers to start pulsating so as to get rid of the 
excess energy (19), and in this way account for the Type II short-period variables. 

It may be remarked that Type II systems seem to be completely free of 
interstellar gas and dust. Any massive stars that condensed when these old 
systems were formed should have burnt themselves out, and there is no matter 
available to form new massive stars, either directly, or through accretion by 
moderate masses. Hence it is not possible that the very luminous red stars 
should be of large mass; the high luminosity must be due to composition rather 
than mass, and a fair-sized mixing zone of high » surrounding the core is one 
possible explanation. 

If, as suggested above, a Type I giant is built on the Hoyle—Lyttleton model, 
implying that a fair degree of mixing takes place in the first stages of evolution, 
while a Type II star begins by consuming only its convective core, then there 
should exist a systematic difference in angular momentum between the two 
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populations. It is perhaps suggestive that Type I systems are confined to the 
galactic plane, and so have condensed from matter with comparatively high 
vorticity, while the typical Type II systems, the globular clusters, form a roughly 
spherical system around the galactic centre. Also, the theoretical differences 
that must exist between the evolutionary tracks of stars which do and do not mix, 
respectively, offer attractive possibilities for explaining the two population types. 
If, however, the conditions found above for rotational mixing to be possible are 
too stringent for stars of either type, then we are faced with the problem of 
explaining both types as being due to the evolution of stars which initially burn 
out just their convective cores. The differences between the two types would 
then presumably be due to the differences in composition reported by Greenstein 
and Schwarzschild (23). It seems worth while working out the consequences 
of all the different assumptions—mixing or no mixing in early or late stages of 
evolution—to see if the observations definitely rule out all but one set of 
possibilities. In particular, if evolutionary theory is found to demand mixing 
at some stage or other, then the conditions found in this paper will need to be 
satisfied; in this way we may gain some information about the distribution of 
angular momentum in initial condensations, or the presence or absence of stellar 
magnetic fields. ‘This information could conceivably act as a censor on some 
theories of stellar origins. . 

So far we have discussed evolutionary problems in terms of the positions of 
stars in the Hertzsprung—Russell diagram. The mass-luminosity diagram is in 
general not so useful, because of the difficulties in estimating stellar masses 
accurately, but is of equal theoretical importance. Recently some Soviet 
astronomers, in particular Fessenkov (10), have drawn attention to discrepancies 
between the close fit of most stars observed to the mass-luminosity curve, and 
theoretical predictions from the mass-luminosity relation for homogeneous stars. 
It is well known that a high value of « throughout a star leads to a high internal 
temperature, and to a consequently high flux from core to surface. Thus in a 
star obeying the simple Kramers opacity law, the relation is Lxp*5M*> (2), 
and so a homogeneous star containing a high proportion of helium will be of far 
higher luminosity than one of the same mass consisting of hydrogen. To avoid 
this difficulty, Fessenkov proposes that the stars emit corpuscular radiation at a 
rate proportional to their luminosity, so as to preserve the closeness of fit with 
the observed mass-luminosity curve. Other writers, notably Struve (9), also 
support the idea of corpuscular radiation, so that the initial stellar condensations, 
which they believe to be of high mass, may be able to lose enough mass so as to 
move down the main sequence into the dwarf regions. However, whatever the 
arguments in favour of or against this last line of evolution, it is seen at once that 
if the natural tendency of a star is to develop into a non-homogeneous state, and 
burn out initially only a fraction of its mass, then Fessenkov’s argument is no 
longer valid. The mass-luminosity relation depends on assuming that energy 
is flowing from the core to the surface, and that the perfect-gas equation of state 
holds throughout the star. But if a star begins by burning out less than half 
its mass, as seems likely, there will not be any serious divergence from the 
mass-luminosity curve while the star still possesses a convective core ; and when 
the mixing zone is completely burnt out, and the star contracts to a shell-source 
state, it is likely that degeneracy will set in over a large part of the core, and so 
there is no reason to expect the ordinary mass-luminosity relation to hold. What 
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relation replaces it cannot be decided until more theoretical work has been done 
but there is no reason to demand corpuscular radiation (by a mechanism as yet 
undetermined) so as to avoid theoretical difficulties. 

13. Conclusions.—(1) Very slight departures of the distribution of matter from 
spherical symmetry lead to strong thermal currents within stars (‘‘.-currents’’). 

(2) No continuous mixing between core and envelope can take place in a 
uniformly rotating Cowling model star; the y-currents induced by the rotational 
currents choke the flow near the core, unless the angular velocity of the star is 
so great that it is past the limit for rotational stability at the surface. 

(3) A star rotating rapidly enough near the core will build up for itself a 
mixing zone of definite mass. If the inner regions of the star are constrained by 
a magnetic field to rotate uniformly, then the mass of the mixing zone is very 
sensitive to increases of the angular velocity above a certain limit. In a star in 
which the angular velocity decreases outwards, however, the mass of the mixing 
zone is likely to increase slowly with increasing angular momentum. 

(4) The mass of the mixing zone is determined by a parameter which varies 
as the star’s dimensions change. ‘Thus the efficiency of mixing increases as the 
star contracts ; in particular, a star which is unable to mix in its initial state with 
a convective core may possess a mixing zone after contraction into a shell-source 
State. 

(5) The effects on the Hertzsprung—Russell diagram and on the mass- 
luminosity curve depend to a great extent on the assumptions made about the 
initial angular velocity distribution of the star, and the existence and strength of 
stellar magnetic fields. ‘The general tendency of the work is to make rotational 
mixing much less probable than has beenthought. The solution of the shell-source 
problem is urgently needed, to see whether both Type I and Type II giant 
sequences can be explained without use of the Hoyle—Lyttleton giant model and 
rotational mixing. The arguments of Fessenkov and others for corpuscular 
radiation, based on discrepancies between the observed mass-luminosity relation 
and the theoretical predictions for homogeneous stars, become much less 
compelling when it is recognized that it is very unlikely that stars remain 
homogeneous. 
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APPENDIX 


It was pointed out in Section 2 that at the edge of a star, near surfaces of 
discontinuity in composition, and at the interface between convective and 
radiative regions, the theory of meridional circulation so far developed gives 
infinities for the circulation velocities. It is concluded that the assumptions of 
negligible viscous force and inertial force, underlying the theory, must be invalid 
near certain singular surfaces. If the circulation velocity is expanded in powers 
of the perturbing parameter (27, then in a first-order treatment only the viscous 
force is retained, as the steady-state inertial terms are non-linear in the circulation 
velocity. ‘The theory of flow near the singular points will now be given for the 
different cases. 
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1. Motion near discontinuities and near the surface.—To the first order in {2?, 
the steady-state equations of motion are 


0d oP ) 2 a 
Ps, + por Q? sin 19 F +n( ve, - —— U - 35) =°. (1) 


0¢ a oP zy 2 dv, 

Po 59 + por?{2? sin 6 cos 6 — Btn? ve wang + BH) (2) 
where ¢ is the gravitational potential and 7 the coefficient of viscosity. These 
equations are cross-differentiated and subtracted one from the other, leading to 

Op ddy © hs 2v, 2coté 2 Oe 
060 dr +135 (¥e,- =) aaa ai %— 3290 


_ Op dpg a . 0 — Ue 2 dv, 
= 35 ar + a iQ sin @ cos 6 +7 = r Vve— zante t 2 30 (3) 


We consider first the motion near a surface of discontinuity of » within the 
star. The type of solution sought is one in which the effect of viscosity is felt 
only within a thin layer near the discontinuity. The vertical velocity at the 
barrier must be zero; provisionally, therefore, we adopt for v, an expression of 


the form 

v,=pP{1—e%"”) =a(r) Po, (4) 
where k is a constant, 7 the radius of the sphere containing the barrier and pP, is 
Sweet’s vertical velocity component, proportional to 2°. The constant & is 
positive below the barrier, negative above. It is shown in the analysis that | k| 
is so large that |r —7| need be only a small fraction of the radius for (4) to reduce 
to Sweet’s expressions; hence variations in p, ¢9, py, etc. over distances 
comparable with |1/k| are negligible. 

By the equation of continuity, 


ee sin 6 cos 6. 


oa I 
eel 2por dr 


Then the terms involving the viscosity in (3) are respectively 
, d 20 Id (por?x)’ \’ 
—<-— 2 ie eee 
3 sin 8 cos 6 = - ~ (bur a)’ — > 4 pe ( aay \ (6) 


_sinsoun {4 ih (72a) — 5 + aie’ (7) 


and 


The dominant term in these two expressions is the one with the highest 
derivative in r, as this has the highest number of & factors. With this term alone 
retained, (3) reduces to 


a a 
aA — = 0 dr - = r2Q)2 sin @ cos 6+ 4nr*k'p er 7) sin 0 cos 6. (8) 


The viscous term modifies 0¢/0@ and dp/d@. Let y=0¢/00 have added to it 
a new term x». Then dp/0@ has a new term given by Poisson’s equation 


op 
Vv? X2= —47 7G ik (9) 
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Hence yz satisfies 


Ls (V2x2) + x2P9 = — dnr*htp e"-” sin 8 cos 0. (10) 


x2 Must become negligible as soon as |r —r|>| 1/k|; a solution of (10) satisfying 
this condition is 


Ae" sin 0 cos 0. (11) 
Again dropping all but the dominant terms, A is seen to be 


A=- —< nr*k*p, (12) 
$o 


We now require the new terms that must be added to (div H). From the 
expression given by Sweet (x) for (div H)/00, it follows that the largest 
contribution comes from the term 


L 1 # (eT, 
qari T,' aF\ 90 )° (13) 


where 7, is the modification in 7 due to the viscous terms. From (1) and (9) 
the dominant term in 07,/00 is 
a PB 
k(r- ‘ S 
— a ** ™ sin 6 cos @. (14) 
0 
This gives for ¢ (div H) ” 
F I k®npr*T, 
477 poly 2b 
The viscous heating terms in the energy equation are of the second order in 
the circulation velocities, and so are dropped. Hence we have 


Pol — $0) _ 3) 


(n+1)° 
so that (15) must be put equal to 


(n— 3) 
— Pol — te) an) '® 


for the change in the vertical component, —p e*’-”, is due to the modification 
in (—div H) which in turn is caused ultimately by the viscous terms in the 
hydrodynamic equations. 

Finally, k is given by (15) and (17) : 

a Sipo (do )? ° Ty 3(n— 3) (18) 
v° LT, ~° (n+r)’ 

where v is the kinematic viscosity 7/p). In most stars the radiative kinematic 
viscosity 2a7*/15«cp” (« here being the opacity) is somewhat greater than the 
molecular viscosity (18); both are of order unity (in c.g.s. units), except near the 
surface of the star. Hence (18) gives |k]|~10~, prov ided y is not too near the 
surface, so that the thickness of the boundary layer is only about a kilometre, 
fully justifying the approximations made in determining (18). 

In fact, the velocity function near the barrier will not be as simple as (4). 
As presented, the theory satisfies the condition of zero flow across the barrier, 
but in addition there are conditions of continuity of x2, x2) and 7; across the 
distorted discontinuity, while the jump in p, is fixed by the jump in w. These 
conditions themselves imply a restriction on the viscous terms in (2). In order 
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— div H, (16) 


ep . 3 sin 6 cos 8, (17) 
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to satisfy all these conditions it is necessary to assume a generalized form of (4), 
involving several exponential terms, and to take explicit account of conditions 
on the other side of the barrier. This greatly complicates the theory; we are 
concerned merely to illustrate what happens near the discontinuity, and to obtain 
an order-of-magnitude estimate for the thickness of the boundary layer, and so 
the accurate analysis is not attempted. 

Before the steady state is reached, material is accelerated up to the limiting 
velocity, so that both viscosity and first-order inertial terms should be included 
in equations (1) and (2). However, it is easy to see that the time for the system 
to reach the steady state is very small. The energy per unit mass that would be 
convected away per second by vertical motion in the absence of a barrier is of 
order (LM). (Q?r3,.GM). The presence of a barrier causes horizontal motion, 
and this energy is used to accelerate the material until the velocity gradient is 
large enough for viscosity to be important and to stop the acceleration. Hence 
the time to reach this state is given roughly by 


” L (2273 
v*~t . (ai ° Gin) ; (19) 
v is here the horizontal component, which from (5) is ~|k|rp. Taking 
r=5x 10! cm, and p to be ~10~1°(Q)? cm/sec, t is found to be a small fraction 
of a second. 

It is pointed out in the text that a finite vertical velocity at the surface implies 
an infinite horizontal component there, which again is physically impossible. 
The difficulty is resolved in a similar way: material is initially accelerated 
horizontally until the viscous forces come into play. The radiative viscosity, 
which far outweighs the molecular at the surface, is proportional to the 
temperature; also, as density and temperature vanish at the surface, a slightly 
more complex theory than that above is required. For a star obeying Kramers’ 
opacity law, it is found that the velocity vanishes as 

v,=pP,{1 —exp(—h/h)*2?}, (20) 
where h is the distance from the surface. Thus both vertical and horizontal 
components vanish at the surface. The thickness of the boundary layer is of 
order h, which for the Sun is about 2°5 x 108 cm and is not much larger for more 
massive stars. This is again small enough to justify the approximations made 
in arriving at (20). 

2. Motion near a convective zone.—In fact, the velocity distribution (20) will 
not hold near the surface of stars of solar order, because of the existence of the 
Uns6éld convective zone, extending in the Sun from a depth of about 10° cm 
up to the photosphere; within a convective zone meridional circulation is 
irrelevant. However, as the polytropic index approaches the adiabatic value #, 
the vertical and horizontal components of velocity are by (16) and (5) 
proportional to 1/(7—#) and 1/(m—¥)* respectively. Similar infinities arise at 
the surface of the convective core, and must be eliminated by a technique similar 
to that used above. 

The velocity function assumed is 


v,=l.s.P,=p(n—3).s. P,, (21) 


where p is Sweet’s function. As a convective zone is approached / remains finite 
and of the same order as p at points far from the convective zone; s is a rapidly 
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varying function of r which becomes large but finite at the convective zone, and 
behaves like 1/(n—¥%) far from the singularity. 

It is sufficiently accurate to assume that (—) varies linearly with (r—,r,), 
where r, is the radius of the convective zone. ‘Thus for (n—%) we write 
(dn/dr), .(r—r,), and assume for s the functional form 


I 
*= (dnjdr),.(r—r,) ° 
with the condition 
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Then as above, dp./00 is given by 


ry 1. sin@cosé 
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and 0T,/00 by 


aT __ Ty dbs 
00 -~—sopy 80” 
where only dominant terms have been retained. 
The energy equation (16) reduces to 


{ Faz eee hay a >. aa ert, 


t=1 n=1 n! t=1 ‘n=7(n—7)! 


where y is written for 

L =!) I 7 (n+t) 

4n° - Ty] * (~40') * ©p¢? * (dn dr), 

On equating coefficients, we obtain an infinite homogeneous set of linear 
equations which, together with the non-homogeneous condition (23), yield an 
infinite determinantal equation for yk’ =u: 


o= Sy 4° 


) 2-8) (8) fed 
(-8) 6-8) (48) (8). 


In a rigorous treatment it would be necessary to show the existence of only 
one physically relevant root of this determinantal equation as the number of 
z-terms taken increases to infinity. This will not be attempted. Inspection of 
the first few members of the sequence of determinants, however, shows that 
there is a root of order —10. Near the convective core y is roughly 10-* c.g.s. 
units, so that |k|~3 x 107° c.g.s. units. The thickness of the layer is then about 
$ km, and the vertical and horizontal components of velocity are 3 x 10~* (Q)? 


and g x (2)? cm/sec respectively. Near the Unséld convective zone | k| is about 
10-6 c.g.s. units, giving 10 km for the thickness of the layer, and 10-4 (Q)? and 











No. 6, 1953 Rotation and stellar evolution 745 


3 (2)? cm/sec for vertical and horizontal components respectively. It is 
remarkable that even in a slowly rotating body like the Sun the very slight 
thermal disequilibrium set up by the rotation can lead, near an unstable region, 
to velocities of an observable magnitude. 


Department of Mathematics, 
The University, 
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PHOTOMETRY OF THE GALACTIC CLUSTER NGC 6025 


A. R. Hogg 
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(Received 1953 November 3)* 


Summary 


The southern galactic cluster NGC 6025 has been examined photo- 
electrically in two colours. A catalogue of magnitudes and colours of 
84 stars in the region has been obtained. The colour-magnitude array is 
examined and compared with curves obtained for the Pleiades by Binnendijk 
and by Eggen. 

These comparisons give distance moduli of 9™-4 and 9™-6 respectively 
(neglecting absorption). Owing to a lack of independent identification of 
cluster members the results do not provide a positive test of the fine structure 
suggested by Eggen, but an examination of the stars in the central region 
yields nothing against the fine-structure interpretation. 





1. Introduction.—The galactic cluster NGC 6025 is included in the catalogues 
by Trumpler (1) and by Collinder (2), in both of which the position, taken from 
Franklin Adams charts, is given as «(1900) = 15" 55™-2 and 6(1g00) = — 60° 13’ fF. 
Collinder describes the cluster as “‘ a loose aggregation of stars of the Pleiades 
type’; Trumpler places it in class 113p, i.e. a detached cluster with but little 
central condensation, composed of moderately bright and faint stars numbering 
less than 50. ‘Trumpler estimates the modulus m— M=g9"-5, but considers this 
to be less satisfactory than most of his results as it is derived from the spectra 
of only five dwarf stars classified as B3 to B8 in the Henry Draper Catalogue. 
The giant branch is completely missing. Mrs C. A. Riecke (3) obtained a 
modulus of m— M=8"-4 using absolute magnitudes derived from an estimate of 
hydrogen line intensities. The estimated dimensions of the cluster, as sum- 
marized by Collinder (2), range from 10’ (Melotte) to 19’ (Raab). No 
detailed investigations of the spectra or radial velocities have been located nor 
does the cluster appear to have been the subject of photometric study. 

In the present work the cluster has been examined in two colours photo- 
electrically, with the object of obtaining a colour-luminosity array from which a 
distance might be derived and a classification of members according to the system 
of Eggen (4) obtained. 

2. Photometric technique.—The observations were made with an II-stage 
EMI photomultiplier tube attached to the Newtonian focus of the Reynolds 
30-inch reflector of the Commonwealth Observatory. The photomultiplier fed 
into an amplifier designed by Kron (5). Records were obtained with a 0-1 
milliamp recording meter. Blue and yellow filters were used and the star light 
was admitted through a focal-plane diaphragm subtending usually 50” of sky. 
The transmissions of the filters were measured by the National Standards 

* Received in original form 1953 May 19. 


+ Values of 154 54™-6 and —60° 10’ have been derived by plotting positions of CPD stars 
and reducing to the equinox of 1900. 





No. 6, 1953 Photometry of the galactic cluster NGC 6025 747 


Laboratory, Sydney, and, together with the makers’ data for the response of the 
cell, are shown in Table I. 
TABLE I 
Percentage transmission of filters and cell response 
A(A) F 3653 F 4937 Cell A(A)  F 3653 ~F 4937 Cell 
4000 ‘ . 54°5 5400 fohe) QI°5 58:1 
4200 . : 59°5 5600 oor) g1° 42°6 
4400 i ‘ 63°8 5800 oho) gI- 30°8 
4600 ° . 65:2 6000 oo gi° 16°5 
4800 : . 64:8 6200 oer) g1°4 79 
5000 ; : 62°7 6800 0-0 90°7 | are) 
5200 ‘ . 61°8 7000 o'o go°4 fohe) 


Generally the method used was to expose the cell to a star for I minute 
through the blue filter, close the photometer for } minute, open through the 
yellow for I minute, close for } minute and finally open for 1 minute through the 
blue filter. The performance of the photometer was controlled by a radio- 
active luminous source and sky corrections were applied. Absorption corrections 
for colours were applied with an average coefficient throughout, based on 
results obtained by Woolley and Gascoigne (6) in the course of a spectrophoto- 
metric programme. The adopted values were k (blue)=0™-34(sec z) and 
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Fic. 1.—Chart of NGC 6025 region. Stars are marked with the last three numbers of the 
CPD designation. The broken line encloses the “ central region’? as defined by “* star counts”, 
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k (yellow) =0™-24(sec 2), the difference o™-14(sec z) being used to correct the 
observed colours to a standard zenith distance such that sec z= 1-1, correspond- 
ing to that of the cluster at culmination. The use of such averaged coefficients 
may not be entirely satisfactory but as the measures were conducted during a 
spell of unusually cloudy weather it was decided to omit special measures of 
absorption on each, more especially as the corrections to be applied were usually 
less than 0™-03. 

3. Programme and results.—The observational programme aimed at determin- 
ing colours and magnitudes of all CPD stars between the limits «(1875) = 15" 51™ 
to 15" 54™ and 6(1875)= —59° 50’ to —60° 15’. Certain regions were examined 
more closely to obtain measures of fainter stars. The stars are noted on the 
chart shown in Fig. 1. Four stars in the cluster were standardized by connection 
to working scales termed 6 and y which had been derived in the course of unpub- 
lished photoelectric measures of the Harvard Fri region. The 6 scale was 
related to the Harvard Photographic Scale (HPg, with the correction given in 
Harvard Bulletin 718) by the expression 


pPg=6b+0-13—0-37Cpe=HPgto-itrs.d,, (1) 
using the results from 24 stars. In this expression the photoelectric colour 


Cpe is measured as the magnitude difference of the yellow and blue readings and 
is related to the difference b—¥y by the equation 


Cpe=b—y— 0-43. 
Photovisual magnitudes were obtained by the comparison of 25 stars in the 


Fi region which were available from the Woolley, Gottlieb and Przybylski 
measures (7) and gave the result 


pIPv=y+0-35—0-12Cpe=IPv + 0:05 s.d. (2) 

The colours were obtained as the difference between (1) and (2), which reduces to 
pCi=0-75Cpe + 0-207. 

The results are shown in Table II, which gives the following particulars :— 


Col. (1) CPD designation of star. (Stars not listed in CPD have been lettered 
and referred to an adjoining CPD star. ‘The lettered stars are not in 
order of right ascension.) 

(2) Magnitudes in the 5 system used for the F1 region. 

(3) Dispersions (standard deviations) of 6 in hundredths of a magnitude. 

(4) Colours, Cpe, measured directly in thousandths of a magnitude. 

(5) Standard deviations of Cpe. 


(6, 7) Number of observations, 7, on n nights. 
(8) Magnitudes pPg. 
(9) Colours reduced to International scale pCi. 
(10) Positivn in cluster, c= central, o= outer. 
(11,12) HD spectral type and number when available. 


4. Colour-magnitude relations.—The colour-magnitude relation is shown in 
Fig. 2 where the coordinates are pPg and pCi. ‘Two alternatives present them- 
selves in interpreting these results. The generally accepted method is to regard 
the cluster members as possessing an inherent scatter in the colour-magnitude 
relation, a scatter which is sometimes referred to as the “‘ cosmic dispersion ’’. 
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Colours and magnitudes in the NGC 6025 region 
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Mean sd. r n PFS X 1000 
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* This reading includes the light from two stars. 

+ This star is near 6329 the light from which may be responsible for some of the apparent 
scatter in the readings for 6332. 

t Small deflections and large scatter. 
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TABLE II (cont.) 

CPD b 1000 Cpe Obsns pCi ; 
No. Mean - s.d. Mean s.d. r pPg X 1000 Posn op 
6366 11°97 + 773 28 11°81 + 787 
6368 11°13 320 18 11°38 — 033 
6370 11°05 320 10 11°30 — 033 
6371 11°33 ’ 1200 11°02 > 1107 
= 
6539 11°78 : 243 
6539a 13°17 ogo 
6539b 13°31 ‘ 1550 
6539¢ 12°81 055 
6540 11°71 095 
6541 10°87 — 344 
654I1a 12°53 — 50 
6541b 14°52 oes +1340 
6541C 11°63 + 020 
6541d 11°94 ae — 408 
6549 11°66 + 050 
6553 10°87 — 380 
6556 11°24 — 220 
6557 8-81 — 583 
6558 10°81 — 375 
6562 9°64 — 493 
6562a 13°70 +1120 
6562b 12°58 160 
6562c 13°66 140 
6563 9°30 515 
6568 10°33 413 
6575 11°40 248 
6576 8-17 285 
6578 10°12 395 : 
6579 11°69 - 895 45 
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The points for the cluster members then fall within a band in the array. Such 
results are typified by the figures of Binnendijk (8) for the Pleiades. In recent 
years Eggen (4) has produced a considerable amount of observational data 
leading him to conclude that the observed bands in such a diagram possess a fine 
structure. Certain observers do not agree with this result, e.g. Johnson and 
Morgan (9). ‘The present observations do not lend themselves to provide a 
precise test of fine structure, mainly because of a lack of independent evidence of 
which stars are actually members of the cluster. However, as will be shown 
later, the results may be interpreted either in terms of a fine structure (FS) or of 
a band (B). In the sequel both methods are dealt with, without any implication 
that either is correct. 

Factors tending to cause a colour-luminosity array to depart from a close 
regular relationship between the two quantities are :— 


(a) Observational error dispersion of the results. 

(6) Distance between the back and front of the cluster. 
(c) Inclusion of close double stars in the array. 

(d) Existence of localized reddening within the cluster. 
(e) Inclusion of non-cluster stars within the array. 

(f) Existence of the so-called ‘‘ cosmic dispersion ’’. 


(a) Error dispersion of the observations.—The standard deviations given in 
Table II average about + 0™-025 in colour and +0™-03 in magnitude on the Cpe 
and 6 scales. ‘The equivalent figures on the pCi and pPg scales are + 0™-02 and 
+0™-03 respectively. ‘These of course refer to the internal consistency of the 
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observations only, and are presumably due mainly to instrumental behaviour, 
errors in scale reading, and atmospheric variations. As in the present instance 
the colour parameter is the more important one in determining departure from a 
curve, a criterion has been adopted that any result departing more than about 
o™-025 in Cpe from a suggested curve, except in a few cases where the s.d. differs 
greatly from its mean, would be unlikely to belong to that curve. The dispersions 
of the present results are rather greater than is usual in photoelectric work at 
Canberra (viz. 0™-o1 or better in magnitude), but in order to make the most 
use of available telescope time, much of the work has been carried out on nights 
that have not been of the first quality and this has increased the dispersions. 

pPg 








tT tT T qT T ' T T T 
0-4 -0-2 0-0 +0-2 +0-4 10-6 +0-8 +1-0 41-2 -1-4pc i 
Fic. 2.—Colour-magnitude diagram for NGC 6025. The full curve is taken from Binnendijk’s 
observations for the Pleiades displaced on magnitude scale by 4:6(3) units. The broken line shows 
Binnendijk’s curves adjusted by 0™-04 and 4™-41 in colour and magnitude respectively. 


(b) Distance between back and front of the cluster.—For spherically shaped 
clusters of relatively small angular diameter 0, the apparent magnitude difference 
between equally luminous stars situated at the nearest and farthest points of the 
cluster is approximately 26 (@ in radians) or 0-0006@ (@ in minutes of arc). 
In the present case 6 does not exceed 30’ and the difference is 0™-o17. This 
would be quite negligible in the type of diagram found. 

(c) Inclusion of close double stars.—The observed result for an unresolved 
binary system is displaced from the result for the brighter component in a colour- 
magnitude diagram by an amount dependent on the magnitude difference between 
the components. Thus for zero magnitude difference there is a displacement of 


SF 
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0-75 unit in apparent magnitude without change of colour (if the system is a 
member of the cluster). Holmberg (10) shows for the Hyades that a maximum 
colour displacement of 0™-05 occurs when the magnitude difference is 1™-5—2™-0, 

No catalogues have been found to give any specific information of any spectro- 
scopic or other close double stars in the cluster. It is possible that some of the 
stars later marked as non-members of the cluster might in actual fact be close 
binary systems belonging to it. 

(d) Existence of reddening in the cluster.—The results were examined to deter- 
mine whether any special areas within the cluster were occupied by stars that 
were systematically redder than might be expected from the general trend of the 
curve in Fig. 2. No clear evidence of any areas of differential reddening was 
obtained although it is possible that an area to the south-east of the cluster may be 
slightly affected in this way. 

(e) Inclusion of non-cluster stars in the array.—This may be the principal 
cause of the scatter observed in Fig. 2. Few observations are available which 
would enable any given star to be positively identified as a cluster member. 
Mr Wood, Government Astronomer of Sydney Observatory, has kindly made 
available some unpublished results of proper motions in the cluster region. He 
considers that the cluster stars have proper motions of something less than one 
second per century, relative to the stars used as comparisons, and that, at the 
present stage, it is not possible to identify members positively by their motion. 
Certain stars of larger proper motion may however be excluded as members of 
the cluster. In these circumstances the results have been treated statistically 
with the object of ascertaining whether the probable number of field stars in the 
central region of the cluster agrees with the number of observations in the same 
region departing from the system of curves determined by Eggen for the Pleiades. 
For this purpose star images on an astrographic plate of the region, taken by 
Wood, were counted. ‘The plate showed stars to about the same magnitude as 
those measured photometrically in the central region. The number of stars per 
square of 10’ side in regions definitely outside the cluster was 6-5 + 2:8 (s.d. per 
single square). It was taken that any square for which the residual was more than 
twice the s.d., i.e. any square containing more than 12 stars, was definitely within 
the cluster. The measuring réseau had been adjusted so that the bulk of the 
cluster was contained in one square which was found to contain 26-5 stars. An 
adjoining square contained 14-5 stars. It was obvious that both squares were 
associated with the cluster but a more detailed examination was needed to deter- 
mine what proportions of the square containing 14-5 stars was covered by cluster 
stars. Smaller regions of the square were examined and this led to the limits 
shown in Fig. I being adopted as representing the central region of the cluster. 
A check was obtained by comparing the density of stars over all the area (viz. 
130 square minutes) within these limits and density in the 26-5-star square 
which was clearly a cluster area. ‘The full area of 130 square minutes would be 
expected to contain 34°5 stars on the 26-5-star basis. ‘The actual number of 
stars counted on the plate was 32, showing the adopted area to be substantially 
correct. Further, the maximum number of field stars to be expected in the area, 
allowing a dispersion of 1-5 times the standard deviation as likely, would be 12-5, 
a figure which will be referred to later. 

Interpretation in terms of colour-luminosity band (B).—The colour-magnitude 
diagram for all the observed stars is shown in Fig. 2, where the results are plotted 
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in pPg and pCi coordinates. The full curve there fitted to the observations is 
taken from Binnendijk’s results for the Pleiades cluster, displaced on the magni- 
tude scale by 4-6(3) units. As the modulus for the Pleiades is 5-0, this is 
equivalent to assuming a modulus of 9-16 for NGC 6025. The colour coordinate 
has been left unchanged. Although the fit is not very satisfactory the diagram 
can be used to examine the ‘‘ cosmic dispersion” of the colour-luminosity 
relation. For this purpose stars whose colour falls within 0™-05 of the line in 
Fig. 2 have been regarded as cluster members. This gives a cluster membership 
of 46 stars, but if the permitted residual is taken as 0™-1 then seven more stars 
are to be included in the cluster. The s.d. of one of the 46 points from the line 
amounts to +0™-029 in colour. As the (internal) s.d. of an observation is 
+0™-o019 the “ cosmic dispersion” is calculated as +0™-o21. This may be 
compared with the figure of +0™-059 given by Holmberg for the case of the 
Hyades cluster, where the stars have been selected by proper-motion criteria. 
However, the present figure is dependent largely on the scatter permitted in the 
selected stars. Its main use is to indicate that the criterion adopted in the present 
case is not less strict in the mean than that employed in the known case of the 
Hyades cluster. The 46 stars to be regarded as cluster members by this method 
are listed in Table III. One of these stars, viz. 296 (i.e. CPD —60° 6296) is 
reported by Wood as having relatively large proper motion. 


TaBLe III 
NGC 6025 
Cluster members by band interpretation 

286c *319 *344 
*286e "sa% 345 
*204 323 *349 
296 *325 350 
*207 326 *355 
*298 *332 300 
303 *334 362 
305a *335 *368 
*310 *336 *370 
31 *338 539 
*316 *339 539¢ 
"907 343 541 
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* Stars marked thus are given as cluster members by both band and fine-structure 
interpretations. 


A more satisfactory fit for the brighter stars is obtained if both the colour and 
magnitude coordinates are adjusted, by 0™-04 and 4™-41 respectively. ‘This 
assumes a modulus of 5™-0+4™-41=9™-41 and allows for a discrepancy of 
o™-o4 in the two colour scales, which might have arisen observationally or have 
been brought about by the Pleiades being slightly reddened relative to NGC 6025. 
This adjustment is shown by the dotted line in Fig. 2. Neither curve gives a good 
representation of the results over the whole range of magnitude. 

Interpretation by fine-structure curves (F'S).--The colour-magnitude array 
for the central region as defined by star counts, plotted on the pPg and pCi scales, 
is shown in Fig. 3. ‘The curves determined by Eggen were superposed on the 
points to give the best apparent agreement. ‘This occurred when the modulus 
for the cluster was taken as 9™-53 and the colour scale was moved so that the 

51* 
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NGC 6025 results were bluer by 0-04 units than Eggen’s figures. Examination 
of this diagram (see also Table IV) shows that, allowing an extension of Eggen’s 
‘* bright blue dwarf ’’ sequence, as shown by a broken line, twenty-one stars fall 
in the system, and eleven fall outside its limits. Of these eleven stars it is possible 
some would be spectroscopic binaries and would, if their components were of 
appropriate magnitude, account for some of the anomalies. Holmberg mentions 
that some 10 out of 49 known cluster members of the Hyades are given as 
spectroscopic binaries. If this proportion is maintained in the present cluster 
then about 4 of the 11 observed anomalies might be accounted for in this manner. 
Accordingly the number of field stars to be accounted for would be between 
7 and 11, which figures are within the expected maximum of 12-5 derived above. 
It is therefore concluded that there is nothing in the present results to provide 
positive evidence against the concept that the actual members of the cluster 
NGC 6025 would give a colour-luminosity array which could be represented by 
the curves drawn by Eggen for the Pleiades. Further, certain of the results 
which fall notably away from the lines in Fig. 3 are from stars which have proper 
motions thought by Wood to be large enough to show non-membership of the 
cluster. These stars are as follows :— 


—59° 6549 and 6556 
— 60° 6296, 6308, 6327 and 6341 


They are to be regarded as non-members of the cluster both by photometric and 
proper-motion criteria. One star, —60° 6334, belongs photometrically to the 
cluster but has a relatively large proper motion. ‘The stars filling the gap 
between the blue dwarf and the bright blue dwarf lines (see Fig. 3) do not show 
proper motion, but, of course, whilst relatively large proper motion indicates 
non-membership of the cluster, the converse, that relatively small proper motion 
indicates membership, is not necessarily true. Finally the bright star — 60° 6348, 
which falls well away from the lines, is noted as having a spectrum which shows 
the lines Hf and Hy in emission, and this may account for its relatively blue 
colour. Twenty-seven stars, denoted by asterisks in Table III, are given as 
members both by Band FS criteria. These 27 stars include — 60° 6334 which is 
reported by Wood as having a relatively large proper motion. 

5. Distance and dimensions.—Collinder’s catalogue lists a variety of figures for 
the angular diameter of the cluster, ranging from 10’ given by Melotte to a 
figure of 19’ from Raab. ‘The present star-counting results suggest a figure of 12’ 
for the central region, but including the outer stars which fall on the lines in Fig. 2 
the cluster may be regarded as extending over an area about 28’ in diameter as far 
as stars down to about the twelfth magnitude are concerned or about 30’ to cover 
the stars listed in Table IV. 

Collinder, from an estimate of the integrated magnitude of the cluster, gives 
a figure of 1020 parsecs for the distance, neglecting absorption. Trumpler’s 
photographic modulus of 9™-5 from spectroscopic parallaxes corresponds, for 
zero absorption, to a distance of 800 parsecs, whilst Mrs Riecke’s modulus gives 
480 parsecs. The fine-structure interpretation of Fig. 3 suggests a modulus of 
g™-5(3), agreeing excellently with Trumpler’s result and corresponding to a 
distance of 805 parsecs. Adopting the present figures of distance and diameter 
leads to an estimate of some 6-2 parsecs for the overall dimensions of the cluster, 
with a possible reduction on account of absorption. 
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TABLE IV 
NGC 6025 interpreted by fine-structure hypothesis 
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Fic. 3.—Colour-magnitude diagram for the central region of the cluster. 
those given by Eggen for the Pleiades, taking the modulus for the cluster as 9™+53 and adjusting 
the colour scale by the addition of 0,040 to Eggen’s figures. 


The following table, which applies to stars brighter than 12™, for which the 
observations are relatively complete, summarizes the results of the “ fine struc- 
ture ’ (FS) and “ band ” (B) methods of examining the measurements. 


TABLE V 
NGC 6025 


Stars brighter than 12™ 


Angular diameter (minutes) 
Distance (zero absorption) parsecs 
Linear diameter (zero absorption) parsecs 
Cluster members 
Dwarfs 
Blue dwarfs 
Bright blue dwarfs 


Total stars 
Density (stars per parsec’) 
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Summary 

‘The methods applied in this paper have been described in previous 
communications from the Norman Lockyer Observatory. The distance 
modulus, based on CI, is used for the recognition of cluster members and the 
determination of the cluster distance. For the clusters NGC 6531, 6546, 
6494 corrections for absorption were obtained, after estimations of the CE. 
This was not possible for the clusters NGC 6469, 7127, 7128. The clusters 
NGC 6469 and 7127 are assumed to contain giants. For the NGC 6544 a 
general description is given. 





The methods of surveying the galactic clusters, described by G. Alter in 
previous contributions from the Norman Lockyer Observatory (1), have been 
applied to seven galactic clusters, with the results given in this paper. ‘The 
material is arranged as in the communications mentioned. 

The clusters of the present paper are located in parts of the sky including 
striking objects, such as conspicuous nebulae, obscuring matter and tufts of filmy 
light. 

The mean deviations of magnitudes of the NPS on the set of plates of the clusters 
NGC 6531, 6546 was +0™-12 on the photographic scale, and +0™-18 on the 
photovisual scale. ‘The mean error of the CI is therefore +0™-21. ‘The corre- 
sponding numbers for the set of plates of the clusters NGC 6469, 6494 were 
+o™-18, +0™-16 and +0™-22 respectively; while for the clusters NGC 7127, 
7128 they were + 0™-11, + 0™-13 and +0™-17 respectively. The order of exposure 
in all cases was: NPS—Cluster. 

NGC 6531 
R.A. 17"58™-6, Dec. —22° 30’ (1900-0) 

This cluster shows a strong concentration in the centre, where, on the Sidmouth 
plates, the stars are almost indistinguishable. Shapley (2), Raab (3), ‘Trumpler 
(4), Collinder (5), C. A. Ricke (6), Hayford (7) and Zug (8) have published results 
of investigations on this cluster. 

The star counting gave 109 stars as the total number within a cluster region of 
13'6 diameter. After correction for background the number of cluster stars is 57. 
The accompanying figure shows the star densities per square degree within the 
cluster region after correction for background. 

Table I gives the CE of some cluster stars after a comparison of CI computed 
from the spectral types, published by Zug and in the HD Extension (1949), with 
the CI found in this investigation. 
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In Table II the coordinates, the magnitudes and the distance moduli of the stars 
within the cluster area are shown. The distance modulus m—M was obtained 
from the CI corrected for the CE givenin TableI. For the other stars, of unknown 
spectral type, the assumption of a mean CE + 0-43 was made. 

The mean distance modulus of the cluster group is 9°34. This value, 
according to the formula 

m— M,=5logr+a,,.7—-5 (1) 


(where a,,,, =0™-35 per kiloparsec was assumed) corresponds to a distance 680 + 60 
parsecs. 





r 


8’ 10° 12” 





Fic. 1.—Star density per square degree above background density. 
I, NGC 6531. II, NGC 6546. III, NGC 6494. 


TaBLe I 

Colour excess (NGC 6531) 

CE Star CE 
+089 +010 

0°39 ; 0°39 

O'12 O'ol 

0°20 1°32 

o' 33 55 O*44 

1'07 K 1°00 


Therefore the cluster diameter is d= 13'-6 = 2°7 parsecs. 
The results of previous observations are 


Distance Diameter 
Raab 715 pe Ig-20’ 4:2 pe 
Shapley 570-910 10 26 
‘Trumpler 98o IZ 34 
Collinder 1590 10 8646 
Rieke 1350 


Of the 57 CI stars of Table II, 44 belong to the cluster and the same relation, 
when applied to 109 stars obtained by counting, gives a total of 84. We may 
correct the result in the following way. Out of 15 stars fainter than 12-00, 8 
belong to the cluster group and the same relation, when applied to 67 stars fainter 
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than 12-00, gives 36 cluster stars. When we add the 36 cluster stars, brighter than 
12-00, we get 72 cluster stars. Eleven stars of a nearby area of the cluster size are 
to be found lying in the (m— MM) cluster range. After applying this background 
correction, we have 61 as the final number of cluster stars. ‘This number agrees 
with the one obtained statistically. The dispersion iso= +8. ‘The magnitude 
distribution (after correction for background) is given in the following table: 


mM, 8 9 10 II 12 »E2°5 
Number of stars I 4 4 6 18 28 

A comparison of the magnitudes and CI of Table I] with that given by Zug 
showed appreciable differences. 


TABLE II 
Coordinates, magnitudes and distance moduli of 57 stars (NGC 6531) 


Zero point=star number 46, unit=1':0, X=east, Y=south. Coordinates of the assumed cluster 
centre X= —1':-6, Y=o-o 
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NGC 6546 
R.A. 18"01r™-2, Dec. —23°19' (1900-0) 
This cluster is a large one, but without appreciable concentration. It is 
mentioned in the catalogues by Trumpler (4) and Collinder (5). The star 
counting gave 107 stars within a region of 13’-6 diameter. Correction for the 
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background density results in 50 cluster stars. ‘The accompanying figure gives 
the star density per square degree within the cluster region. 

Table III shows the results of a comparison of the measured CI of 8 stars with 
the CI computed from the HD Extension (1949). The mean CE is +0™-12. 


TABLE III 
Colour excess (NGC 6546) 


Star CE Star CE 
60 —0'07 gI +0'10 
61 +O°13 96 —o'1g 
75 +O°15 100 + 0°29 
89 +O’ol 102 +051 


In Table IV the distance modulus m— M is given after the correction for CE 
givenin TableIII. Forthe stars which do not appear in Table III, the assumption 
of a mean CE +0™-12 was made. ‘The stars 60 and 8g, which are the brightest 
of the cluster and are also of late spectral types (M2 and Ko respectively), were 
regarded as giants and the value M =o was assumed for them. 


Tas_e IV 
Coordinates, magnitudes and distance moduli of 23 stars (NGC 6546) 


Zero point=star 96=HD 165554. Unit=1’-0, X=west, Y=north. Coordinates of the assumed 
cluster centre X= —2:°7, Y=+0°5 
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The mean distance modulus of the cluster group is 9-90 corresponding to 
r=830+60parsecs and d=13'-6= 3:2 parsecs. 
The results of the previous observations are 


Distance Diameter 
‘Trumpler 3000 pe 13’ Ill pe 
Collinder 1800 12 6-3 


Trumpler’s estimate was not a direct one. 

Among the 23 CI stars, 13 belong to the cluster and the same relation when 
applied to 107 stars, obtained by counting, gives a total of 61 stars. ‘This became 
59 after correction for background and agrees with the number of cluster stars 
obtained statistically. The dispersion iso= +8. 
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The magnitude distribution (after correction for background) is given in the 
tollowing table: 
mM, 9 >12°5 
Number of stars 2 48 


NGC 6469 
R.A. 17546™9, Dec. —22°19’ (1900-0) 


This cluster is fairly irregular and not well detached from the environs. _ It is 
mentioned by Shapley, Raab, Trumpler and Collinder in their papers. 

Table V gives the coordinates, magnitudes and distance moduli of stars within 
the cluster region. In the first part of the table the absolute magnitudes were 
taken from the main sequence. Here appear two ranges in the distance modulus 


TABLE VA 
Coordinates, magnitudes and distance moduli of 21 stars (NGC 6469) 


Zero point=star number 1=HD 162694. Unit=1’-0, X=west, Y=south. Coordinates of the 
assumed cluster centre X= —o0°8, Y=3°7 
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which indicate two cluster formations. ‘The mean distance modulus of the first 
group is 6-4, its distance 190 + 20 parsecs ; the corresponding values for the second 
group are, respectively, 10:5, 1200+ 130parsecs. But it is not probable that 
there are two clusters in exactly the same area; so the nearer group is assumed to be 
a group of giants. ‘This assumption is supported by the fact that in HD Extension 
two stars appeared with late spectral types ; it is corroborated also by Lundmark (5) 
who maintains that the cluster is a Praesepe cluster superposed on a Pleiades cluster. 
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Using the giant values of M for the stars of the first group, the second part of 
Table V gives the redistribution of m— M. Star 40 has been included in the giant 
group as it appears in HDEasaGostar. The contrary happens with the star 24, 
which appears as A. 

The mean value of m—M is 11-02, and consequently r = 1600 + 170 parsecs, 
d=12'-4=5°8 parsecs. cy 

In the case of this cluster, where no CE is available, the absorption has not been 
considered and the distances are computed according to the formula 


m—M,=5logr—5. 
The results of previous investigations were 


Distance Diameter 
Shapley 1820-2880 pe "Io pe 
Raab set vin 
‘Trumpler 2370 ; I0"4 
Collinder IIIO 3-6 
Lundmark 1560 45 


The surrounding of the cluster by a dark nebula suggests the necessity, in the 
case of the above distance, of correcting for absorption; but no spectral types are 
available. 

The star counting gave 164 stars within a diameter of 12’-4. The surrounding 
of the cluster by a dark nebula prevents a true background correction. The 
number of CI cluster stars is 20. The magnitude distribution in the cluster is the 
following : 

mM, II 12 13 
Number of stars 2 9 9 


NGC 6494 
R.A. 17251™-0, Dec. —19°0’ (1900-0) 

This cluster is a bright, large one. Results of investigations on this cluster 
have been published by Shapley, Raab, Trumpler, Collinder and Klauder— 
Lambrecht (9). 

The star counting gave the total number of stars in the cluster region as 333, 
and the cluster members as 149 within a cluster region of 27’-2 diameter. Our 
figure shows the star density per square degree within the cluster region. The 
density distribution given agrees quite well with that given by Klauder and 
Lambrecht. 

Table VI gives the coordinates, magnitudes and distance moduli of the stars 
within the cluster region. The cluster range contains 95 out of 165 CI stars. 
The mean value of m—M for the cluster group is 8-34. This, according to the 
formula (2), corresponds to a distance of 470 parsecs. 

Correction for absorption can be estimated for the 17 cluster stars for which the 
spectral types are available from HD Extension (1949). ‘Table VII gives both 
m—M,andCE. ‘Theaveragem— M is 9-63 and the distance, according to formula 
(1), is 770 parsecs. The mean CE is +0:52. No assumption can be made for the 
CE of the other stars. First, because the number of cluster stars with known CE 
is relatively too small and no investigation of the value of CE in the other stars is 
possible. Second, it is not possible to apply mean CE as the cluster covers quite 
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TaBLe VI 
Coordinates, magnitudes and distance moduli of 165 stars (NGC 6494) 


Zero point=star number 152=HD 163426. Unit=1’-0, X=west, Y=south. Coordinates of the 
assumed cluster centre X= —o'1, Y= —2°4 
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TaBLe VI (Continued) 
Coordinates, magnitudes and distance moduli of 165 stars (NGC 6494) 


Zero point=star number 152=HD 163426. Unit=1’-0, X=west, Y=south. Coordinates of the 
assumed cluster centre X= —o-1, Y= —2-4 
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Tasie VII 
Distance modulus and colour excess (NGC 6494) 
Star m—M CE Star m—M 


122 8- +0°33 124 9°7 
100 8-: +O'4I 55 9°9 
152 8:5 +0°53 130 10°0 
168 9° + 0°57 163 10°0 
82 9°3 +0°25 164 10°0 
153 9°4 + 0°60 125 10°5 
103 9°5 +O°41 225 10°5 
63 9°6 +068 233 It'l 
161 9°6 +0°66 


a big area in which sudden changes of absorption are obvious even at a first glance 
at the area. ‘This is also borne out by Klauder’s and Lambrecht’s investigation of 
the density distribution. 
‘Therefore we get from Tables VI and VII respectively : 
r=470+40parsecs, d=27':2=3-7 parsecs 
and (corrected for absorption) 


r=770+50parsecs, d=27':2=6rI parsecs. 
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The results of previous investigations were 


Distance Diameter 
Shapley 870-1380 pe 5° 10-0pe 
Raab 400 E 5°5 


Trumpler 660 5+ 
Collinder 957 


The CI method gives a different result for the star number from the statistical 
method. Among the 165 CI stars, 95 belong to the cluster modulus range, and the 
same relation, when applied to the 333 stars obtained by counting, gives, after 
correction, 169 cluster stars. ‘This number is different from the number 149 
obtained statistically. ‘The dispersion is c= +13. The difference is due to 
deviations in CI. As the stars in the area of this particular cluster appear in 
chains, some fainter stars located in the haloes of bright stars appear brighter than 
they are. ‘The modulus range, and hence the star number, depends on a statistical 
behaviour, so that a few faint stars having larger deviations are not included within 
the right range of m—M in Table VI. 

The magnitude distribution in this cluster, after correction for background, is 
in the following table : 

m 10 II 12 : >13°5 


v 


Number of stars 12 17 24 g6 


NGC 6544 
R.A. 18"01™-2, Dec. —25° 01’ (1900-0) 

This cluster is a very faint object mentioned in Shapley’s catalogue of open 
clusters, but, according to Collinder (5), it ‘‘seems to be a globular cluster or 
anagalactic nebula”. NGC 6544 appears on Sidmouth plates as a very compact 
conglomeration, not resolved, of very faint stars. The unresolved central part is 
of about 1-5 diameter. Around that there are about seven stars in a diameter of 
about 2’°7. —————— 

NGC 7127 
R.A, 21° 40™5, Dec. +54° 09" (1900-0) 

This cluster is a small and poor one, mentioned in Shapley’s catalogue. It is 
located in the middle of a small area, extremely poor in stars. On account of its 
particular situation the star counting was not done in the usual way, but a compari- 
son of the cluster’s region with some neighbouring regions of the same size was 
made. Within a diameter 2’:8 of the cluster’s region there are 15 stars. After 
the correction the number of cluster stars is 8. 

Table VIII gives the distribution of m— M and it is given in two parts, of which 
the first shows the stars separated in two ranges, and the second is a rearrangement 
of the first, treating the first group of stars as giants (as was done in the case of 
NGC 6469). In this table we see that the giants fit well into the dwarf range. 

Table VIII shows that the mean distance modulus is 10-28 in case A and 10°55 
in case B corresponding to 

Case A r=1120+120pe, d= 
Case B r=1300+120pe, d= 


‘S=1-0pe. 
‘S=IIpe. 


>/ 
2 
2! 
z 


Since spectral types are not available it is not possible to suggest a correction 
for the interstellar absorption. 
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The data in Shapley’s catalogue are 

r= 1820-2880 parsecs, d=1'-0=0°8 parsecs. 

The number of cluster stars is 10 (among them 3 giants), against 8 obtained by 
counting. 

The magnitude distribution is given in the following table : 

m,, 10 II 12 13 14 15 
Number of stars I I 5 I I 


Tasie VIII 
Coordinates, magnitudes and distance moduli of 15 stars (NGC 7127) 


Zero point=star 13. Unit=1’-o, X=east, Y=north. Coordinates of the assumed cluster centre 
X=—o'5, Y=+04 
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NGC 7128 
R.A. 21"40™7, Dec. +53°15' (1900-0) 


This cluster is a small one mentioned in the catalogues of Shapley, Trumpler 
and Collinder. It is located within a small area of low star density, like the 
previous cluster. So the star counting and the correction for background were 
done in the same way as with NGC 7127. Within a diameter 3'-1 of the cluster’s 
region there are 16 stars. After correction the number of cluster stars is 9. In 
Table IX we find the mean distance modulus of the cluster group (without 
correction for absorption) to be 10-99, corresponding to 

r=1570+170parsecs, d=3':I =I-5 parsecs, 
The results of the previous investigations were 
Distance Diameter 
Shapley 1820-2880 pe 2" 7 pe 
‘Trumpler 3050 2 34 
Collinder 6670 "3 594 


Trumpler’s distance was not measured but computed. 
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The number of cluster stars after reduction is 9, coinciding with the number 
obtained by counting. 
The magnitude distribution is given in the following table: 
m.,. 12 13 
Number of stars I 2 


TaBLe IX 
Coordinates, magnitudes and distance moduli of 14 stars (NGC 7128) 


Zero point=star 15. Unit=1’-0, X=east, Y=north. Coordinates of the assumed cluster centre 
X=—o8, Y=+06 
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papers which proved most helpful to me. I am grateful to Professor F. J. M. 
Stratton and the I.A.U. for the grant which enabled me to carry out my work at 
the Norman Lockyer Observatory. I should like also to express my thanks to 
Mr D. L. Edwards, Director of the Norman Lockyer Observatory, and to Mr D. R. 
Barber for their interest and helpful discussions throughout the course of my work 
at the Norman Lockyer Observatory. 


Norman Lockyer Observatory, 
Sidmouth : 
1953 April 23. 
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SOME REMARKS ON THE INTERPRETATION OF 
APSIDAL-MOTION CONSTANTS IN CLOSE BINARY SYSTEMS 


Zdenék Kopal 


(Received 1953 October 26) 


Summary 


The aim of the present paper has been to establish explicit formulae for 
the computation of the constants k, of apsidal motions in close binary systems, 
invoked by the j-th harmonic distortion of their components of arbitrary 
structure. If the degree of central condensation of such components is high, 
and their (equilibrium) density distribution p(a) remains analytic at the 
surface a=a,, it is shown that 


_ j+1—7;,a) 


Je 2j+2n,(a;) ; 


where the quantity 7,(a,) is expressible by a series on the right-hand side 
of equation (14) within the domain of its convergence (the existence of the 
latter defining the ‘‘ stars of high condensation ’’). 

If the density distribution characterizing our stars is such that not all 
derivatives of p(a) are finite on the surface, it is shown that the values of k; 
may be effectively approximated by means of the equation 


(2j+ 1)a,"/* *k; = 3(j+2) i Da® da, 


where D= p/p, p(a) denoting the mean density interior to a. The derivation 
of this formula reveals that the error of this representation of k; should 
diminish with increasing value of 7 as well as with increasing degree of 
central condensation. 

A numerical application to the polytropic family of models discloses that, 
if the ratio of the central to the mean density of our configuration exceeds 
(say) one hundred, the above approximate formula should furnish the values 
of k, correctly within the limits of errors affecting their empirical determina- 
tion in most practical cases, while the question of errors of our approximate 
values of ky or ky will scarcely ever arise. 





Ever since Russell’s pioneer work on the subject (1) it has been generally 
recognized that the observed rate of apsidal advance in close binary systems 
exhibiting eccentric orbits represents a royal road to certain empirical 
information concerning the internal structure of the constituent components. 
In more precise terms, the rate of apsidal advance arising from j-th harmonic 
distortion of either component is known to be proportional (to the first order in 
small quantities) to the constant 


oe I —7;(4;) 
*  2g+2nf(a,) ’ 


(1) 
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where 7,(a,) denotes the surface (a= a,) value of a particular solution of Radau’s 
differential equation 


d nica 
as! +n(n-1)+ 6D +1)=i(j +9), (2) 


satisfying the initial condition 

(0)=j—2, (3) 
where D=p/p, p(a) denoting the density of the respective undistorted star at a 
distance a from its centre, while p(a) stands for the mean density inside a 
volume of radius a. The j-th harmonic distortion of figure of the respective 
components due to its own axial rotation and to the tidal effect of its mate— 
which causes the apsidal advance—can, in turn, be shown (2) to be proportional to 

2 

eee (4) 
J +741) 
The values of k; (or A;) in specific eclipsing systems can be established by an 
approximate analysis of the observational data (3); and their interpretation 
represents so far the only available empirical check of the theories of the internal 
structure of main-sequence stars. 

The foregoing definition of the constant k; should enable us to predict the 
rate of apsidal advance in close binary systems caused by mutual distortion of 
their components—at the price of a numerical integration of equation (2)— 
provided that the equilibrium density distribution p(a) of the distorted bodies 
is a priori known. It throws, however, very little light on the following converse 
question of considerable astrophysical interest: namely, what is the exact amount 
of information an empirical knowledge of k, entitles us to draw concerning the internal 
constitution of the respective stars ? ‘This question will be the subject of our present 
discussion. 

If we regard this problem from a completely general point of view, any attempt 
at a comprehensive answer encounters serious difficulties. In special cases, 
however, our problem can be solved in a closed form. Thus, if the function D(a) 
in (2) is set equal to zero (which would correspond to an infinite density 
concentration), equation (2) can be solved to yield 

n(a)=j+I, (5) 
which, in turn, leads to the well-known results of k; =o and A;=1. Moreover— 
and this is a fact of salient importance—numerical integrations of (2) have 
revealed that, for models of high central condensation, the values of k;—though 
finite—approach zero rapidly; and the values of A; tend asymptotically to one. 
Since the density concentrations of actual stars are very likely high, the whole 
astrophysical interest is largely confined to this particular case. 

The point of departure of the present investigation will, therefore, be a 
search for the reason why the values of k; as defined by equation (1) diminish 
so rapidly with increasing degree of central condensation. In doing so we shall 
be led to construct new asymptotic relations between k,; and p(a), which will 
not only enable us to compute the values of k,; for centrally condensed stars in a 
more expeditious manner (i.e. by simple quadratures rather than by way of 
numerical integration of a differential equation), but will also render the 
dependence of k; on p(a) much more transparent than it was in our original 
formulation of the problem. 
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In order to do so, let us depart from Radau’s equation (2) which is quadratic 
in the dependent variable. One way of removing this non-linearity is to replace 
n(a) by a new variable y(a), defined by 


_ ady p 
"= yda_ 5’ (6) 
which transforms (2) in 
d*y — fIF+t) 
i a +a(a)} (7) 
where we have abbreviated 
2(dp/d 34 
se)= LES). 32. (8) 


pa"da 


~0 
In order to solve this equation, let us set 
y(a) =a’*! exp [ u(a) da, (9) 


where u(a) denotes a new variable defined by a non-homogeneous differential 


equation 


du K 
SF outs. Smee 
qq tht a g(a) (10) 


of the Riccatti type, where we have abbreviated 
K=2(j+1). (11) 


Suppose that we assume now a solution of (10) to exist in the form 





2 u,(a) 
u(a)= X = (12) 
i ) K' 1 
of a series in descending powers of A, an insertion of (12) in (10) and the fact 
that the ensuing series must vanish identically determines the successive values 
of u,(a) as: 


Up = ag, } 
u,= —atly , | 
, 
Uy = —a{u,' +u,"}, t (13) 
" i-1 
u; = ee D> uth 
\ k=1 : J 


in terms of the function g(a) depending on the structure of the respective 
configuration, and of its successive derivatives. If we insert now (12) in (9), 
differentiate the latter logarithmically and insert in (6), we find that, on the surface, 


: a ua,) 
a,)=j+1+ + ———7—|. 
nf 1) Pi ; 9 2th +1) tt 


(14) 

For stellar configurations of high central condensation, the quantity g(a) as 
well as its derivatives are likely to be very small for a=a,, and tend in fact to 
zero as the condensation becomes infinite. Moreover, the quantity K whose 
powers constitute the denominators of successive terms on the right-hand side 
of (14) is numerically equal to 6, 8 and 10 for the second, third and fourth 
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spherical harmonic distortion, respectively—a fact which should assist the 
convergence of our expansion. ‘This is why, for centrally condensed configura- 
tions, »,(a,) departs so little from its limiting value of 7+1—the less so, the 
greater j—and, consequently, why the constants k; are numerically as small 
as was heuristically established by previous investigators.* 

An application of the foregoing equation (14) to the actual calculation of k; is, 
however, inseparably connected with the problem of convergence of a sequence 
of the values of |u,(a,)|._ Whether or not such values—small as they are— 
diminish actually with increasing 7 so as to render the infinite series on the 
right-hand side of (14) convergent depends, of course, on the particular 
properties of the respective stellar models, and cannot be answered a priori in 
full generality. For some models this convergence may indeed be so rapid that 
the series on the right-hand side of (14) may be terminated after a very few terms. 
Other models may, however, possess such special features as to make a solution 
of (10) of the form (12) actually non-existent. This is, for instance, true of the 
well-known polytropic family of models whose successive density derivatives at 
the surface are either zero or infinite.t In the face of such a situation the need 
evidently arises for broadening the basis of our analytical approach to the 
problem at issue so as to encompass even such recalcitrant cases. 

In order to do so, let us return to Radau’s equation (2) and inquire as to the 
possibility of constructing such solutions of it, which, in addition to fulfilling the 
initial condition (3), is required to assume at the surface (a=a,) a specific value 
n(a,)=A (say). If we revert from 7(a) to the variable y(a) introduced by the 
transformation (6), this is tantamount to a search for characteristic solutions 
of the following linear boundary-value problem: 


L\y|=g(a)y, (15.0) 
y(0)=0, —-y'(a,) =Ay(44), (15.1) 
where 
d? j(j+1) 
in? or (16) 


and A is a prescribed parameter. 
The operator L is evidently self-adjoint. Its Green’s function appropriate 
for the given boundary conditions will, therefore, be symmetrical and of the form 


I 





A/ yak J j Y £3+1),3+1 > 
G(a, €)= +1 {E97 4+ C€'*}a (€ >a), (17) 
in which we have abbreviated 
aZHC, = J+ ayr_ (18) 


: jJ- a,A + I 2 
which, by (1), can be seen to reduce to 


Y—-1 = 2q 2+] 
C;*=2a,”)*1k,. 





* For a recent comprehensive summary of numerical determination of apsidal-motion constants 
for different stellar models, cf., e.g., L. Motz, Ap. ¥., 115, 562, 1952. 

t As is well known, for complete polytropes p=ppy", where y denotes the Emden function 
corresponding to a polytropic index n. In consequence, the k-th derivative of p will contain a 
term multiplied for 4"—*); and as y vanishes on the surface y("—*)(a,) will become zero or infinite 
depending on whether n 2 k. 
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With the aid of the foregoing Green’s function the differential boundary-value 
problem may be rewritten in the form of the following integral equation, 
revealing that 


(2) + 1)y(a)= —a’*IC, |. a(E)y(6)E'*? dé 
<a! | g(e)y(e)e'* dé 


— ai |" gy db. (19) 


Now we have already seen that, for configurations of high central condensation, 
k; is apt to become a very small quantity. If so, however, then the constant C; 
(being inversely proportional to k;) must necessarily be large. The three integrals 
on the right-hand side of (19) are all quantities of the same order of magnitude ; but 
only the first one is multiplied by C;. Therefore, to a first approximation, 
equation (19) should be essentially equivalent to 


“ay 


(2j+ 1)y(E)= — a AC, | gE) (ENE 4; (20) 


but since the limits of the integral on the right-hand side are now constant, the 
solution of this equation is clearly of the form 


y(a)= Aa’, (21) 


where 4 is a constant. Since, however, equation (20) is homogeneous in the 
dependent variable, its validity requires that the value of C, involved in it be, 
not arbitrary, but given by the equation 


v—] _. 27+1b — 
C,t=24a,3+1k, = — 





“ay dp 
3 | I op £25+1 dé, (22) 


2j+1)9 pda” 


which thus expresses the desired parameter &; in terms of a simple quadrature of 
the ratio p’/p. 

It may be noted that the integral on the right-hand side admits of the 
foregoing simplification. Since the mean density p(a) inside a sphere of radius a 
is defined by 


a | ye 
= a pa” da, (23) 
J0 


so that 


a dp 
PP i (24) 


a resort to partial integration reveals that 


“Id, — . oe 
; ae da= 3| ‘ " Dia da—2(j +2) | _ Da da. (25) 
For centrally condensed configurations D <1 and, therefore, D?<D. In such 
cases, the first integral on the right-hand side of the foregoing equation is likely 
to be very small in comparison with the second and may, therefore, be neglected. 
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Doing so we arrive at last at the final result of this investigation in the form 
- 
(2) + 1)a,”“4k,=3(j+2) | Da® da, (26) 
“0 


expressing the apsidal-motion constants k; of centrally condensed configurations 
in terms of a quadrature of the ratio D=p/p, weighted by a®’. It transpires, 
therefore, that for such configurations the values of k; are primarily determined 
by the behaviour of the ratio D throughout the interval 0 <a<a,, with the 
contribution of the outer parts of stellar interiors becoming relatively the more 
important the higher the value of /. 

It is interesting, in this connection, to recall that Motz (4) has recently 
examined a suggestion (credited to Russell) that the apsidal-motion constant k, 
is determined essentially by the radius of gyration H of our configuration which 
is known to be defined by 

8a (% 


H?= om], pa‘ da, (27) 


where denotes the total mass of the configuration. On the basis of a series of 
numerical integrations Motz has established the existence of an empirical 
relation between k, and H/a, which appeared to be quasi-linear in their 
logarithms ; but some models deviated conspicuously from it. The outcome of 
our investigation throws some more general light on the question of the 
dependence of k; on the density distribution of the respective models; and 
discloses, in particular, that, to a first approximation, k; should be proportional 
to an integral of (p/p)a®’—rather than of pa®’—taken between the limits (0, a,). 

In order to test the performance of our approximate formula (26) in a 
practical case, we set out to approximate with its aid the apsidal-motion 
constants k,, for j=2, 3, 4, appropriate for the polytropic gas spheres 
characterized by the index »=3 and 4. In this case equation (26) assumes the 
particular form 


(2j+1)x,24k,=(j+2) | Ux dx, (28) 
where x, stands for the first zero of the corresponding Emden function y(x) and 
ae aa (29) 


denotes a function which has been tabulated by Comrie and Sadler (5). 
The outcome of our computation is summarized in the following table, in 
which the constants k; approximated with the aid of (28) are compared with those 


k; 
n=3 n=4 

_ ©°'01207 0-001 100 
iin (00144) (0-001 19*) 
7 0003524 0°0002319 
ee (000368) (0:000235) 

. 0°001387 0:0000653 
j=4 (0:00140) (0:00007) 


* See footnote on p. 775. 
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(in parentheses) obtained by Chandrasekhar (6) on the basis of numerical integra- 
tion of a differential equation equivalent to Radau’s.* Considering the simplicity 
of (26), the degree of approximation attained with its aid appears to be satisfactory 
for all practical purposes (for one can very seldom establish the observed value 
of k; within an error of less than 5—10 per cent), and exhibits a marked improve- 
ment with increasing values of both 7 and n—1.e. precisely in such cases when a 
numerical integration of Radau’s equation becomes most laborious. 


Department of Astronomy, 
University of Manchester: 
1953 October 26. 
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PHOTOELECTRIC OBSERVATIONS OF OCCULTATIONS AT 
THE CAPE OBSERVATORY 


A. W. F. Cousins and R. Guelke 


(Communicated by H.M. Astronomer) 


(Received 1953 December 1)* 


Summary 


A simple apparatus for observing occultations photoelectrically is 
described. The curves obtained for the occultations of Antares on 1950 
May 4, 1951 July 15, 1952 April 13 and 1953 March 8 are given, reduced toa 
standard form in which the ordinates represent percentages of the 
unobstructed starlight, and the abscissae angular measure across the star. 





Introduction.—In 1949 Professor R. O. Redman and Dr A. D. Thackeray 
called the attention of Dr Stoy to the series of occultations of Antares, then 
about to begin, and suggested that they should be observed in an attempt to 
determine the diameter of the star. It was doubtful whether a refractor of 
24 inches aperture would be large enough to give useful results but the necessary 
subsidiary apparatus was designed and built in time for the first occultation on 
1950 May 4. The observation made on that occasion has been described 
elsewhere (1). Its value is considerably impaired by the presence of a parasitic 
50-cycle ripple on the trace. Observations free from this defect were obtained 
on 1951 July 15, 1952 April 13 and 1953 March 8. It was not possible to observe 
the other occultations of this series because of cloud. 

The photometer is of simple construction. A metal plate with a series of 
interchangeable diaphragms replaces the usual plateholder on the 24-inch 
Victoria Telescope. A photomultiplier tube and amplifier are enclosed in a 
metal box which is hinged to the telescope so that it can be swung aside while the 
cross-wires of the guiding telescope are being adjusted to correspond accurately 
with the centre of the diaphragm. The R.C.A. 931-A photomultiplier tube has 
its own light-tight housing within the box and is provided with a shutter and a 
Fabry lens. A voltage divider for the multiplier is mounted on the tube socket. 

The amplifier circuit is illustrated in Fig. 1. The sensitivity control in the 
anode circuit of the multiplier is variable in steps from 20 megohms downwards. 
A potentiometer provides grid bias for balancing the moonlight and dark current 
from the multiplier tube. No other controls are necessary. 

The photometer was originally operated entirely from batteries, but these 
have been largely replaced by a mains-operated power pack with regulated output 
voltages. A 6-volt storage battery still provides current for the amplifier heater 
and an H.T. battery is connected between the output cable and the Y-plates of 
the cathode ray oscillograph used as a recorder. The purpose of this latter 
battery is to centre the spot on the oscillograph screen. It has been found best to 


use a comparatively low voltage, between 70 and 80 volts per stage, on the 
multiplier tube. 
* Received in original form 1952 December 1. 
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Three twin beam oscillographs have been tried, and in every case it has been 
found necessary to connect the output from the amplifier tube directly to the 
Y-plates of the cathode ray tube to avoid mains interference. Even this is not 
always effective, and, as has already been remarked, there was considerable 
interference on the trace of 1950 May 4. These instruments allow a timing wave 
and clock pulses to be recorded as well as the star trace, and this feature has led 
to the routine use of the instrument for timing ordinary occultations (2, 3). 
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Most of the light reaching the photomultiplier tube is scattered moonlight. 
Only when the moon is very young or the star unusually bright is the starlight of 
comparable intensity. ‘The unwanted light can be reduced by using a smaller 
diaphragm, but the size cannot be reduced too far or vignetting of the image will 
aggravate the already troublesome seeing “noise”. At the Cape it has seemed 
best to use a comparatively large diaphragm. For the reappearance of 
1950 May 4 a diaphragm with an aperture of 2mm was used. For subsequent 
occultations the aperture was I mm (corresponding to 30") in diameter. With 
a reflecting telescope and good seeing, smaller sizes could undoubtedly be used. 

If the objective were perfectly clean, the scattered light would be mainly from 
the sky, but in the case of the Victoria Telescope much the greater portion originates 
in the tarnished surfaces of the comparatively old objective. Any variation in 
the moonlight due to clouds or haze has an immediate effect on the intensity of 
this scattered light. A milliammeter in series with the amplifier valve enables 
the observer to see what is happening and to make the necessary adjustments to 
prevent saturation. 

It is possible to check the sensitivity by moving the telescope slightly and 
“occulting”’ the star by the edge of the diaphragm. Such a check also serves to 
verify the centring of the star. 

The regular (8-inch) guiding telescope serves admirably for disappearances: 
for reappearances guiding presents something of a problem. This was solved 
temporarily for the reappearances of Antares by fitting a small telescope that 
could be inclined sufficiently to guide on a nearby bright star (e.g. o or 7 Sco) 
that was itself not occulted. 
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The Antares observations.—The circumstances of the four occultations 
observed at the Cape were as follows: 


Angle between 


Date Type “8° GMT. PA. normmaltolimbend ,,**** 
Moon : : ("/sec) 
Moon’s motion : 
d h m ° ° 

1950 May 4 RD 16°7 2 21°0 267 14 O'414 
1951 July 15 DD I1°4 16 42°7 43 68 o'1s1 
1952 April 13 RD 19'1 22 16:9 301 12 0°459 
1953 March 8 RD 22°0 o1 56:2 235 52 0°233 


The records obtained on the second and later occasions are shown in Plate 11. 
A record of an occultation of A Scorpii (4™-8 B3), a star of presumably negligible 
diameter, is given for comparison. The timing trace consists of a ripple at either 
50 or 100 cycles per second punctuated once a second by a signal from the standard 
clock. In the reductions allowance was made for the fact that the time and star 
traces are not quite vertically above each other on the film and that the camera 
drive was not always uniform. The response of the equipment was not linear 
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Fic, 2.—Calibration curves: A, 1951 July 15; B, 1952 April 13; C, 1953 March 8. 
(Thick portions actually employed in reductions.) 


and the necessary correction was determined by the following method. The light 
of the moonlit sky was used to provide a source of constant intensity and a 
rotating sector was run in front of the diaphragm to produce a rapid alternation 
of light and darkness on the sensitive surface of the photomultiplier tube. The 
bias was then adjusted so that this constant signal was recorded at different 
heights on the oscillograph screen, the amplitude of the approximately square 
wave so produced being a measure of the sensitivity of the apparatus at that 
particular height. From these results it is possible to deduce a curve of 
deflection against incident intensity. The non-linearity is almost entirely due to 
the amplifier. The results are shown in Fig. 2. It is possible that some portions 
of these curves may be in error by one or two tenths of a millimetre. 

No corrections have been applied to the measures of the trace of the occultation 
of 1950 May 4 which were made on the original film with a long screw measuring 
machine. No calibration is available but as the deflection was small and not too 
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near the saturation point, it is believed that any distortion due to non-linearity 
of response will be unimportant. 

The observed curves reduced to a standard form of representation in which 
the ordinates represent percentages of unobstructed starlight, and the abscissae 
angular measure across the star, are illustrated in Fig. 3. 

The overall time constant of the apparatus was about one millisecond. The 
speed of response is therefore adequate to avoid distortion from this cause. 
The trace of A Scorpii (Plate 11) confirms this. 
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(a) Occultation of Antares, 1950 May 4. (6) Occultation of Antares, 1951 July 15. 
(c) Occultation of Antares, 1952 April 13. (d) Occultation of Antares, 1953 March 8. 

The short vertical lines indicate the standard error cf a plotted point of the corresponding 
intensity. 


The most serious cause of distortion of the curves is the seeing. Seeing noise 
is directly proportional to the amount of starlight. It is therefore most 
conspicuous when the star is clear of the Moon and is absent when the star is 
occulted. Variations in the scattered light are normally too slow to produce 
distortion. Shot noise merely broadens the trace. 

The seeing noise has been determined by measuring the trace from the 
unocculted star at intervals of 0-or second and computing the r.m.s. variations 
over periods of 0-1 second. Slower variations have little effect in distorting the 
curve. Peaks occasionally reach, but rarely exceed, twice the r.m.s. value. The 
uncertainty due to known causes and expressed as the standard error of a plotted 
point is indicated by short vertical lines in the figures. No satisfactory criterion 
is available for the occultation of 1950 May 4. 
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The effective diameters deduced from the curves on the assumption of a 
uniformly illuminated circular star disk and a smooth occulting edge are as 
follows : 


1950 May 4 0-045 
1951 July 15 0”:035 
1952 April 13 0”:045 
1953 March 8 0”-040 


The thanks of the authors are due to Dr R. H. Stoy for his continued interest 
in the observations and for his help in carrying them out, to Mr B. Offen who 
constructed the mechanical parts of the photometer and to Dr David S. Evans 
for helpful discussions. 


Royal Observatory, Dept. of Electrical Engineering, 
Cape of Good Hope: University of Cape Town. 
1953 November 20. 
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Occultation traces observed at the ¢ ‘ape 
(1) Antares, 1951 July 15. D.D. 
(2) Antares, 1952 April 13. RD. 
(3) Antares, 1953 March 8. RD. 
(4) 4 Sco, 1951 September 7. D.D. 
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Occultations of Antares as observed at Pretoria. 
(1) Antares, 1952 April 13. Top line, 1 second pulses; centre line, occultation trace; 
bottom line, 100 ¢ sec timing trace. 
(2) Antares B, 1952 April 13. Same scale as (1). 
(3) Antares B, 1953 March 8. 
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OBSERVATIONS OF OCCULTATIONS OF ANTARES 
WITH THE RADCLIFFE REFLECTOR 


David S. Evans*, ¥. C. R. Heydenrycht and ¥. D. N. van Wyk 
eyaenry J 
(Communicated by H.M. Astronomer at the Cape) 
(Received 1953 December 21) 


Summary 
The occultation of Antares on 1952 April 13 was observed photoelectrically 
and traces for both components were obtained. At the occultation of 
1953 March 8 only the trace for the companion was obtained. 





A special prediction of the circumstances of the occultation of Antares on 
1952 April 13 was kindly supplied by H.M. Nautical Almanac Office. These 
were: 

Disappearance: 21" 21™078 (U.T.) at P.A. 52°-4. 

Reappearance: 21"57™ 348 (U.T.) at P.A. 349°2. 

Age of Moon: 19:1 days. 

Time interval between passages of terminator and dark limb: 4 min. 
Star 116 B Sco, 9’ from dark limb in P.A. 68° at moment of reappearance. 


The reappearance was observed using the star 116 BSco (HD 148760) as a 
guide star. 

The photocell used was an E.M.I. tube type 5060 placed at the Newtonian 
focus of the Radcliffe reflector. A diaphragm was placed in the focal plane and 
the cone of rays was allowed to expand so as to form a circular patch on the 
cathode. ‘The tube was operated at 100 volts per stage and the maximum 
collector current drawn was just under 70 micro-amperes. These figures are 
both far below the ordinary safe operating limits. The dynode supply had a 
stability of better than one part per thousand for periods of up to 10 minutes, 
with a ripple of less than 10 millivolts. The signal was brought to the amplifier 
by a coaxial cable of low capacity, about 12 micro-micro-farads per foot. The 
amplifier circuit with built-in vacuum tube voltmeter is shown in Fig. 1. 

A Hathaway Type S8—D.B. oscillograph was used for recording. Three 
galvanometers recording on bromide paper 10 inches wide driven at 10 inches 
per second were used. One recorded the signal output from the photomultiplier 
tube and amplifier: one carried a series of one second pips from a Muirhead 
phonic Motor Clock driven from a 1000 c/s tuning fork oscillator: the third 
recorded a sinusoidal voltage of 100 c/s from a Hewlett Packard Audio Oscillator. 
The internal time marker of the Hathaway recorder is inaccurate and the marks 
were used only for convenience in measuring the trace 

The galvanometer used for recording the occultatio. was of the moving coil 
type with a d.c. sensitivity of 6-4 cm/milliampere/metre and a natural frequency 
of 580 c/s. The sensitivity in the recorder was 2-3 cm/milliampere. For the 


* Royal Observatory, Cape of Good Hope. 
+ National Physical Laboratory, S.A. Council for Scientific and Industrial Research, 
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4“ 
occultation the amplifier was adjusted for a bias of — 5-0 volts on each tube and 
125 volts on the anode of the photomultiplier tube. The d.c. supply voltage 
on the amplifier was 320 volts and the plate current from each tube 21-5 milli- 
amperes. ‘The coarse sensitivity control was set in the 200K (R;) position. 
Under these conditions the d.c. response was found to be linear up to 15 milli- 
amperes into the galvanometer. At minimum gain the current gain was found 
to be 0-65 milliamperes;volt. The theoretical analysis of the circuit yielded 
the figure of 0-67 milliamperes volt. 

Incident 

a 
Ad 


























H.T. 
1100-!700V. 


Fic. 1. 


Resistors 

R,=50 K2; R,=100 KQ, wire-wound potentiometer; Rz=10 KQ wire-wound; Ry=50 KQ 

wire-wound; R,=200K2 wire-wound; Re=1 K2 = wire-wound potentiometer, pre-set; 

R,=Rg=5 K2 wire-wound; Ry=1 K2 wire-wound potentiometer; Ryy=Ry=5'5 K2 

wire-wound; R,,=200 K2 wire-wound; Ri3=50 K2 wire-wound potentiometer; Ryy=75 KQ 

wire-wound; Ri;=R,;=10 KQ selected; Ryg=2 KQ potentiometer; R\z=2 KQ potentiometer, 
pre-set; Ryg=Ryy=12 KQ; Ry =215 KQ; Roe=5 K2Q; Ro3=100 KQ. 

Meters 
M,: 0-100 pA; M,, M3: microammeters, not used; My, M,: 0-530 mA; M,: o-5 mA. 

Supply 

XX: 6-volt filament supply. 

Minimum gain was used during the occultation. For the same settings the 
linearity test curve of Fig. 2(a) was obtained. The frequency response curve for 
the whole installation is shown in Fig. 2 (8). 

At the occultation everything went smoothly. Antares was about 3 hours 
east of the meridian at about 45° altitude. The seeing was poor. ‘lhe sky was 
partly cloudy, but the phenomenon took place in a gap in the clouds. ‘There 
may have been a little haze or high cirrus over the position at the time, but if 
so, it was very thin. 
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The trace shows both components. At this occultation the track of the star 
relative to the lunar limb made an angle of 58°-4 to the normal at the point of 
emergence. ‘The rate of advance of the limb parallel to itself was taken as 
o”-241/second. The effective wave-length of response for the equipment as a 
whole on the blue companion was estimated as approximately 4oooA. ‘The 
measures of the trace for the companion are compared with the theoretical 
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Fic. 2 (a).—Voltage increment on signal grid of amplifier. 
(Unit = 0-42 volts.) 
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diffraction pattern for a point source at this wave-length in Fig. 3. ‘lhe 
interval between the times of reappearance of the two components was 3°59. 
This result, together with the corresponding intervals for the 1953 observations 
at the Cape and at the Union Observatory, gives a separation between the 
components of 3”-16 and a position angle of 274° 38’. The residuals are :— 


” 


Pretoria, 1952 0022. 
j Cape, 1953 0:036. 
Johannesburg, 1953 0°025. 
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Curvature of the limb -~as neglected and the residuals are the residuals in time 
converted to angle at the rate of advance of the lunar limb parallel to itself 
appropriate to each occultation. 

The measures of the trace for the primary component are shown in Fig. 4. 
Seeing fluctuations up to 7:5 per cent (3°75 per cent on either side of the mean) 
are observed to affect the unobscured light of the star. Seeing fluctuation 
amplitudes are proportional to intensity. Hence if the figure of 7-5 per cent 
be taken as the maximum for the unobscured light, the corresponding maximum 
displacements due to seeing at smaller deflection can readily be computed. 
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Fic. 3.— Measured points compared with the theoretical diffraction pattern for a puint source. 
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kc. 4.—The vertical lines enclosed in the circles show the maximum possible displacement 
due to seeing at the given level. 


These values are shown by the lengths of the vertical lines within the circles 
in Fig. 4. It will be noticed that, even if at each level the observed points 
are displaced by these amounts, the curve corresponding to the passage of a 
straight edge over a uniform circular disk cannot be produced. The best fit 
to the “‘straight-line”’ part of the curve gives a diameter of 0”-042. This curve 
should be compared with that obtained at the Cape on the same occasion. 
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At the occultation of 1953 March 8 apparatus which was largely identical 
with that used on the previous occasion was employed. Owing to an unfortunate 
accident only the trace for the companion was obtained. The measures are 
compared with the theoretical diffraction pattern in Fig. 5. The traces are 
all reproduced in Plate 12. 

At this point it may be convenient to make a comment on the trace obtained 
with the Radcliffe reflector at the occultation of 1950 June.* It is now suspected 
that, between the last check on this point and the occultation itself, the build-up 
of lunar background light may have produced a slight saturation effect. The 
degree of saturation can only have been slight, however, for the equipment 
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Fic. 5.— Measured points compared with the theoretical diffraction pattern for a point source. 


responded with a deflection not noticeably different from that produced by the 
star on test runs. This will have affected the upper part (shoulder) of the curve 
obtained, but the deviations from the standard form of curve appropriate to the 
passage of a straight edge over a uniform disk still remain at the lower part (toe) 
of the curve. 


Our thanks are due to the Director and Staff of the Radcliffe Observatory 
for the facilities and assistance during the observations, and to the S.A. Council 
for Scientific and Industrial Research, Pretoria, for permission for some of us 
to participate in this investigation. We are also indebted to the other members 
of the Electrotechnology Sub-division of the National Physical Laboratory for 
their invaluable assistance both before and at the time of the observations. 


Royal Observatory, National Physical Laboratory, 


Cape of Good Hope: Pretoria. 
1953 November 30. 


* M.N., 111, 64, 1951. 














MEAN AREAS AND HELIOGRAPHIC LATITUDES OF 
SUNSPOTS IN THE YEAR 1949 


Royal Observatory, Greenwich 


(Communicated by the Astronomer Royal) 


(Received 1953 October 26) 


The following results are in continuation of those given in M.N., 113, 262, 
1953, and are derived from the measurement at Herstmonceux of photographs 
taken at the Royal Observatories of Greenwich and the Cape and at the Kodai- 
kanal Observatory, India. From May 2, the Greenwich photographs were all 
taken at the new site at Herstmonceux Castle. 

Table I gives the mean daily areas of umbrae, whole spots and faculae for 
each synodic rotation of the Sun included in the year 1949; the means for the 
year as a whole are included in Table II, which summarizes the yearly values since 
the last sunspot minimum of 1944. 








TABLE I 
Mean Daily Areas 
Projected* Corrected for Foreshorteningt 
No. Rotation Days - A vy A al 
of | Commenced Photo- e Whole, . Whole 
eu UT. graphed Umbrae Spots Faculae Umbrae Spots Faculae 
1948-49 
1275 Dec. 30°63 27 490 3006 1726 363 2228 2072 
1276 «Jan. 26:97 28 837 4839 2229 592 3485 2678 
1277. Feb. 23:31 27 809 4877 2141 603 3667 2622 
1278 Mar. 22-63 27 422 2355 1793 341 1930 2277 
1279 Apr. 18-91 27 389 2083 2077 280 1520 2504 
1280 May 16°14 27 291 1573 2134 220 1221 2650 
1281 June 12°35 28 382 2070 2435 286 1558 2972 
1282 July 9°55 27 477 2926 2177 343 2107 2646 
1283 Aug. 5°76 a9 500 3014 2107 373 2242 2527 
1284 Sept. 2-00 27 537 3160 2300 392 2288 2736 
1285 Sept. 29:27 28 444 2674 2182 310 1887 2606 
1286 Oct. 26°56 27 407 2653 2242 302 2001 2723 
1287 Nov. 22°86 27 365 2311 2214 273 1739 2665 


* Expressed in millionths of the Sun’s disk. 
+ Expressed in millionths of the Sun’s hemisphere. 


Table III gives for each rotation in the year 1949 the mean daily area of the 
whole spots (corrected for foreshortening) and the mean heliographic latitude 
of the spotted areas for both northern and southern hemispheres. 

The mean heliographic latitudes of the entire spotted area and the mean 
distance from the equator of all spots are also tabulated. ‘The mean values for 
the year 1949 are included in Table IV. 
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TABLE II 
Mean daily areas 

No. of Days Projected* Corrected for Foreshorteningt 

.. Photo- Without Whole Figg te Whole 2 

Year graphed Spots Umbrae Spots Faculae Umbrae Spots Faculae 
1944 366 157 30 160 284 23 126 344 
1945 365 14 102 560 774 78 429 940 
1946 365 ° 389 2381 1794 291 1817 2188 
1947 365 ° 558 3559 2326 405 2637 2894 
1948 366 ° 419 2618 1849 314 1977 2331 
1949 365 ° 485 2873 2140 356 2129 2597 


* Expressed in millionths of the Sun’s disk. 
+ Expressed in millionths of the Sun’s hemisphere. 











Tas_e III 
Spots North Spots South 
of the Equator of the Equator Mean 
= A = ia A ~ Mean Distance 
_— Mean —— Mean Latitude from 
No. Rotation Helio- : Helio- of Entire Equator 
Daily ; Daily i " 
of Commenced graphic graphic Spotted of all 
Rotation U.T. Ares Latitude Area Latitude Area Spots 
1948-49 
1275 Dec. 30°63 1312 19°96 g16 11°86 +6°88 16°63 
1276 Jan. 26-97 1258 14°57 2226 10°87 —1 ‘68 12°21 
1277 Feb. 23:31 2328 13°22 1338 11°Q5 4°03 12°76 
1278 Mar. 22°63 1215 13°28 715 15°33 + 2-69 14°04 
1279 Apr. 18:91 808 13°25 712 14°21 +0°39 13°70 
1280 May 16°14 812 10°98 409 9°93 + 3°98 10°63 
1281 June 12°35 709 10°32 849 10°21 —0-'92 10°31 
1282 July 9°55 768 15°04 1338 15°26 —4°21 15°18 
1283 Aug. 5°76 1702 13°39 540 8=14'II +6°77 13°57 
1284 Sept. 2:00 605 12°39 1683 13°13 —6°38 12°94 
1285 Sept. 29:27 1344 13°46 543 13°06 -+5§ 83 13°34 
{ 1286 Oct. 26°56 1314 16°74 688 8-18 + 8-18 13°80 
1287 Nov. 22°86 1128 16°93 611 9‘90 7-50 14°46 
TABLE IV 
Spots North Spots South 
of the Equator of the Equator Mean 
= A _ P A a Mean Distance 
Mean Mean Mean Mean Latitude from 
Ye Days .1,,  Helio- pings Helio- ofEntire Equator 
| car Photographed _—Daily ~ — hic Spotted _ of all 
! graphec graphic graphic p 
: Area Latitude Area Latitude Area Spots 
: 1944 366 42 19°00 83 22°81 —8-70 21°53 
Old Cycle 7 4°17 6 7°61 —I1'l4 5°72 
j 1944 3664 I. Cycle 35 22°18 77 24°02 —9'63 23°45 
f 1945 365 121° 20°13 309 20°26 -8-92 20°22 
4 1946 = 365 1127 20°74 690 18°79 5°73 20°00 
1947 365 992 ~=—«1658 1645 17°86 —4°91 17°38 
1948 366 936 14°92 1041 13°53 —o:06 14°19 


1949 365 1178 14°30 95st 1214 2°49 13°33 
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Tables II and IV are in continuation of similar tables in Monthly Notices: 
for the years 1874 to 1888, 49, 381, 1889; 1889 to 1902, 63, 465, 1903; IQOT to 
1914, 76, 402, 1916; I913 to 1924, 85, 1007, 1925; 1923 to 1933, 94, 870, 1934; 
and 1933 to 1945, 110, 501, 1950. 

The rotations in Tables I and III are numbered in continuation of 
Carrington’s series. 

The chief features of the record for the year 1949 are as follows :— 

(1) The high values for the mean daily area of sunspots and of faculae— 
respectively 2129 and 2597 millionths of the Sun’s hemisphere—exceed the 
corresponding values for 1948, and indicate a prolongation of maximum 
conditions from the notable peak of 1947. 

Inspection of the run of mean daily areas of sunspots for each synodic rotation 
of 1949, given in Table I, shows very high values for Rotations 1276 and 1277, 
i.e. about 10 months after a similar peak in 1948, which in turn came about 
12 months after the great peak of the cycle in 1947. The high level maintained 
by faculae throughout 1949 is also seen from Table I. 

A short survey of the sunspots for 1949 has been published in Monthly 
Notices, 110, 169, 1950. ‘The largest group of the year, with a mean area for its 
disk passage of 2213 millionths of the Sun’s hemisphere, crossed the central 
meridian on January 22 in latitude 23° N. 

(2) On no day in 1949 were sunspots and faculae absent from the Sun’s 
disk. 

(3) The mean weighted latitude of all spots was 13°-3, ie. 4° nearer the 
equator than for the maximum year 1947 (Table IV). The northern hemisphere 
had the larger mean area (1178 millionths as against 951 millionths for the 
southern hemisphere), the average latitude of the latter being 2° nearer the 
equator than the former. 

(4) The number and distribution, northern and southern hemisphere, of 
spot groups of (a) two days or longer; (5) one day only, are as follows :— 


(a) (5) 

Northern spots 284 69 
Southern spots 229 41 
Total 513 110 


The range in latitude for the longer-lived spots was from 0° to 35°. The latitudes 
of the one-day spots ranged from 0° to 32°. 

An appended table gives the mean daily areas of sunspots and faculae 
(corrected for foreshortening and expressed in millionths of the Sun’s hemisphere) 
for each calendar month of 1949. 


Monthly Mean Daily Areas of Sunspots and Faculae 


1949 
Month Spots Faculae Month Spots Faculae 
Jan. 2240 2206 July 1752 2825 
Feb. 3640 2691 Aug. 2265 2436 
Mar. 3277 2493 Sept. 2185 2806 
Apr. 1882 2453 Oct. 1815 2622 
May 1217 2437 Nov. 2313 2744 
June 1606 2858 Dec 1486 2616 


Royal Greenwich Observatory, 
Herstmonceux Castle, Sussex : 
1953 October 24. 
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THE EFFECT OF WINDS AND BODILY TIDES ON ‘THE 
ANNUAL VARIATION IN THE LENGTH OF DAY* 


Yale Mintz and Walter Munk 


Summary 

The effect of winds is to vary the length of day by 0-7 milliseconds 
between February and August. The bodily tides result largely in a semi- 
annual variation of 0-4 ms. The combined effect of winds and tides is not 
far out of line with astronomical measurements since 1950. The fluid core 
of the Earth hardly participates in these variations in angular velocity and 
there is a half-year slippage of 23 metres between core and mantle at the 
equator. Eventually the astronomical measurements may provide a basis 
for determining the Love number h. 


* The tull text of this paper is published in M.N., Geophys. Suppl., 6, No. 9, 1954- 
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ERRATA 
M.N., 103, No. 3, 1943: 
George Alter, Galactic Absorption and Star Distribution. 


P. 164, line 8 from bottom, 


. 6b a 
for B= loge ar,’ 
1 0b 
read B = loge ae . 
P. 165, line 2, 
3 
for G=e Sloe [@py%o—12(1 —e-i4gfo)], 
3 = 
read G=exp ‘iiss [ape%o— 12(1 ==@ “eo)} . 


M.N., 113, No. 3, 1953: 
Obituary Notices, L. J. Comrie. 
P. 298, 1. 17 from bottom, for Pukekoke read Pukekohe. 


P. 296, 1. 20 from bottom, for ten million read two million. 
M.N., 113, No. 5, 1953: 
R. Wilson, The ionized helium series originating from the fifth quantum level. 


P. 560, line 16, for m=n- (£4/F) read m=n"Z (A/F). 
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